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CORRECTIONS 


Volume 35, 1957 
Page 572. In equation (6), ‘“‘+2(a1:+<a2) { —a;(a;+a2)}~2Ht sin x’’ should read “+2(a:+¢@2) 


{ —a;(a;+1)}~2H¢ sin x”. 


Volume 38, 1960 
Page 229. On page 229, the formula for g (equation (19)) should read 
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Page 332. In line 15 up ‘‘(Wilson, Rose, and Pomerantz 1959)’’ should read ‘‘(Wilson 
Rose, and Pomerantz 1960)”’. 


Page 333. In line 8 up “Wilson, B. G., Rose, D. C., and Pomerantz, M. A. 1959. Can. 
J. Phys. 38, 328.” should read “Wilson, B. G., Rose, D. C., and Pomerantz, M. A 1960. 


Can. J. Phys. 38, 328.” 


Pages 1114-1133. Two errors in the Appendix have been pointed out to the author of the 


paper “Scattering of radio waves by an ionized gas in thermal equilibrium’’. These errors do 
not occur in the corresponding analysis of Pines and Bohm, which is correct within its range 
of validity. 

The first error consists in the omission of the second order term in the expansion of the 
logarithm leading to equation (58). The second occurs in the final integration leading to 
equation (73), where the lack of independence of Fourier components with wave number 
vectors of equal magnitude and opposite direction was not taken into account. 

The corrected forms of equations (6), (7), (73), and (76) are obtained by replacing T in 
the original equations by T/2. The correc ted form of equation (8) is obtained by replacing 
? in the original equation by 2d?. The results obtained for the frequency spectrum of scattered 


power are not affected. 
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TE SURFACE WAVES GUIDED BY A DIELECTRIC-COVERED 
METAL PLANE! 


D. Morris AND A. G. MUNGALL 


ABSTRACT 


Investigations of the phase velocities of TE mode surface waves are described. 
The surface waves were excited over a sand-covered metal plane and the phase 
velocities of the first three TE modes determined as a function of the sand depth, 
at a frequency of about 9300 Mc/sec. A phase comparison system was used for 
the measurements. The simultaneous existence of two modes with different 
velocities, predicted theoretically for certain sand depths, was found experi- 
mentally. 

The variation of relative surface wave power carried above the dielectric layer 
with thickness of the layer is also discussed. 


INTRODUCTION 


A previous paper (Mungall and Morris 1959, subsequently referred to as 1) 
described the determination of the phase velocity of surface waves in the 
first three 7M modes propagated over a plane dielectric-covered metal 
surface. This paper will deal with the determination of the phase velocity 
of transverse electric modes of surface waves propagated over the same 
guiding structure. Most other experimental investigations have been confined 
to the dominant mode of propagation over the dielectric-coated metal plane, 
which is a 7M mode. However, Rich (1955) has given the equations for the 
field components of 7‘E modes and diagrams of the field lines associated with 
the first 7E mode, and the use of this mode has been suggested in connection 
with the design of surface wave antennas (Hansen 1957). In I equations 
relating the phase velocity of the waves to the dielectric depth are given, 
and the relations are shown graphically for both 7M and TE modes for a 
particular value of dielectric constant. 

The measurements described in this paper serve to check the correctness 
of the theory of 7E mode surface wave propagation over a sand-covered 
metal plane. 


Theoretical 

The equations defining TE surface wave propagation over a perfectly 
conducting plane covered with a layer of loss-less dielectric, as given in paper 
I, are as follows: 

‘Manuscript received August 24, 1960. 
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2 7 
(1) v= 42 Ko— a] 


and 
w? 2 } 
(2) B2 cot Bal = 5; (x2—1)—B2 


where v = surface wave phase velocity, 


w = angular frequency, 
ko’ = dielectric constant of layer, 

2 = phase constant of propagation in the dielectric, 
c = velocity of electromagnetic waves in free space, 


l = thickness of dielectric layer. 
The magnetic permeabilities of both the conducting plane and the dielectric 
layer are taken equal to that of free space. 
Equations (1) and (2) can be combined to give a dimensionless relation in 
terms of relative dielectric thickness and velocity, in the form given by 
Hansen (1957) or as follows (Mungall and Morris 1960): 


v a oe "| 
(3) fm icah-iy nmV'—1 


where L = 1/Xo, 
V = o/c, 
Xo = free space wavelength. 


From equation (3) the minimum dielectric depth for 7E mode propagation 
(occurring at V = 1.0) is L = 1/(4(k2’—1)4), ice. corresponding to a cutoff 
wavelength of A, = 4/(xs’—1)}. 

For a given value of V, there is more than one value of Z which will satisfy 
equation (3). This can be interpreted as indicating the existence of higher 
order modes. The nomenclature followed in this paper is to describe the mth 
solution of equation (8) as the 7E, mode. Thus the lowest TE mode is TF. 
The existence of this mode structure was examined experimentally, and the 
results are compared with theoretical predictions in Fig. 5 

It is also of interest to consider the power flow in the air above the dielectric 
for the various modes and also for different velocities within a given mode. For 
TE modes the proportion of the total surface wave power carried above the 
dielectric layer can be shown to be 


(4) P, = ae a (KV — 1) Ee 


P4+P,° ee , ee , 
(x3—1)V 414(4: aE «{-(45%) I 


where P, and Pg, are respectively the power carried in the air and in the 
dielectric. Calculations from (4) for the case x2’ = 3.0, representative of dry 
sand, give the results shown graphically in Fig. 1. The relative power carried 
above the dielectric falls rapidly with decrease of phase velocity, particularly 
for the higher order modes. Experiments showing qualitative agreement with 





this will be described. 
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Fic. 1. Fraction of surface wave power carried above dielectric layer as a function of 
relative phase velocity for first three TE modes propagated over a sand-covered metal plane. 


EXPERIMENTAL 


1. Technique 

The experimental method used was similar to that described in I except 
for the technique of exciting and detecting the surface waves. The method 
depends on the phase comparison of two signals arriving at two mixers (both 
powered by the same local oscillator klystron) via two different paths, but 
both being derived from the same klystron. One path is over the surface 
wave structure and the other consists of a fixed length of waveguide. The 
heterodyne outputs from the two mixers are amplified and applied to an 
oscilloscope, and phase changes are measured by the use of a variable phase 
shifter in the waveguide path. All measurements were made at a frequency 
close to 9300 Mc/sec (Ao = 3.2 cm). 

The simple technique described in I for exciting and detecting 7M mode 
surface waves, i.e. propagation over a metal sheet followed by the surface 
wave structure, could not be used in the present investigation since it is not 
possible to propagate a 7‘E wave close to a metal sheet. Instead, the equip- 
ment shown in Fig. 2 was used. This equipment was contained in a wooden 
box lined with microwave absorbing panels. The signal path was as follows. 
Horn 1 transmitted a wave with its electric vector vertically polarized over 
a metal sheet, path A, at velocity c, to horn 2. The wave received by horn 2 
then travelled along a waveguide twist to horn 3 from which it was radiated 
with horizontal polarization. Part of the energy was launched as a TE surface 
wave and travelled at a velocity v over the dielectric-coated conducting plane, 
path B, to horn 4. The received signal then entered one of the mixers. The 
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Fic. 2. Schematic diagram of equipment for exciting and detecting 7’E mode surface 
waves. 


microwave relay assembly, consisting of horns 2 and 3 and the connecting 
waveguide twist, was mounted on a carriage to facilitate accurate motion in 
the direction of wave propagation. 

The following procedure was then used. The microwave relay was positioned 
initially at one end of its travel and the setting of the phase shifter in the 
waveguide path required to bring the phase pattern on the oscilloscope to a 
straight line was noted. The microwave relay was then moved in 1-cm steps 
to decrease the path over the guiding structure and the corresponding phase 
shifter readings noted. A change of phase took place when the carriage was 
moved if the velocity of the wave propagated over the sand-covered metal 
sheet was less than that of free space. Direct coupling between the trans- 
mitting and receiving horns at the ends of the box was very small, due to 
the difference in polarization of the transmitted and received waves. It was 
reduced to negligible proportions by interposing a sheet of microwave absorbing 
material across the box just above the microwave relay assembly. 

The total path length employed was 200 cm. The surface wave guiding 
structure consisted of a thick aluminum plate, 90 cm long and 60cm wide, 
covered with a layer of dry sand. It was found that, to obtain reproducible 
results, the sand had to be well compacted and possess a smooth surface. A 
removable scraper was devised to satisfy these requirements. The scraper had 
a rounded lower surface, and its carriage ran on ball races on rails bolted to 
the underside of the aluminum plate. The uniformity of the gap between 
the scraper blade and the table as the scraper was moved from one end to 
the other was better than +0.2 mm, after shimming of the rails. The sand 
layers produced by this rounded scraper blade were well compacted and had 
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a smooth surface, and the sand depths were probably uniform to the above 
accuracy. The dry sand had a dielectric constant of 3.0 at the frequency 
used, 9300 Mc/sec. 


2. Results 

Most of the graphs showing phase change as a function of the path length 
over the guiding structure were similar to those obtained for the 7M modes 
and shown in I. With the horns described for exciting the surface waves, an 
accompanying free-space component was always present to a greater or lesser 
extent and the graphs usually showed the results of interference between the 
surface wave and the free-space component. In general, when a surface wave 
predominates over the free-space wave the plot of phase change as a function 
of path length is monotonic, approximating a straight line. The interference 
results in a periodic oscillation of phase about the average slope of the graph. 
When the surface wave interferes with a much stronger free-space component 
the phase changes periodically as the path length is varied. The phase velocity 
of the surface waves can be determined from the average slope of the mono- 
tonic curves and from the periodicity of the oscillatory phase curves. Since 
the proportion of the surface wave power carried in the air above the dielectric 
decreases with decreasing phase velocity (see Fig. 1), if the receiving horn 
was positioned above the dielectric it received a very small proportion of the 
surface wave power mixed with a relatively large amount of the free-space 
power for the low velocity waves. Hence, for better detection, in this case 
the receiving horn was lowered until it rested on the metal table and the 
sand layer was tapered into the mouth of the horn. Even with this technique 
it was not possible to get the amplitude of a low velocity surface wave greater 
than that of the free-space component. However, the mouth of the horn 
extended somewhat above the sand surface and so still received some of the 
free-space component. Also the launching efficiency of the exciting horn was 
probably poor. This horn had to be movable and so it could not be positioned 
below the surface of the sand. Hence, probably only a small proportion of 
power was launched as a surface wave at the lower wave velocities. 

When the sand depth was such that two surface wave modes could propa- 
gate, the presence of both could be inferred from the phase change graphs 
obtained when the receiving horn was first set just above the sand surface, 
and was then set below the surface and the sand tapered into the mouth 
of it. Such graphs are shown in Figs. 3 and 4, which were both obtained for 
the same sand depth, 19.0 mm (ZL = 0.59). Figure 3 shows the phase curve 
obtained when the horn was above the sand surface. The curve shows only 
slight variations due to the interfering free-space wave, showing that the 
amplitude of the latter is much less than that of the surface wave. Hence, the 
amplitude of the interfering free-space wave was probably even less when 
the receiving horn was lowered. From the mean slope of the curve in Fig. 3 
the derived relative phase velocity is 0.92. This is in good agreement with 
the theoretical value of 0.93 for the 7E: mode, computed from equation (3) 
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Fic. 3. Phase change as a function of path length at sand depth 19.0 mm with receiving 
horn above sand surface. TE: mode indicated. 

Fic. 4. Phase change as a function of path length at sand depth 19.0 mm with receiving 
horn just below sand surface. Circles indicate experimental points. Theoretical curve obtained 
for interference between waves in TE, and TE: modes. 


for L = 0.59. The theoretical relative phase velocity of the 7E,; wave at this 
dielectric depth is 0.63. In Fig. 4 the experimental points were obtained with 
the receiving horn below the sand surface. The points are compared with the 
solid curve obtained by calculating the phase change as a function of path 
length when a wave of amplitude unity and velocity 0.93c interferes with a 
wave of amplitude 0.4 and velocity 0.63c. The actual ratio of the amplitudes 
of the waves was not known, but the experimental points tend to follow a 
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curve of the same shape as the theoretical one. Thus, the experimental results 
shown in Fig. 4 indicate the simultaneous propagation of two TE modes, at 
different velocities. 

Figure 5 shows the experimental variation of phase velocity with dielectric 
depth for the first three TE modes compared with the theoretical curves 
computed from equation (3) with x2’ = 3.0. The experimental points denoted 
by circles indicate cases where the phase curves obtained were essentially 
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Fic. 5. Comparison of experimental and theoretical TE surface wave phase velocities as 
a function of sand depth over a conducting plane. Experimental points for A) = 3.2 cm. 


monotonie, i.e. the surface wave component was predominant. The points 
denoted by triangles represent cases where periodic phase curves were obtained, 
indicating that the surface waves were weak. The accuracy of the deter- 
mination of the phase velocity from the periodic curves was usually not 
better than 2-3%. However, the satisfactory agreement between the experi- 
mental points and the theoretical curves appears to be adequate verification 
of the mode structure described theoretically for 7E surface wave propagation 
over dielectric-covered plane metal surfaces. 
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A FABRY-PEROT SPECTROMETER 
FOR AURORAL AND AIRGLOW OBSERVATIONS! 


G. G. SHEPHERD 


ABSTRACT 


A 1-in. Fabry-Perot interferometric spectrometer employing mechanical 
scanning has been constructed to operate at low order. Such an instrument gives 
a resolution similar to that of a grating spectrometer but has, at that resolution, 
a much higher light gathering power. The spectrometer is capable of resolving 
the sodium D lines in the twilight airglow, scanning once per minute. Some 
preliminary results on the variation of D2/D, ratio during twilight are pre- 
sented. This type of spectrometer should be useful for a number of other 
applications. 


1. INTRODUCTION 


The use of the Fabry-Perot interferometer in conjunction with photographic 
spectrographs for high resolution studies is well known. In recent years, it 
has been developed into a scanning photoelectric instrument, mainly by 
Chabbal (1953, 1958a, 6). At the Bellevue Conference on Interferometry* in 
1957, more than thirty papers were given on the use of the Fabry-Perot 
spectrometer, indicating the increasing interest in such a device. In these 
papers, and elsewhere (Jacquinot 1954) it is made clear that the Fabry-Perot 
spectrometer has a much higher light gathering power than has a grating 
spectrometer operating at the same resolution. That is, if G is the light 
gathering power, A the area of the selective element, 2 the solid angle it 
will receive, and 7 the transmission of the instrument, then G = AQT. Fur- 
ther, if R is the resolving power of the instrument, then GR is a constant for 
any spectrometer, but this constant has a much larger value for the Fabry- 
Perot spectrometer than for a grating spectrometer, if both have the same 
area, A. 

For this reason, instruments of this type would be very useful for observa- 
tion of faint sources such as the airglow and transient sources such as the 
aurora. Nevertheless, they have not as yet become very widely used, perhaps 
because they are considered to be complicated technically. Applications of 
interferometry to the airglow and aurora have so far consisted mainly of line 
width measurements at high resolution (Armstrong 1959; Wark 1960). In this 
paper, the design and construction of a relatively simple 1-in. (2.5 cm diameter) 
mechanical scanning Fabry-Perot spectrometer is described. An instrumental 
half width of 2.5 A is obtained at an order of 170. The passbands are 36 A 
apart so that an interference filter can be used as a monochromator. This 
l-in. instrument was intended as a prototype for a larger one, but the results 
of some preliminary observations of the sodium twilight airglow which will 


be described demonstrate that it is itself a useful instrument. 
1Manuscript received August 26, 1960. 
Contribution from the Physics Department, University of Saskatchewan, Saskatoon, Sas- 


katchewan. Supported by National Research Council grant No. A-822. 
*]. phys. radium, 19, 13-434 (1958). 
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2. PRELIMINARY CONSIDERATIONS 


Figure 1(a) shows a schematic arrangement of a Fabry-Perot spectrometer. 
An etalon of spacing ¢ is placed in a gas of refractive index n, followed by a 
lens L, forming a circular fringe pattern at its focus, which is coincident 
with the diaphragm. A small central portion of the fringe pattern passes 
through the diaphragm, is collected by the lens Lz, and falls on a photo- 
multiplier, or some other detector. The transmission of the etalon is a function 


— so —4 





(a) (b) 


_ Fic. 1. (a) Schematic arrangement of a Fabry-Perot spectrometer. (6) Passbands of a 
Fabry-Perot interferometer (Airy function) for a finesse equal to 20. 


of t, m, 1 (the angle of incidence), and o (the wave number of the radiation). 
For fixed ¢, m, and i, the transmission vs. o graph represents the passbands 
of the instrument, shown in Fig. 1(d). The separation between passbands will 
be denoted by Ao, and is equal to 1/2nt. The ratio of this quantity to the half 
width of a single passband, w, is called the finesse, and will be denoted by N. 
The finesse obtained depends on the quality of the surfaces, the reflectivity 
of their coatings, and the size of the diaphragm used. Typical values range 
from 10 to 30. 

The multiple passbands of the Fabry-Perot can cause some difficulty. How- 
ever, if one is studying a simple, isolated spectral feature, a wide passband 
filter will suffice. In other situations where nearby features are present, or 
worse still, a continuum, some form of monochromator is essential. At the 
low orders used here the passbands are 20-40 A apart and an interference 
filter will isolate a single passband. At higher resolutions it is necessary to 
use several etalons in series (Chabbal 19588). 

Two methods of scanning are available. The first is by varying the refractive 
index of the gas surrounding the plates, and the second is by mechanical 
motion of one of the plates. There is a third method, tilting the etalon (varying 
i) but this method has some limitations and will not be considered here. It 
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can easily be shown that using the refractive index method, the scanning range 
is independent of the order, and spacing #, and that for one atmosphere pressure 
variation in air, the scanning range is 1.5 A at 5000 A. For mechanical scan- 
ning, a motion of \/2 causes a scan through one order, which for ¢ = 1 cm 
is only equal to 0.2 A at 5000 A. On the other hand, for ¢ = 0.1 mm, a motion 
of \/2 gives a scan of 20 A. Obviously mechanical scanning is advantageous 
at low order if a large scanning range is desired. 


3. MECHANICAL DESIGN 


A cross-sectional view of the mounting is shown in Fig. 2. The Fabry-Perot 
plates (A) are 1 in. in diameter and } in. thick. They are inexpensive plates, 
intended for use in an undergraduate laboratory, rather than for research 
purposes. The upper plate is held in the upper holder (B) by a three-point 
support. Each point consists of an open vee-shaped piece with the two ends 
touching the etalon plate, and the apex connected through a steel ball to a 
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Fic. 2. Cross-sectional view of the mechanical scanning mounting for the 1-in. Fabry- 
Perot spectrometer. A, Fabry-Perot plates; B, upper (fixed) plate holder; C, lower (movable) 
plate holder; D, air space cover; E, coarse adjustment screw (3); F, fine adjustment screw 
(3); G, invar spacer (3); H, air space; I, annular membrane 0.045 in.; J, air inlet. 


spring-loaded rod, carefully ground into a reamed hole in B. This support 
mechanism holds the upper plate firmly in any desired orientation. The three 
coarse adjustment screws (E) are used to slide the upper plate in its bearing 
pieces to the desired position and are then withdrawn from contact with the 
plate. The lower plate is held in a support which is connected to the lower 
holder (C) through an annular membrane (I), all of which is machined from 
one piece. A cover plate (D) encloses an air space (H), sealed with two ‘‘O” 
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rings. Air pressure produced in H causes elastic deformation of the membrane, 
giving a motion to the lower plate which is linear with pressure. The ‘‘O”’ 
ring connecting the lower plate support to D does not interfere with the 
deformation of the steel membrane. The required thickness of the membrane 
was calculated from an equation given by Timoshenko (1930), which served 
as a guide for the first machining. After initial tests had been made using 
the etalon, the membrane was further machined to 0.045 in. thick, giving a 
displacement of \/2 for a pressure change of 10 cm Hg. This scanning mech- 
anism is a variation of that used earlier by Chabbal and Soulet (1958) in 
which the membrane is deformed by a cam and lever arrangement.* B and 
C are separated by three invar spacers (G), and compression of these by the 
spring and nut devices (F) provides the fine adjustment of the etalon. The 
length of G is chosen to give thermal compensation of the etalon spacing. The 
mounting is machined of SAE 3140 steel, which is a 1.3% nickel alloy, and 
was chosen for its resistance to creep. The apparatus is considered satisfactory, 
being relatively easy to adjust, and stable for periods of several hours or 
more. The scanning motion has been tested over seven orders and the adjust- 
ment found to remain good through that range. 

This mounting is fixed on teflon supports to the bottom of a cylindrical 
cast aluminum alloy cell having 1-in. thick walls. Four Chromalox strip heater 
elements are clamped to the outer wall. The large thermal conductivity of 
the aluminum gives a reasonably uniform temperature over the inner wall. A 
temperature-sensitive resistance, consisting of Balco wire wound on a 0.2 in. 
diameter lucite rod, is inserted in the wall of the casting. A commercial 
regulator designed for domestic application gives adequate temperature 
control. 

The pressure drive for the Fabry-Perot is a simple arrangement, consisting 
of a 2.5 in. diameter sylphon bellows, driven through a cam by a synchronous 
1 r.p.m. motor. The pressure variation produced by the bellows was surprisingly 
linear with displacement. Some non-linearity is produced because the change 
in volume is not negligible compared with the total volume and this was 
allowed for in the cam design. Unfortunately, it was not possible to locate 
the sylphon bellows in the temperature-controlled cell. This did not affect 
the stability unduly, although it was necessary to operate the scanning 
mechanism for about one hour before the etalon was adjusted. 

The detector was an E.M.I. 9502A photomultiplier, cooled with dry ice, 
and was followed by a preamplifier of the type described by Hunten (1953). A 
further triode stage provided attenuation and high frequency cutoff, and 
coupled the preamplifier to a G-10 Varian recorder. 

The interference filter was placed in a box above the etalon, outside the 
cell. (The optical axis of the system is vertical, which is considered to improve 
the etalon stability.) A tilting mechanism is provided for the interference 
filter in order to shift the peak of its passband to the desired wavelength. 

*During the preparation of this manuscript it was learned that Karandikar (1959) has used 


a similar mechanism, employing liquid pressure, not applied directly to the steel membrane, 
but indirectly through another membrane in contact with the plate support itself. 
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The entire apparatus, with the exception of the electronics, has been located 
outdoors on the roof of the Physics building. It performed quite satisfactorily 
here, and runs were taken at temperatures down to —25°C with no great 


difficulty. 


4, OPTICAL CONSIDERATIONS 


When a particular diffraction grating has been selected, the manner in 
which it can be used to study a given spectral feature is quite limited. This 
is definitely not the case for the Fabry-Perot spectrometer, for which several 
parameters can be independently varied. For example, the resolving power 
of a grating is fixed, but that of the Fabry-Perot is variable and can be 
increased without limit. While the variety of adjustments introduces technical 
complications, this flexibility gives the Fabry-Perot several advantages. 

The first problem chosen to test the instrument was the observation of 
sodium emission in the twilight. Because of the low brightness of this emission, 
the highest possible light-gathering-power was desired. Thus the lowest possible 
resolution which would still separate the D lines was to be used, about 3 A. 
Because of the Rayleigh scattered sunlight accompanying the emission, a 
monochromator was essential; but the low order permits the use of an inter- 
ference filter. A filter of 20-A half width and 70% transmission was available* 
and this was considered more useful than a narrower band-pass filter having 
lower transmission. If the Fabry-Perot passband is to be 3 A in width and 
the finesse is NV, the separation between orders is 3N A. This determines the 
etalon spacing that must be used. This also illustrates that even at low 
resolution, a high finesse is desirable. 

To determine the finesse of the etalon, nine-layer dielectric coatings were 
applied and the localized fringes observed in collimated monochromatic 
light. By scanning, the contours of the plates can then be observed as described 
by Chabbal (1958c). Two pairs of 1-in. plates were available and different 
combinations of these were tried, giving a variety of contours. One combina- 
tion was definitely superior to the others, showing a spherical curvature of 
finesse about 20. To investigate the microirregularities of the surfaces a single 
line was scanned in the usual fashion, but using only a small portion of the 
etalon and a very small diaphragm. In this way, all other contributions to 
the line width were negligible, and the observed profile gave the finesse of 
these defects, again about 20 in this case. 

To give GR its maximum value, the finesse associated with the coatings 
should be similar to that of the defects and sc five-layer coatings of ZnS 
and cryolite, with a finesse of about 23, were applied. Some measurements 
were made with three-layer coatings, but the gain in transmission was not 
great enough to warrant the loss of resolution. 

The last parameter to be determined is the size of detector diaphragm, but 
this is readily adjustable and so several can be tried. After some preliminary 
twilight runs this was left at 1.5 A. 


*This was kindly supplied by Dr. D. M. Hunten of the Physics Department, University of 
Saskatchewan. 
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The final adjustments were at an order of 170, giving passbands 36 A apart, 
an instrumental line width of about 2.5 A, and thus a resultant finesse of 
about 14. Under these conditions it was estimated that the interference filter 
would allow the secondary passbands a transmission equal to 7% of the 
primary, and the tertiary passbands 2% each, giving 18% altogether. This 
was considered a satisfactory compromise with the resolution. 








F.P. half width - 3A 
I.F. half width - 20A 
D,/D, = 1.33 


5870 5890 5910 MA) 


Fic. 3. Synthetic sodium twilight spectrum using Fabry-Perot of half width 3 A, and 
interference filter of half width 20 A. Dotted curve: response to white light only; solid curve: 
white light plus sodium emission. 


A synthetic twilight spectrum is shown in Fig. 3, where the white light 
intensity has been arbitrarily chosen equal to that of D» (as seen by a single 
passband without the interference filter) and the D2/D, ratio fixed at 1.33. An 
instrumental passband of 3A was used here. It is seen that the intensity 
drops only to about 50% between orders because of the secondary passbands 
moving in. The sodium emission lines appear to ride on top of the white 
light maximum, making them a little difficult to separate. As will be seen, 
the actual spectra differed somewhat from this. 


5. SODIUM TWILIGHT OBSERVATIONS 

The first observation with the instrument of twilight sodium was made on 
January 29, 1960. Following some subsequent adjustments, 13 runs were 
obtained between February 20 and May 26, 1960. The procedure followed in 
making a run was similar to that of Hunten and Shepherd (1954) except that 
only the zenith was observed. Figure 4 shows some selected 1-minute scans 
from a morning run on February 21. It should be noted that the sodium 
lines were not well centered in the scan. The first is of a low brightness source 
(Shepherd 1954) and should give the white light profile of Fig. 3. It does 
not, in the sense that the peaking effect of the interference filter is not as 
great as was expected, although it is present. The next scan was taken in 
early twilight, with the unrefracted shadow line at 5 km, and shows the two 
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Feb 21,1960, A.M. 


Fic. 4. Selected scans from a sodium twilight run for various heights of unrefracted shadow 
line. Photomultiplier voltages are indicated for each scan. 


Fraunhofer absorption lines. The next scan is taken at 58 km, and shows the 
emission lines quite well resolved. The last is taken with the shadow line at 
120 km, and shows a trace of emission, which must be nightglow. This is the 
only occasion on which the nightglow was observed. 

The procedure used in analysis was also similar to that of Hunten and 
Shepherd (1954) except that the two D lines were treated separately. In 
addition to giving the variation of sodium brightness during twilight these 
results can also give the D./D, ratio, which has not been extensively measured 
so far. A summary of the previous results is given by Chamberlain, Hunten, 
and Mack (1958). In all cases the ratio obtained is an average over one or 
more evenings. Chamberlain, Hunten, and Mack (1958) have calculated the 
D;/D, ratio as a function of sodium abundance (atoms/cm? vertical column) 
for a solar depression of 6.5° and zenith angles of observation of 0° and 75°. 
Donahue and Stull (1959) have calculated D2/D; as a function of solar 
depression for several abundances, but only for a zenith angle of 75°. 

Only for the later runs was the interference filter adjusted to give equal 
transmission at D,; and Ds». For the earlier runs, a correction factor was 
applied by noting the wavelength of the maximum in the low brightness 
source spectrum and using the derived interference filter characteristics. Appli- 
cation of these correction factors gave D2/D, ratios which were much too 
large. The correction was obviously being overestimated and in fact, the 
most reasonable values appeared to be the uncorrected ones. Under these 
circumstances the absolute value of D:/D, for the early runs could not be 
considered reliable and instead it appeared more useful to concentrate in the 
analysis on the variation of D2/D, during twilight. 

For each run, D2,/D, was plotted against Z, the unrefracted shadow height. 
The scatter of plotted points was disturbing at first, and some smoothing 
was attempted, including averaging several runs together. But it was con- 
cluded that these attempts served only to average out the true variations. On 
re-examining the original plots, it was seen that they could be grouped into 
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several characteristic plots. From the calculations of Donahue and Stull (1959) 
one would associate each of these with a different sodium abundance. For an 
abundance of 1 X 10° atoms/cm? column their plot from Z; = 30 km to 80 km 
is essentially horizontal, showing only a slight maximum at the center. At an 
abundance of 5X 10°, the low Z, end of the plot has fallen, making the curve 
unsymmetrical. At 10X 10°, the high Z; end has also fallen, making the curve 
symmetrical again with a central broad maximum. At 20X10°, the high Z, 
end falls even more and for higher abundances, the maximum disappears 
completely, with D2/D, decreasing with increasing Z, through the entire 
range. Eleven runs had variations which could be included by one of these 
cases, but the other two runs did not fit at all, with D2/D, increasing rapidly 
with Z,. The grouping of nights that resulted is given in Table I. For each 


TABLE I 


Grouping of runs according to characteristic D» 








/D, variation, 1960 


Runs Associated abundance 








March 18, May 3, May 5, May 17, May 26 1X 10° atoms/cm? column 


April 18, May 7 5X 10° atoms/cm? column 
February 20 10 X 10° atoms/cm? column 
February 21, March 23, April 11 20 X 10° atoms/cm? column 


April 8, May 6 


group, one run is shown in Fig. 5, along with the calculated curves of Donahue 
and Stull (1959). This is done for abundances of 110°, 510°, 10109, 
and 20X10° atoms/cm? column. Unfortunately, the calculations are for a 
zenith angle of 75° and the runs were made at 0° but the difference only 
becomes important for large abundances. Because of the uncertainty in the 
absolute value of D2/D,, the calculated curves were shifted to obtain best 
visual agreement with the results and this shift is indicated by the arrow 
at the right of each plot. This uncertainty of the absolute D./D, ratio is not 
just that introduced by the instrument (as mentioned above), but includes 
that of the intensities at the bottom of the Fraunhofer lines for D,; and Ds. 
This latter quantity is still not agreed upon to within 10%. Points indicated 
by a cross are those of Chamberlain, Hunten, and Mack (1958) and have 
been interpolated from their Table 3. (The results in this table appear to 
have been shifted from the corresponding plot in their Fig. 3.) Since the 
differences between the shapes of the curves of different calculated abundances 
are slight, it is unfortunate that the experimental scatter is so large. This is 
particularly true for the points at the ends of the curve which have much 
to do with determining the shape; but nevertheless the experimental variations 
are considered to be real. However, the accuracy does not warrant the con- 
clusion, for example, that the abundance on March 23 was 20X 10° and that 
on February 20 was 10X10°. Rather, it merely suggests that the abundance 
on March 23 was greater than on February 20. Support for this is given by 
the fact that in Table I the runs near the seasonal maximum, February and 
March, have been assigned the higher abundances as compared with the runs 
during May. 
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May 6, 1960. P.M. 
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Fic. 5. Variation of Ds/D, ratio during twilight for a number of runs. Solid curve is that 
calculated by Donahue and Stull (1959) for the abundances shown (atoms/cm? column), 
shifted by the amounts indicated. Crosses are points calculated by Chamberlain, Hunten, 
and Mack (1958), unshifted. 


The plot at the bottom of Fig. 5 is for May 6, which like April 8, could 
not be fitted into any group. A smooth curve has been sketched through the 
points. The increase of D2/D, at the end of the run could perhaps be explained 
as a result of nightglow, which has so far not been considered. One would 
have to further assume that the-ratio of twilightglow to nightglow intensity 
was much lower on this evening than for the others for which the effect is 
not seen. The measurements are not adequate to consider this point in detail. 


6. CONCLUSIONS 

A 1-in. mechanical scanning Fabry-Perot spectrometer has been constructed 
and operated at low resolution with a free spectral range (between orders) of 
36 A. By using an interference filter as monochromator, sodium twilight 
spectra have been obtained with an instrumental half width of 2.5 A. This 
can be compared with the 10-A half width for a 4X4 in. grating spectrometer 
operating at roughly the same light gathering power. It is concluded that 
such an instrument is very suitable for studies involving a spectral range of 
about 30A or less. With a scanning monochromator this range could be 
extended considerably. 
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Observations made of the D./D, ratio during twilight gave variations 
which were different from night to night, and which resembled the variations 
for different sodium abundances as calculated by Donahue and Stull (1959). 
Two of the 13 runs were exceptions and did not fit into any of the above 
patterns. Absolute values of D:/D, were not obtained with certainty on all 
runs, but are in general agreement with what is expected. 

Future measurements should have the aim first of eliminating any uncer- 
tainty in the absolute value of D./Dj,, second of reducing the scatter in the 
measurements, and third of making simultaneous absolute brightness measure- 
ments. This will provide two independent abundance values during the same 
run and give a crucial test to the calculations. In particular, it should remove 
the uncertainty in the ratio of residual intensities at the bottom of the Fraun- 
hofer lines. The absolute calibration of a Fabry-Perot is more difficult than 
that for a grating spectrometer because of the uncertain instrumental profile, 
but these difficulties can be overcome. 

While it would be possible to make these improvements on 1-in. instru- 
ments, it has been decided instead to concentrate on a larger 6-cm instrument, 
for which the 1-in. was intended as a test version. It is anticipated that this 
new instrument will produce results an order of magnitude better than the 
present one. 
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ON THE OPTICAL AND ELECTRICAL PROPERTIES OF 
SEMITRANSPARENT COPPER FILMS! 


R. S. ApHAv? 


ABSTRACT 


Simultaneous measurements were made of the optical constants m and k, 
reflectance R, transmittance 7, and of the d-c. conductivity o, of semitransparent 
copper films deposited onto quartz by evaporation in vacuo. All depositions were 


carried out at a pressure p ~ 10-5 mm Hg and with a deposition rate D = 2 
A/sec. The optical constants were determined by an ellipsometer with a photo- 
multiplier tube attachment at the wavelengths 435 my, 546 my, and 630 mug. 
The reflectance values were computed from the optical constants. The trans- 
mittance values were measured at these wavelengths by a unicam spectrophoto- 
meter. The d-c. conductivities were measured by a kohlrausch bridge. These 
measurements were carried out in the thickness range 60 A to 455 A. The dis- 
agreement of the results with the Garnett theory and the Fuchs theory is 
discussed and conclusions regarding the structure and the degree of porosity 
are drawn. 


Sennett and Scott (1950) have used the Garnett theory to explain the 
appearance of a maximum in the absorption ‘A’ against the thickness curve, 
and have shown that the agreement between theory and experiment is perfect 
for silver. Reynolds and Stilwell (1952) have shown that there is a close 
agreement between theoretical and experimental values of the d-c. con- 
ductivity of thin evaporated films above a certain thickness. Givens (1955) 
has stressed the need of pressure p < 10~-> mm Hg during deposition of films 
of copper. In the present investigation, simultaneous measurements are made 
of the optical constants m and k, reflectance R, transmittance 7, and of the 
d-c. conductivity o, of semitransparent copper films deposited onto quartz by 
evaporation in vacuo. The disagreement of the results with the above-men- 
tioned optical and electrical theories is discussed and conclusions regarding 
the structure and the degree of porosity are drawn. 

The semitransparent copper films in the thickness range 60 A to 455 A 
were prepared by vacuum deposition at a pressure p = 10-> mm Hg and at 
a deposition rate D ~ 2 A/sec. The optical constants » and k were deter- 
mined at the three wavelengths \; = 435 mu, A» = 546 my, and A; = 630 mu, 
with the help of an ellipsometer with a photomultiplier tube (RCA 931 A). 
This wavelength range was selected since pure copper is reported to have an 
absorption band with its long wavelength limit between the range 550 mu- 
600 mu. The new formulae obtained by Ditchburn (1955) were used for com- 
puting the real term » and the phase-shift term k of a complex refractive 
index. The reflectance R was computed from the values of m and k. The 
transmittance 7 in the same wavelength range was directly measured by a 
spectrophotometer. The absorption A was obtained by assuming 7+R+ 1 

1Manuscript received March 8, 1960. 


Contribution from the Physics Department, Gujarat College, Ahmedabad, India. 
2Present address: Dominion Observatory, Ottawa, Canada. 
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= 100%. Hence it includes the light scattered by the films, which is of the 
order of 1% over the thickness range and the wavelength range under con- 
sideration. The d-c. conductivity was measured by means of a Kohlrausch 
bridge. The method used in measuring the film thickness was essentially that 
described by Tolansky and his collaborators (1948). 

The Garnett theory assumes that (a) the film may be represented by 
spherical metal crystallites of diameter that is small as compared to the 
film dimensions; (b) these spherical particles are distributed at random in 
space and the density of the film is characterized by the volume factor g 
giving the fractional volume occupied by the metal. A film of metal of optical 
constants and k (= nk) exhibits apparent optical constants m, and k; 
(= m,k,) related to n, k, and q by the formulae, 


nik, = 3qb/[(1—ga)*+4q°b"] 


a 
[(1—ga)"+4q°b"] 


where g is the volume factor and naturally 0 < g < 1. The constants a and 


° ° 
kin, = 2 


b are given by 


(k°—n*+1)(k*—n*—2)+4n*k® 


= (F —n?—2)°+4k 


’ 


b = 3nk/[(k°—n?—2)?+4k'n’]. 
The product mk; gives a measure of absorption of the film. According to 
Sennett and Scott (1950), the curve ,k; against q possesses a maximum 
if A > 0 where 
K = 2(k®—n?)'—9(k?—n?)?+ (12+8n?k?) (k? —n?) —4(3n2k? +1). 


For the wavelengths used, the values of A are greater than zero and are 


shown in Table I. 


TABLE I 
Theoretical values of n, k, a, 6, and K for the three wavelengths 
\ = mu n k a b RK 
435 1.15 1.70 0.92 0.38 2.45 
546 1.07 2.25 1.21 0.27 137.2 
630 0.59 5.78 1.09 0.01 6.510 


Hence the theory expects a maxima for the absorption against the thickness 
curve for the three wavelengths under consideration. Table II shows the 
theoretical values of 2,k; assuming that g ranges from 0.6 to 1. 

The experimental results of transmittance 7), reflectance R, and absorption 
A for these copper films at the three wavelengths, are in the Table III. The 
absorption 4 is negative for a film of thickness 60 A for wavelengths \. = 546 
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TABLE II 
Theoretical values of mk, at the three wavelengths for various 
q values 
Volume nik, 
factor ee 
q A, = 485 mp Ac = 546 my Az = 630 mu 

0.6 1.66 2.80 0.164 
0.7 1.93 5.6: 0.375 
0.8 2.08 3.47 1.40 
0.9 2.06 3.12 33.9 
1.0 1.95 2.45 25.0 





TABLE III 


Transmittance 7, reflectance R, and absorption A for thickness range 60 A to 455 A, for the 
three wavelengths 

















Thick- Ai = 435 my Ae = 546 mu A; = 630 mu 
ness — —- —- <a — - - 
in / T% . R% : A% T% R% AG T% R% A%N% 
60 67.0 32.0 1.0 75.0 25.8 —-0.8 71.0 31.3 —2.3 
105 47.0 29.1 23.9 55.0 ea 17.3 52.0 28.3 19.7 
165 34.2 32.4 33.4 42.0 31.5 26.5 38.0 32.7 29.3 
212 24.8 39.5 35.7 31.6 33.0 35.4 26.0 34.8 39.2 
286 14.9 42.2 42.9 21.0 37.6 41.4 15.7 43 .6 40.7 
350 11.4 44.0 44.6 16.2 39.1 44.7 11.0 45.8 43.2 
455 6.5 47.6 45.9 10.8 41.5 47.7 6.7 56.2 37.1 





mu, and \; = 630 mu, because below the depth of penetration (of about 150 A) 
the substrate partly modifies the reflection and also because the film structure 
becomes more granular. 

Fuchs (1938) made an analysis of the problem of electrical conductivity of 
thin metallic films with the help of the electron theory of metals. He obtained 
expressions of the form, 


a/oo = 4/3M log(1/M) for thin films, 


where o = electrical conductivity of film, 

oo = electrical conductivity of bulk metal, 

M =t/I, 
where ¢ is the thickness and / is the mean free path of electrons in the metal 
at a given temperature. 
In this case, l 


I 


420 A, 


ll 


0 


0.5618 X 10° (ohm-cm)~!. 


Dingle (1950) does not assume that there is a diffuse scattering at the surface 
of the film, and modifies the formulae by multiplying it by the ratio (1—p)/ 
(1+) where p is the probability of an electron being scattered elastically 
at the surface. 

The experimental and theoretical results of electrical conductivity of the 
same films are shown in Table IV. 
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TABLE IV 


Experimental and theoretical values of d-c. conductivity of the same films in the thickness 
range 60 A to 455A 











Thick- Resistance in: Experimental Theoretical 

ness Vacuum Air —_—_ 
inA (ohms) (ohms) k=t/l 1+log (o0/c) k=t/l  1+log (¢0/c) 

60 © © — — — — 

105 10,000 10,000 0.25 4.74 0.01 2.42 

165 7.00 7.90 0.39 1.78 0.10 1.67 

212 4.76 5.60 0.50 1.72 0.20 1.47 

286 3.01 3.25 0.68 1.66 0.50 1.27 

350 2.39 2.58 0.88 1.64 1.00 1.16 
455 1.81 1.89 1.08 1.63 2.00 1.08 





The theoretical values of mk, against volume factor g as given in Table II 
are plotted in Fig. 1. 

The experimental values of absorption A against thickness as given in 
Table III are plotted in Fig. 2. 

Experimentally, therefore, absorption maxima occur for wavelengths above 
\ = 550 mu, while theory suggests maxima above \ = 435 my. The explana- 
tion of the differences is to be found in the intrinsic differences in structures 
of experimental and theoretical films. The results of Sennett and Scott (1950) 
with silver show that the degree of fit with the theory is better for films 
deposited at higher rates of evaporation. The results of Krautkramer (1938) 
show that the degree of fit is better for films deposited at a higher substrate 
temperature. In the present investigation the rate of deposition is about 


° 


2 A/sec, for all the films and the substrate temperature is about 30° C. Halli- 
day, Rymer, and Wright (1954) have examined thin films of copper obtained 
by deposition im vacuo, by an electron diffraction method. Their results show 
that thin copper films consist of crystallites in the form of disks, of diameter 
90 A and thickness 15 A. Electron micrograph of copper films obtained by 
Sennett and Scott (1950) shows that the aggregation of copper crystallites 
does not appear as regular and as well defined as in the case of silver. 

Figure 2 indicates that the maxima occur at about 320 A thickness for a 
wavelength of A; = 630 mu. Figure 1 indicates that the maxima occur for 
the same wavelength at a g value of 0.9. Hence the porosity of a film of 
thickness 320 A is of the order of 10%. Again, the absorption A falls to zero 
at a thickness of about 60 A for the wavelength under consideration. The 
indications are that a film of copper of 60 A in thickness is almost trans- 
parent. This could be on account of an oxide film of 60 A in thickness formed 
on all the copper films immediately after exposure to air. Garforth (1948) 
observed that an oxide layer of 70 A was formed on the copper films even at 
a pressure of oxygen of 10-? mm Hg in about 7 hours at room temperature. 

The above observations were further carried out over a period of 1 day. A 
typical result at A; = 630 my, over a period of 18 hours is shown by a dotted 
line (Fig. 2). As the oxidation proceeds, the absorption maxima become less 
marked and move towards shorter thicknesses. The absorpticn also falls to 
zero at a greater thickness, giving a higher intercept on the thickness axis. 
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Experimental curves of absorption A against thickness ¢ in A. 
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The initial 60 A intercept moved to 85 A over a period of 18 hours, indicating 
that an oxide film of 85 A is formed over these films. Hence the absorption 
against thickness curve provides a method of estimating the thickness of 
overlying oxide film. 

In Fig. 3, the experimental and theoretical curves of (1+log(eo/c)) do not 
coincide. Above a certain thickness, they run almost parallel indicating the 
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Fic. 3. (1+log(oo/o)) against thickness in A. 


porous structure of experimental films, even for thicknesses greater than the 
mean free path of electrons in metal. Hence any formulae for d-c. conductivity 
of thin films must take into account the volume factor g. In the present 
investigation, the critical thickness at which electrical conductivity is almost 
established is found to be about 120 A. This is almost double that for silver, 
as reported by Oppenheim (1956). This again depends upon the crystalline 
property of bulk metal and upon the numerous conditions of preparation of 


the film. 


The author is very much indebted to Dr. Y. G. Naik for suggesting the 
problem and to Mr. D. K. Gijare for helping in the techniques of preparation 
of thin films. Thanks are also due to Dr. D. M. Dennison for many valuable 


suggestions and encouragement. 
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THE DECAY OF Pm"! 


AGpDA ARTNA AND MARGARET E. LAw 


ABSTRACT 


The 52.8-hour activity of Pm"® has been investigated using a high resolution 
beta spectrometer, a lens type coincidence spectrometer, and a scintillation 
spectrometer in conjunction with a multichannel analyzer. The beta spectrum 
was found to consist of two groups with maximum energies of 1.072+0.002 Mev 
and 0.786+0.004 Mev, and intensities of 97.1+-0.4% and 2.9+0.4% respectively. 
A gamma ray of energy 285.7+0.3 kev was found to be in coincidence with the 
0.786-Mev beta group. No other gamma rays with intensities greater than 0.1% 
were found. The K conversion coefficient for the 286-kev transition was measured 
to be 0.075+0.008. This together with the values of 6.50.7 and 4+1 obtained 
for the K/L and L/M conversion ratios respectively indicate that this transition 
is M1 in character with less than 10% E2 admixture. 


INTRODUCTION 


The 52-hour activity of Pm!** has been investigated previously by Rutledge, 
Cork, and Burson (1952), and by Kondaiah (1952). Rutledge et al. reported 
a strong 285-kev gamma ray in coincidence with a 1.05-Mev beta group, and 
a weak 1300-kev gamma ray. However, Kondaiah found that the decay con- 
sisted of a single beta group of end point 1.05 Mev leading to the ground 
state of Sm!4°. More recently, Chapman, Grace, Gregory, and Sowter (1960) 
have carried out nuclear alignment studies on Pm!4*. They also detected a 
285-kev gamma ray associated with the 52-hour activity, but found it to be 
in coincidence with a weak beta group of end point 0.770 Mev. Since there 
is obvious disagreement between these three reports it was decided to investi- 
gate this decay scheme further. 

Pm?*° is the product of the 2-hour negatron decay of Nd!4°. Neodymium 
oxide enriched to 82% in Nd'48 was irradiated for 15 hours in the high flux 
reactor at Oak Ridge, and was received some 30 hours after removal from 
the reactor. The sample also contained 4% Nd!4* and 7% Nd", producing 
the 11-day Nd!47 and 27-hour Pm!®*! activities. In addition, a small amount 
of samarium was present, giving the 47-hour Sm!" activity. 


EXPERIMENTAL 

A. Gamma-ray Spectrum 

The gamma-ray spectrum of the sample was investigated using a scintilla- 
tion spectrometer consisting of a DuMont 6363 photomultiplier and a 3-in. 
Nal crystal, in conjunction with a RCL 256 channel analyzer. Figure 1 shows 
a typical spectrum taken 7 days after removal of the sample from the reactor. 
No higher energy radiations were seen with intensities greater than 10% of 
the 540-kev gamma ray. 


1Manuscript received August 8, 1960. ; 
Contribution from the Department of Physics, McMaster University, Hamilton, Ontario. 
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COUNTS PER SECOND (xi0%) 





CHANNEL NO 


Fic. 1. Gamma-ray spectrum taken 7 days after removal of the sample from the reactor. 
The energies of the peaks are given in kev. 


The spectrum was followed over a period of 20 days to obtain the half 
lives of the various peaks. The majority of these could be accounted for by 
the known gamma rays resulting from either Nd'47, Pm'!, or Sm** decay. 
The 286-kev gamma ray was found to have a half life of 54+2 hours and is 
therefore associated with the Pm!'** decay. No other gamma rays were found 
to decay with this half life. 


B. The 286-kev Transition 

(1) Coincidence Experiments 

The beta spectrum in coincidence with the 286-kev gamma ray was investi- 
gated with a double long-lens coincidence spectrometer of the Gerholm type 
(Gerholm 1955). One spectrometer (No. 1) was set on the 286-kev K con- 
version line and the beta spectrum scanned with the other (No. 2). A Bell, 
Graham, and Petch (1952) fast-slow coincidence circuit was used, set at 
8 musec resolving time. The sources were prepared by dissolving the neo- 
dymium oxide in nitric acid, and evaporating to dryness a small droplet of the 
solution on a backing of V.Y.N.S. film. The resulting sources were approxi- 
mately 2mm in diameter and 300 ug/cm? thick. 

In all, three complete experiments were done. Each experiment lasted 4 
days, during which time the beta spectrum was scanned approximately 20 
times. Figure 2 shows the Fermi plot of the coincidence spectrum obtained 
from one of the experiments. The mean end point of the beta spectrum in 
coincidence with the 286-kev K line was found to be 0.788+0.009 Mev. 

From the singles counting rate in spectrometer No. 1, set on the 286-kev 
K peak, it was possible to obtain an accurate determination of the half life 
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Fic. 2. Fermi plot of the beta spectrum in coincidence with the 286-kev K conversion 

line. 


of the Pm!** activity. The counting rate obtained had to be corrected for the 
27-hour Pm"! and the 11-day Nd!47 contributions. The latter was measured 
after all the shorter-lived activities had died, while the former was estimated 
from the intensity of the gamma rays associated with Pm!*' as measured 
with the scintillation spectrometer. Neither correction amounted to more than 
5% of the total. In addition, it was estimated from the internal conversion 
spectrum that 1.5+1.0% of the activity was due to 47-hour Sm'**. However, 
because of the similar half lives of Sm'5* and Pm'*’, the half-life plot was 
insensitive to this small amount of impurity, and over a period of 4 days 
the correction for this was negligible compared to the statistical errors. The 
resulting half life of Pm'4* is 52.8+0.3 hours. 

(ai) External and Internal Conversion Experiments 

The internal and external conversion lines of the 286-kev transition were 
studied using a 50-cm Siegbahn type beta-ray spectrometer (Johns ef ai. 
1953). Line sources for this instrumenr ¢ prepared as described above 
(B (i)), the source backing in this case being aluminum-coated mylar. These 
sources were 2.5 cm long and had varying widths and thicknesses depending 
on the resolution and counting rates desired. 

The A internal conversion coefficient, ax, was measured directly by com- 
paring the number of A conversion electrons with the gamma-ray intensity 
from the same source. The number of gamma rays was obtained from the 
external conversion peak according to the equation V, = Vy...XfXk where 
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Ny.e. is the photoelectron peak height, f is the photoelectric yield for a given 
gamma-ray energy and for a given radiator material and thickness (this 
factor has been discussed by Johns et al. 1954), and k is a factor depending 
on the source radiator geometry. ax is then given by 


ap = Ner/N, _ Nex/Np.e.Xf Xk. 


The source radiator geometry factor, k, was measured by repeating the 
experiment in the same geometry with the 412-kev Au’ transition, for which 
ax = 0.028 (Wapstra et al. 1958). Figure 3 shows the internal and external 
conversion peaks for both these transitions. The gold source consisted of 
0.52.5 cm? rectangle of gold-coated mylar (30 ug/cm? Au). The promethium 
source of the same dimensions had an approximate thickness of 500 ug/cm?. 
The continuum underlying the external conversion peaks is mainly due to 


100 
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2050 2100 450 1500 
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Fic. 3. The internal (A) and external (B) conversion lines of the 412-kev gamma ray of 
Au'® (i) and the 286-kev gamma ray of Pm" (ii). 
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high energy betas which have penetrated the radiator backing used as the 
beta stopper. In the case of the promethium peak, some of the continuum 
is also due to Compton electrons accompanying the higher energy gamma 
rays of Pm!'*!. In order to check the reproducibility of the geometry, the 
experiment was carried out with two different radiators of 2.2 mg/cm? and 
2.9 mg/cm? uranium. The K conversion coefficient was found to be 0.075 
+0.008, the error being mainly due to the uncertainty in the 286-kev external 
conversion peak height. The photoelectric yield factor for this instrument, f, 
is known to better than +5%. 

The internal conversion lines were studied using a 0.2 X2.5 cm? beta source, 
approximately 50 ug/cm? thick. The spectrometer resolution under these 
conditions was 0.5%. Figure 4 shows the K, L, and M peaks. The K/L and 
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_ Fic. 4. Internal conversion peaks of the 286-kev transition. Superimposed on the Z peaks 
is the normalized outline of the K peak. 


L/M ratios obtained from these were 6.5+0.7 and 4.0+1.0 respectively. In 
addition, from a comparison of the profiles of the A and ZL lines it could be 
concluded that the Z peak consisted mainly of the Z; component. An upper 
limit of 0.1 could be set for the Z,+Z2/Z; ratio. 

The energy of the transition as measured in both external and internal 
conversion experiments is 285.7 +0.3 kev. The standards used for this measure- 
ment were the thorium F line and the 411.77-kev gamma ray of Au! for 
internal conversion, and the 316.46 line of Ir'®?, and the 411.77 line of Au'®’ 
for external conversion. 
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C. Beta-ray Spectrum 

The beta continuum was scanned with the high resolution spectrometer. 
The end point of the spectrum as obtained from a Fermi plot is 1.073+0.002 
Mev. Subtracting from this the energy of the gamma ray, the end point of 
the second group is 0.786+0.002 Mev. Since only one gamma ray was found 
to belong to this decay the beta spectrum was assumed to consist of only 
these two groups. 

The method of Fermi analysis normally used to obtain branching ratios 
could not be used here because of the presence of impurities, in particular 
Sm!53, which have end points similar to that of the inner group. This type 
of analysis is very sensitive to such small corrections. Instead, the intensity 
of the 0.786-Mev beta group (6;) was calculated in terms of the number of 
conversion electrons (x), and either the number of betas in the ground 
state transition (.Vgo), or the total number of betas (.Vg7), according to the 
equations 

by = Nex(1+ar)/axNpot+axNex(1+ar), (i) 


b, = Nex(1 +ar)/axNer. (ii) 


Figure 5 shows a Fermi plot of the total spectrum. The number of ground 
state betas could be found by reconstructing the spectrum as shown in the 
insert (curve B). The second curve in the insert (curve A) shows the total 
beta spectrum after the Nd!47, Pm'®!, and Sm!5* contributions had been 


1000 2000 8000 
Hp 
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, 


200 300 400 500 600 700 600 900 1000 100 
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Fic. 5. Fermi plot of the total beta spectrum. The insert shows the total beta spectrum 
(A), and the spectrum of the high energy beta group (B) reconstructed from the straight 
line of the Fermi plot (C). 
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subtracted. Since curve B is an ideal spectrum with no source thickness 
effects, equation (i) gives an upper limit for 5;. The spectrum was scanned 
three times. The mean value of 5; obtained from equation (i) was 0.033 +0.004, 
and from equation (ii) 0.028+0.003. A similar analysis of the singles spectrum 
measured in spectrometer No. 2 of the coincidence spectrometer gave results 
in agreement with these. 

In addition, the relative intensity 5; could also be calculated from the 
data of the coincidence experiments, using the following equation: 


by =A coine(2)/A totai(2) X (1 tar) /ax . w(1) 


where 4 coinec2) iS the area under the beta continuum in coincidence with the 
286-kev K peak; A tota2) is the area under the total beta continuum, and 
w(1) is the transmission of spectrometer No. 1 set on the 286-kev K peak. 
Using values of a7 = 0.089+0.009, ax = 0.075+0.008, and w(1) = 0.017 
+0.002, a mean value of 0.029+0.006 was obtained for d). 

From these results the intensities of the 1.072-Mev and 0.786-Mev beta 
groups were found to be 97.1+0.4% and 2.9+0.4% respectively, leading to 
log ft values of 7.1 and 8.1. 


SUMMARY 
Table I compares the measured values of the conversion coefficients and 
ratios for the 286-kev transition with the theoretical predictions for several 
multipolarities. From these results it was deduced that the gamma ray is 
>90% M1 with possible E2 admixture. 


TABLE I 
Conversion data for the 286-kev transition 


Theoretical conversion coefficients and ratios (k = 0.56, Z = 62) 
. ene ———-——-— Experimental 








El E2 E3 M1 M2 values 
aK 0.014 0.050 0.16 0.078 0.33 0.075+0.008 
QL 0.0018 0.011 0.084 0.010 0.058 0.011+40.002 
K:L ae 4.4 23 7.9 §.7 6.5+0.7 
Ly:Le:L3 14:1:1 23 -6:1 1:2.5:1.5 74:1 7:h.7¢31 LL, >L24+L; 
2.8 2.8 2.8 3.0 4+i 


L:M 3.0 


Figure 6 shows the proposed decay scheme for Pm'**, which is essentially 
in agreement with that of Chapman ef al. (1960). The only discrepancy is 
in the branching ratio of the beta groups. The value of 1.8+0.3% obtained 
by Chapman et al. for the relative intensity of the inner group, is lower than 
the value measured in this work. The spin of the ground state of Sm'*® is 
known to be 7/2 (Bogle and Scovil 1952). From the M1 character of the 
286-kev transition it can be deduced that the excited state has negative 
parity with spin 5/2, 7/2, or 9/2. Chapman et al. (1960) have ruled out the 
7/2 possibility from their alignment studies. There is no evidence from 
Coulomb excitation experiments (Heydenburg and Temmer 1955) for collective 
motion in this nucleus. It would therefore seem most likely that the 286-kev 
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B- 0.786 MEV 
2.9 % 
Log ft 8.1 
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97.1 % 
Log ft 7 0.2857 MEV $/2-, 9/2~ 


M,>90 % 
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Fic. 6. The decay scheme of Pm", 


level is due to particle excitation with spin 5/2-. The log ft values and the 
allowed shapes of the two beta spectra (Figs. 2 and 5) indicate that these 
are both first forbidden non-unique. This leads to a positive parity for the 
ground state of Pm'**. The possible spin assignments are 5/2, and 7/2 if the 
spin of the 286-kev level is 5/2, and 7/2, or 9/2 if the excited level has spin 
9/2. The shell model prediction for the spin of Pm?*® is 5/2. However, recent 
unpublished work at Berkeley (quoted by Chapman et al. 1960) has shown 
that the spin is 7/2. 

In comparing this decay scheme with those of neighboring odd neutron 
nuclei one might expect a level at approximately 100 kev. Sm!4755, Gd?5!37, 
and Sm!5%g9 have levels at 121 kev, 108 kev, and 66 kev respectively. Since 
no evidence was found for a gamma ray of this energy with intensity comparable 
to that of the 286-kev transition, this level, if it exists in Sm?4%, must be 
very weakly fed. 
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PROPAGATION OF ELECTROMAGNETIC WAVES ALONG A 
THIN PLASMA SHEET! 


JAMEs R. WAIT 


ABSTRACT 

It is shown that a thin ionized sheet will support a trapped surface wave. The 
effect of a constant and uniform magnetic field is to modify the phase velocity 
and polarization of the surface wave. The essential features are illustrated by 
numerical results for selected values of the electron density, collision frequency, 
and gyro frequency. The effect of locating the plasma sheet near and parallel to 
a conducting plane is also considered. In this situation other modes of a waveguide 
type are possible in addition to the surface wave. 


INTRODUCTION 


The propagation of electromagnetic waves in a plasma medium is currently 
receiving a great deal of attention. In the theoretical approach to such prob- 
lems, it is usually assumed that the medium is infinite in extent and homo- 
geneous in its macroscopic electrical properties. The complications which 
arise when the medium is inhomogeneous are considerable. Nevertheless, it 
is desirable to consider the influence of boundaries if experimental data are 
to be interpreted. In particular, such boundaries or interfaces can support 
surface waves which propagate in a manner very different from that in homo- 
geneous media. It is the purpose of this paper to consider the surface waves 
which may exist in a thin ionized (plasma) sheet in the presence of a steady 
magnetic field. Following Poeverlein (1958), the thickness, 5, of the sheet is 
taken to be very small so that the currents, J, flowing in the sheet are purely 
transverse. 

FORMULATION OF BOUNDARY CONDITIONS 


Starting with the relevant equations (Ratcliffe 1959) for an electron plasma 
and harmonic time dependence, the (boundary) conditions relating the fields 
just above and just below the sheet can readily be obtained. If the sheet is 
contained in the plane z = 0, the equations of motion may be written 


(v+iw)J,—wJy = wore, 
wJ2+(vt+iw)J, 


where v is the (constant) collisional frequency of the electrons with the ions, 
wo” is the electron plasma frequency, w, is the z component of the (electron) 
gyro frequency, and e is the dielectric constant of free space. Using a sub- 
script 1 to denote the region above the sheet and a subscript 2 to denote the 
region below the sheet, the boundary conditions 


Ay,— Hey, = — Js } : ea 
A,,—H2, = J, for 2. = 0 


1Manuscript received August 22, 1960. 
Contribution from the National Bureau of Standards, Boulder, Colorado, U.S.A. and the 
Laboratory for Electromagnetic Theory, Technical University, Copenhagen, Denmark. 
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immediately follow from Ampere’s law. To within the same approximation 


E = FE, = Es, 

3 7 - ~ |} for s = 0. 
( ) i 

Using equations (1) and (2), the boundary conditions may be written in the 
form 


nlHiy—H2,] = -—ME,—NE, 





(4) 
nl iz — Ho-z] = — NE,+ ME, 
where 
2 

. — 600 (vtiw) 

(5) M = (v+iw)” en? 
and 

= €w9 bw, : 
(6) N — (v+iw)’+w, 1, 


and » = (u/e)? = 120m. 

The boundary conditions as developed above must be met by any form 
of the electromagnetic waves which can exist in the presence of the sheet. 
The restrictions are that the sheet thickness is infinitesimally small. As a 
consequence, only the component of the steady magnetic field perpendicular 
to the sheet is important. This is the reason the transverse components of 
the gyro frequencies do not appear in the expressions for M and N., 


THE SURFACE WAVE SOLUTION 

We shall now show that such a sheet can support a surface wave. Whether 
such a wave can be excited by a physical source has been discussed elsewhere 
(Wait 1960). The properties of such a wave will be discussed here without 
consideration of the source. Furthermore, the wave is assumed to propagate 
in the x direction parallel to the sheet. The variation in the y direction is 
taken to be zero, so that the disturbance is purely two-dimensional in nature. 

To facilitate the solution of the problem, electric and magnetic Hertz 
vectors are introduced. As will be shown, it is only necessary that these 
have z components which are denoted II and II*, respectively. Sometimes 
these are called Hertz potentials. The field components are expressed in terms 
of these two scalar functions by 





3 Oo 
(7a) E, = as II 

; . on” 
(7b) E, = thors 


by 
N 
ll 
—s, 
=~ 
+ 
la, 
RIS, 
‘niga 
ms 
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(7c) 
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(82) H, =< 2a" 
7 @xdz 
. Ol 
(8d) H, = — ter 
( _) 
(8c) H,=\k +o2 II. 


The boundary conditions at the sheet (i.e., z = 0) are then given by 


oll; 








(9a) ik(II,—T2) = Mo +Niknlhi 

am* a Pe 
— ae ee ed 
(9b) | oe N+ Miknlhi 
(9c) Ii = lp 

aml, dalle 

Cc St nes, gaat hs 
(9d) Oz Oz 


We now assume the following form of the solutions 


(10a) Il, = ae~*°%e~ #8 
(10d) Ie = a’et *O%e— #2 
(10c) nllt = pe *C%g Sz 
(10d) nls = b’et #C%e— B82, 


Since these are to be solutions of the wave equation 
C+ = 1. 


Furthermore, if the fields are to be positively damped for directions away 
from the sheet 


Im C <2. 


The latter two boundary equations, (9c) and (9d), immediately indicate 
that 


a’ = —a and b’ = Db. 


The first two boundary equations, (9a) and (96), then give 





6 _ MC+2 
(11a) re 
and 
(110) © ieee 


a” M420 
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In order that the wave be a solution to our problem, the right-hand side of 
these two equations should be equal. This leads to the following equation 


for C 
(12) (M+2C)(MC+2)+N°C = 0. 
Solutions are 


—(V? 2 9 2 2 of 
(13) C= Cy = AEN eM EN FA) 16M*} 


Normally only one of these solutions satisfies the condition Im C < 0. This 
root is denoted Co; it corresponds to the trapped surface wave which propa- 
gates in the direction of the positive x-axis. The corresponding value of S 
is given by ; 


(14) So = (1—C,°)}. 


PROPERTIES OF THE SURFACE WAVE 
To present the results in a graphical form, certain dimensionless parameters 


are now introduced. They are 


202. ont tet 
cw w w 
and therefore 
‘ 1-—71Z 
M = iB 
(15) ‘ aoa ae 
: ; QO 
°) Ws Begs 
BO -a-izy 


Since the surface-wave fields vary as exp(—ikSox), the phase velocity 
parallel to the sheet is 1/Re Sy relative to the velocity of light c. This quantity 
is plotted in Fig. 1 as a function of B which is proportional to the thickness 
of the sheet. Values of Q, which are proportional to the gyro frequency, are 
shown on the curve. The case Q = 0 corresponds, of course, to zero gyro 
frequency or no magnetic field. Two values of Z are shown. Z = 0 corre- 
sponds to a vanishing collision frequency. The influence of collisions, shown 
by the curves for Z = 0.2, is not great. In general it is seen that the phase 
velocity is less than the velocity of light. It is also seen that an interesting 
transition region occurs near Q = 1. 

The attenuation of the waves in the x direction is —Im kS» in nepers 
per meter. Thus —Im Sy is a measure of the damping of the wave. If Z = 0 
this quantity is zero since Sp is a real quantity. In Fig. 2 this “horizontal 
attenuation” is plotted as a function of B for three values of Z. As expected, 
the damping of the surface wave increases as Z becomes larger. Again there 
is a transition effect for Q near 1. In fact, as Q becomes greater than 1 the 
horizontal attenuation becomes very small even when collisions are con- 


sidered. 
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Fic. 1. The phase velocity of the surface wave relative to c, as a function of the normalized 
thickness of the thin plasma sheet. 
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Fic. 2. The horizontal attenuation of the surface wave as a function of the normalized 


thickness of the thin plasma sheet. 
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Some further light is shed on the problem when the vertical attenuation 
is considered. Since the fields vary in the vertical direction according to 
exp(—ikC\z|) it is seen that a measure of the decay of the surface wave is 
—Im Cy. This quantity is plotted in Fig. 3 as a function of B. It appears that 





VERTICAL ATTENUATION, —ImCy 
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Fic. 3. The vertical attenuation of the surface wave as a function of the normalized thick- 
ness of the thin plasma sheet. 


the wave is highly trapped when Q is small. As Q increases to values some- 
what greater than 1, the wave spreads out on both sides of the sheet. In 
fact, in the limiting case of Q@—, the surface wave degenerates into a 
plane wave travelling along the sheet. In this situation the sheet is effectively 
a perfect magnetic conductor. 

The wave polarization, P, is conveniently defined as the ratio of the com- 
ponent of the electric field parallel to the sheet to the component perpendicular 


to the sheet. That is 


Ey _ Ey a MCo+2. 


ie ie ae 
This quantity is plotted in Fig. 4 as a function of B for the case Z = 0. The 
curves for small finite values of Z would not be very different from those 
shown. It may be observed that for small values of Q the field is almost 
vertically polarized. When Q becomes greater than 1 there is a tendency for 
the polarization to change to horizontal. This change in polarization is quite 
abrupt for small values of B as indicated in Fig. 4. 

These results could indicate that surface waves propagating along air- 
plasma interfaces may be profoundly influenced by the presence of an impressed 





(17) F 
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Fic. 4. The polarization of the surface wave in a plane normal to the direction of propaga- 
tion. 


magnetic field. In fact, this behavior may find application in microwave 
diagnostics of plasma. It is admitted, however, that the model chosen is 


highly idealized. 


AN EXTENSION OF THE PROBLEM 

As emphasized, the previous analysis was restricted to a thin plasma sheet 
with unbounded media on both sides. It would seem to be worth while to 
investigate the case where the plasma sheet is located above a conducting 
plane. Intuitively one would feel that the surface wave of the kind discussed 
above should not be influenced by the presence of the conducting plane if 
the separation was reasonably large, say of the order of a wavelength. On 
the other hand, the space contained by the plasma sheet and the conducting 
plane could be expected to support waveguide type modes. The following 
analysis supports these contentions. 

With reference to the usual Cartesian co-ordinate system, the thin plasma 
sheet is again taken to be located in the plane z = 0. Now, however, there 
is assumed to be a conducting plane boundary at z = —d. As far as the 
harmonically varying fields are concerned, this boundary is assumed to be 
perfectly conducting. The steady or d-c. magnetic field which is applied is 
considered not to be affected by the presence of the conducting medium 
occupying the space z < —d. 

The boundary conditions at the sheet are unchanged, but now an additional 
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boundary condition at z = —d must be met. This latter condition is that 
the tangential components of the electric field should be zero on the plane. 
This suggests writing the Hertz potentials in the following form 


(18) iowa <s for z > 0, 
(19) Il; a ees Om for tas d < Zz < 0, 
(20) = ll = be *%e ** for z > 0, 
(23) gS = B(e® Fe OH) 8 for -—d <2 <0. 
Since dII,/dz and II,* vanish at zg = —d, the boundary condition E, = E, = 0 
at zg = —d is met. Application of equations (9a) to (9d) now leads readily 


to four linear equations involving the coefficients a, a’, b, and b’. Setting the 
determinant of these equal to zero gives 


(22) [(14+V)+MC][(1+V)C+M]+N°C = 0 
where 

B 1 —2ikCd 
(23) V= sea 


It is seen that when 


—ImkCd > 1 


V may be replaced by unity as the equation reduces to equation (12) of the 
previous section. A surface wave does indeed satisfy this condition if d is of 
the order of a wavelength or gréater. Thus, a conducting plane at a reasonable 
distance from the plasma sheet has little influence on the surface wave on 
the structure. 

A first-order correction to the surface wave solution to account for the 
presence of the conducting plane can be obtained by perturbation. For 
example, in equation (22), V is replaced by Vo, where 


: iia 
(24) Vo = —=27kC0d 
l—e . 


and C» is given by equation (13). Solving the resulting quadratic in C leads 


to the first-order solution Cy given by 


—[MP4N?+ (14+ V0)" {LM +N*+ (1+ Vo)P?—4 (1+ Vo)? My) 


fa Po. - 
(25) Co 2(1+Vo)M 


When —Im kCod > 1, Vo may be replaced by unity and the previous solution 


is obtained. 
In addition to a modification of the surface wave solution, equation (22) 


has other solutions. These are the waveguide modes. To illustrate their 
behavior, N is set equal to zero. The equation breaks into two independent 
parts. These may be written in the form 


(26) 1—e ** = —2/(MC), and l+¢ = —2C/M. 
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The first of these equations corresponds to the 7M (transverse magnetic) 
modes while the second corresponds to the TE (transverse electric) modes. In 
general, it appears these transcendental equations must be solved by numerical 
or graphical means. However, if |M| is sufficiently large a simple perturbation 
method may be used to solve for the modes. 
The following substitutions are now made 

2kCd = 2ne+A (n = 0, 1,2,3,...) 
for the 7M modes, and 

2kCd = (2m—1)4r+A’ (m = 1,2,3,...) 


for the TE modes. Then under the condition that A and A’ <1, it is readily 
found that 


€,tkd | 


(27a) ASTM én, e=2 forss me 1253: ..;.. 
and 

(27) Al Moe for m = 1,2,3,... 
Thus 

(28a) Ss re a a for the 7M modes (n = 0, 1, 2,...) 
and 


é li ine £7 
(28) Sn 241 | r 1+5rea) § 


for the TE modes (m = 1, 2,3,...). 
These relatively simple expressions are valid if 
| Mkd| > 1. 


The corresponding results for N #0 are somewhat more involved. Then 
coupling between 7‘E and 7M modes takes place. 

It can be concluded that a thin plasma sheet parallel to a conducting 
plane will support a strongly trapped surface wave whose propagation charac- 
teristics depend mainly on the electrical properties of the sheet. However, in 
addition there are waveguide modes which have to be considered. 
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A MAGNETIC SUSCEPTIBILITY BALANCE AND THE 

TEMPERATURE DEPENDENCE OF THE MAGNETIC 

SUSCEPTIBILITY OF COPPER, SILVER, AND GOLD, 
295°-975° K! 


M. GARBER,”? W. G. HENRY, AND H. G. HOEVE? 


ABSTRACT 
A magnetic susceptibility balance is described which has a reproducibility of 
better than +0.15%. “he temperature dependence of the magnetic susceptibility 


of copper, silver, and gold from 295° to 975° K is reported. 


1. INTRODUCTION 


In an earlier paper, referred to hereinafter as I (Henry and Rogers 1956), an 
apparatus was described which is suitable for making measurements at room 
temperature of the magnetic susceptibility of weakly dia- and para-magnetic, 
machinable alloys with a precision of approximately +0.2%. It was pointed 
out that serious errors may arise in the measurement of the magnetic suscepti- 
bility of metal specimens containing ferromagnetic impurities if any part of 
the specimen is in a fieid in which the ferromagnetic impurities are not reason- 
ably saturated. The minimum allowable field for silver specimens was found 
to be about 6500 gauss. Phases other than solid solutions are in general not 
readily machinable. For this reason a magnetic susceptibility balance capable 
of making precise measurements on cast specimens, cut roughly to length, 
was felt desirable. Some solid solutions such as those formed by the metals 
of the first transition series in copper are unstable at room temperature. As 
it was planned to measure the temperature dependence of some of these 
alloys to determine the effective Bohr magneton value of the solute, an 
apparatus which would maxe measurements above room temperature was 
required. In this paper an apparatus is described which meets all these require- 
ments and which has a reproducibility of better than +0.15%. With a suitable 
Dewar the apparatus could be used for low-temperature measurements. 

The basic components of the apparatus are: a modified Varian Associates 
electromagnet with pole pieces 12 in. in diameter; a Varian Associates power 
supply; a combined mechanical and electronic balance which is a modified 
version of that suggested by Mauer (1954). In order to avoid the necessity 
of accurately machined specimens the pole pieces were designed using the 
electrolytic tank technique to produce a volume of approximately constant 
B,0B,/dx where x is the vertical direction. The furnace consists of two 
globars inside a silver shield. 
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There have been numerous magnetic susceptibility devices described in the 
literature. In general the investigators have been interested in either results 
of much lower precision or only in reasonably high precision relative to, say, 
aroom temperature measurement. Stevens, Cleland, Crawford, and Schweinler 
(1955), who measured the susceptibility of germanium at low temperatures, 
have recently described an apparatus capable of a degree of precision similar 
to that of the present apparatus. The problem of irregularly shaped specimens 
was avoided by making the measurements relative to oxygen. This method, 
which is satisfactory at low temperatures, would not be satisfactory at tempera- 
tures at which oxidation would take place. Also the method loses precision as 
the temperature increases because of the decrease in the susceptibility of 
oxygen. The specimens were placed in the fringing field of an electromagnet 
and no attempt was made to shape the pole pieces. 

In the present paper the results of the measurement of the variation of the 
magnetic susceptibility of copper, silver, and gold from room temperature to 
975° K are presented. There have been only three previous investigations in 
which all three metals were measured (Honda 1910; Owen 1912; Shimizu 
1937). The materials used by Honda and Owen were sufficiently impure to 
cast doubt on their results and, further, no numerical values were given. Ross 
(1934) and Bitter, Kaufmann, Starr, and Pan (1941) found that copper 
becomes less diamagnetic as the temperature increases, whereas Shimizu found 
the reverse. The present work confirms the results of Ross and Bitter et al. 
regarding the sign and magnitude of the temperature dependence. The absolute 
value of the diamagnetic susceptibility of silver is found to go through a 
minimum at approximately 750° K. Up to 750° K the silver becomes less 
diamagnetic; Shimizu found that silver becomes more diamagnetic in that 
temperature range. Rao and Savithri (1942) found the temperature dependence 
of the susceptibility of silver to be less than their limit of error, which was 
approximately +1.0%. The sign of the temperature dependence of gold is 
found to differ from those found by Shih (1931) and Shimizu. The present 
work appears to agree with that of Vogt and Gerstensberg (1959) but they 
give no data and hence a good comparison cannot be made. The differences 
are probably due to the differences in purity of the materials, and in the 
method of handling ferromagnetic impurities. 


2. POLE PIECE DESIGN, CALIBRATION AND CORRECTION FOR 
FERROMAGNETIC IMPURITIES 

(a) Pole Piece Design 

The pole pieces were designed by the method of electrical analogy (Malavard 
1947). Eight electrostatic models were investigated, the best of which is 
described here. References to other pole piece designs are given by Bates 
(1951). 

Assuming the pole faces to be magnetostatic equipotentials, profiles are 
sought in the x—z plane such that along the vertical symmetry axis (x-axis) 


(2.1) B, 0B,/dx = constant. 
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Since the pole pieces are long in the y direction (horizontal and parallel to 
the faces), By = 0 and 0/dy = 0. Equation (2.1) represents the condition 
that all elements of the sample lying in the x—y plane experience the same 
vertical force per cc, provided the magnetization is a linear function of the 
field. A discussion of the more general case in which ferromagnetic impurities 
are present is given in I. For elements of the sample which are not in the 
x—y plane, the complete expression must be used in place of (2.1) (1), 


(2:2) —B, 0B,/dz+B, 0B,/dx = constant. 


As a first approximation the variations with respect to z are neglected and 
attention is confined to (2.1). The relationship between B,(x,0) and x is thus 
given by 


(2:3) B, = (ax+6). 


In order to obtain a relationship between the gap width and x, it is assumed 
that B,(x,0) varies inversely with the gap width and is given by 


(2.4) B, = d/2z, 


where d is a constant and 22 is the gap width. On combining (2.3) and (2.4), 
it is found that 


(2.5) 25 = d/(ax+6)}. 


Equation (2.5) represents the trial equipotential surface which was used as 
a starting point for the electrostatic tank experiments. The constants a and 6 
were determined by applying the following conditions: at the base of the 
specimen (x = O0cm) the gap width is 5.0cm; at the top of the specimen 
(x = 2cm) the field is half that at x = Ocm. The latter condition ensures 
that for the lowest magnet current used the field will be strong enough (6500 
gauss) everywhere over the sample to saturate ferromagnetic impurities and 
vet allow for a sufficiently strong force per unit volume. The relation between 
z and x is 
(2.6) 2 = 50/(8—3x). 


In the electrostatic model, the gap was allowed to narrow to 3.1 cm, being 
limited by furnace design, and below this point the pole pieces were tapered 
apart at an angle of 60° in order to provide a stronger field. The model showed 
that the above design led to a nearly constant force density over the desired 
working length of about 2 cm; that the tapering did not adversely affect the 
desired field over this length, and that replacing the curved shape by straight 
line segments, in order to facilitate machining, would have no significant 
influence. 

A cross section and plan view of the pole piece are shown in Fig. 1. The 
pole pieces were made from castings kindly supplied by the Bureau of Mines, 
Ottawa, Canada. The analysis of the material is: C, 0.016; Mn, 0.20; Si, 0.03; 
S, 0.023; P, 0.015; Cr, 0.03; Ni, 0.06; Mo, <0.01; V, 0.003; Al, 0.035; Cu, 
0.15; Sn, 0.01; Co, 0.003 wt.%. 
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Fic. 1. A cross section and plan view of a pole piece. 
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Fic. 2. The low and high fields and the specimen position. 


As in I, two fields were used; the lowest was determined by the condition 
that no part of the specimen should be in a field of less than 7000 gauss. In 
Fig. 2, the variation of the square of both the low and high fields with x is 
shown. The mean values of the slope of the squares of the fields over the 2-cm 
length are —84.0X10° and —25.4X10° gauss? per cm. The dotted curve was 
that found by the method of electrical analogy. The constant of proportionality 
between the electric field and the field of magnetic induction was found by 
fitting the two curves in the high field region. 

The relationship between the field strength and the gap width is better 
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approximated by B, = d/(2z)*** than by equation (2.4). The pole pieces were 
not redesigned using this information because the degree of the success of the 
method was not known beforehand. This does suggest, however, a self-con- 
sistent method for designing pole pieces. Thus as a first step an approximation 
such as (2.4) would be used, then the actual relationship between B, and z 
would be determined experimentally in the tank. The design expression would 
then be reformulated using the experimentally determined counterpart to 
equation (2.4). 


(b) Field Calibration 

The calibration was made using copper specimens free from interfering 
ferromagnetic contamination, of constant diameter, and varying length (1.6 
to 2.0 cm). The absolute mass susceptibility of the copper had been determined 


earlier (1) to be —0.0859,X10-* c.g.s. e.m.u. The force F, in the vertical 
direction is given by 


(2.7) F, = (m/a’)(x—«/p) 


where m is the mass of the specimen, 1/a’ is the mean value of (2.2) over the 
specimen. The mass susceptibility of the specimen is x; « is the volume sus- 
ceptibility of the atmosphere and p is the density of the specimen. All the 
experiments were done in helium and hence «/p could be neglected in com- 
parison with the mass susceptibility of copper. 

To determine the constants of the apparatus, the bottoms of the copper 
specimens were located in a fixed position relative to a fixed mark on the 
pole pieces, approximately 1.3 cm above the minimum gap. The constants for 
the low and the high field, shown in Fig. 3, for the copper specimens of various 
lengths were calculated from the relationship 
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Fic. 3. The low and high field calibration. 
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_ mass(gram) X85.99 - 
7 force (mv) 


(2.8) 


The constant a@ is related to a’ of equation (2.7); a = —a’g (mg/mv) X10%, 
where (mg/mv) is the constant, milligrams per millivolts, of the balance 
discussed in Section 3 and g is the gravitational constant. The a@ for the high 
field is constant to within +0.1% for specimens with lengths between 1.6 


and 2.0 cm. 


(c) Ferromagnetic Impurities 

When (2.2) is a constant, the extrapolation for ferromagnetic impurities 
should be made against 1/(Bmax+Bmin), Where Bmax and Byin are the maximum 
and minimum fields acting on the base and the top of the specimen at the 
particular setting of the magnet current. This result follows from an argument 
very similar to that given for the short specimen Gouv method, in I. 


3. BALANCE POWER SUPPLIES, AND OPERATION 

The combination mechanical and electronic balance with electronic velocity 
damping is very similar to that described by Mauer (1954). The major modi- 
fication consists of the addition of a balancing bridge circuit to the output 
which converts the instrument to a null-type balance. By making other minor 
modifications the sensitivity of the device has been increased by a factor of 
about 20. 

The modified circuit is shown in Fig. 4. In order to improve the stability 
of the circuit, all the resistors are the manganin wire-wound type, except those 
in the bridge galvanometer sensitivity circuit. The high voltage is supplied by 
batteries; power supplies using voltage regulating tubes were found unsatis- 
factory. The filament current is supplied by a separate, transistorized power 
supply with a stabilization factor of about 2000. 

The light source is operated by a well-shielded, battery-stabilized, d-c. 
supply. The light source and split phototube are completely shielded from the 
room lights. Separate shielded cables from the Mauer circuit to the mechanical 
balance are used, one for the phototube output and supply and one for the 
solenoid and light source. This was necessary to remove pickup in the photo- 
tube bridge circuit of both the 60-cycle ripple in the light source supply and 
the irregular variations of the current in the solenoid. 

A general view of the apparatus is shown in Fig. 5. The balance is a Sar- 
torius SM10, semimicro balance with a sensitivity of 0.01 mg and a capacity 
of 100 g on each pan. The balance is operated at constant load by adjusting 
the weights on the pan to compensate for the weight of the specimen. The 
solenoid, 9.52.2 cm, which consists of two layers of 0.01-cm enameled wire, 
is enclosed in a Pyrex case with tungsten leads to facilitate the evacuation 
of the system. The electrical connections to the solenoid are made through 
loosely coiled vertical springs. The balancing bridge circuit eliminates the 
influence of the springs. The cylindrical Alnico magnet 12.5 X0.9 cm is fastened 
to the base of the chamber. The radioactive sources were found essential to 
eliminate electrostatic forces. To insert a specimen the pin is withdrawn and 
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Fic. 5. A general view of the apparatus. 
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the inner quartz tube removed through the open end at the bottom of the 
outer quartz tube. Access to the weight pan is through a port. 

The solenoid, Alnico magnet system has a sensitivity of about 50 mg per 
ma. Without the bridge circuit a deflection of the end of the balance arm of 
about 0.25 mm is required to produce a current of 2ma. The use of the 
balancing bridge circuit reduces the deflection to the limit of the reproduci- 
bility of the cathetometer, +0.01 mm. This feature is important in applica- 
tions where the rate of change of B, 0B,/dx is significant (see Section 5). In 
Fig. 6 the millivolt, force relationships are shown with and without the 
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Fic. 6. A comparison of the Mauer (1954) circuit, with and without the addition of the 
balancing bridge circuit. 


bridge circuit. It is seen that with the bridge circuit the milligrams per millivolt 
are more constant and the useful range is extended. The variation of Aw/Amv 
without the bridge circuit is presumably due to the change of position of the 
magnet in the coil. Some of the scatter in both cases may be due to the dis- 
turbing of the system when weights were either added or subtracted. The 
results with the bridge circuit were reproduced within experimental error 
after an interval of one year. 

The method of operation is as follows. The balancing bridge circuit is iso- 
lated from the Mauer circuit and the former is balanced. This is not altered 
again. With the balancing birdge circuit connected, the input grid of the 5691 
is grounded and the output of the 12SN7 is balanced. With the input grid not 
grounded and with the beam clamped near the reference position, the photo- 
tube bridge is balanced. The maximum deflection on the Cambridge galvano- 
meter (1765 ohms, 20 mv full scale) in the balancing bridge circuit, as the 
beam passes through the point of maximum lift, is called the reference deflec- 
tion. The beam is then released and the resistance in the battery leg of the 
bridge circuit is adjusted to bring the galvanometer back to the reference 
deflection. The voltage drop across a 500-ohm resistor in the balancing bridge 
circuit is then read with a Cambridge Vernier-type potentiometer. 
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The magnet is a model V-4012A water-cooled (closed system) electromagnet 
made by Varian Associates with a 30.5 cm diameter pole piece. The magnet 
was supplied with a widened yoke and correspondingly increased coil box 
separation to accommodate the pole pieces. The number of ampere turns at 
maximum current is 80,000. The field is set by setting the current. The magnet 
is demagnetized by reversing the current to a value that was found experi- 
mentally to result in an effectively zero field when the current is returned to 
zero. The magnet power supply is a model V-2100, self-regulating power 
supply made by Varian Associates which maintains the current constant to 
+0.001% at 0.8 amp but requires a small amount of manual adjustment at 
2.0 amp. The voltage drop across a Leeds and Northrup 0.1l-ohm, 15-amp 
capacity, standard resistance in series with the magnet current is read con- 
tinuously during a measurement. 

The influence of the magnetic field on the operation of the balance is mini- 
mized by the following features. The balance is placed about 1 meter above 
the center line of the pole pieces and separated from the magnet by two 
steel plates 1X1 m and 6 mm thick. The balance is enclosed in a steel cylinder 
which forms part of the vacuum system; the solenoid, contrary to Mauer’s 
design, is fastened to the balance arm. The two sources of spurious forces which 
might lead to errors are the possible solenoid, magnet interaction, and a 
weight, magnet interaction. Regarding the solenoid, magnet interaction, it is 
found that the current through the solenoid required to balance the force on 
a standard copper specimen, at maximum field strength, is independent of 
the initial current through the solenoid. The initial current may be varied in 
either direction by suitably loading the pan. Regarding the weight, magnet 
interaction, it is found that when a 100-g weight of the same material as the 
other weights is added to both sides of the balance, one on the pan and the 
other inside the steel cylinder surrounding the balance, there is no change in 
the force on the specimen holder at maximum field. 

To maintain precision, it is necessary to take some care in making the 
measurements. The magnet is always cycled in the same way. The zero field 
reading is interpolated to take account of drift. The sample holder is measured 
regularly and boiled if necessary in aqua regia to free it from ferromagnetic 
dust. 


4, FURNACE AND TEMPERATURE CALIBRATION 

A section and a plan view of the furnace are shown in Fig. 7. The heaters 
for the furnace are two globars, manufactured by the Carborundum Company, 
with an active length of 10cm and a diameter of 0.8 cm. These are operated 
on a-c. current. The globars are rubber-mounted and mechanically isolated 
from the furnace assembly; the vibration which does occur does not influence 
adversely the operation of the apparatus. The globars are operated in series, 
and the power consumed at 500° C is about 200 watts. 

The temperature control system is very similar to that described by Berry, 
Henry, and Raynor (1951) and it controls the temperature in the absence of 
the field to +0.5° C. Two chromel-constantan thermocouples in series drive 
the temperature control circuit. As the magnet field is changed an e.m.f. is 
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Fic. 7. The globar furnace. 


induced in the control thermocouples and some manual anticipation is required 
to maintain the temperature within the limits of +1.0° C. The e.m.f. of the 
calibrated, platinum-rhodium, specimen thermocouple is read continuously 
during a measurement on a Rubicon portable potentiometer. The maximum 
temperature difference between any two parts of the specimen, as measured 
by a probe thermocouple inserted in a hole along the axis of the specimen 
occurs at 400° C and is 4° C. 

The temperature calibration of the specimen thermocouple was checked by 
observing the change in the susceptibility of tin, lead, zinc, and antimony 
on melting. These metals all showed an abrupt change at their melting points. 
The melting points were taken to be 232, 327.5, 419.5, and 630.5° C respec- 
tively. It is estimated that the absolute error in the mean temperature of the 


. 


specimen is not greater than +1.5° C. 


5. ERRORS 
In this section the errors resulting from the specimen not occupying the 
same volume as the calibrating specimen are considered. Other sources of 
error are also considered. The expression for the lateral force is given. 
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It is convenient to let B,(x,0) = Bwhere B = B(x). The errors are expressed 
in terms of B® only and its derivatives. 


(a) Positioning Error in the x Direction 

The force F, in the x direction on a specimen in the calibrating position is 
approximately equal to mx,/2/AB*, where x, is the true mass susceptibility, 
m is the mass; / is the length, and AB? is the difference in the squares of the 
fields at the top and the base of the calibrating specimen. If the specimen 
is now displaced by an amount 6x cm, the apparent susceptibility xq is given 
by 


(4.1) Xa = 21F,/mAB? 


where F, is now the new force. It follows that the true susceptibility is given 
by 


AdB’ /d; 
(4.2) X: = i x). 


For the low and the high field (AdB?/dx)/AB? is approximately 0.03 per cm. 
At 975° K there is an error approximately equal to or less than the experi- 
mental error which arises because of the increase in length of the specimen 
holder. No correction, however, has been made to take account of this effect. 

As indicated in Section 3, the bridge circuit reduces the deflection of the 
balance under load to +0.01 mm. The Mauer circuit, if used alone, could 
result in errors of the order of 0.1%. 


(b) Positioning Error in the z Direction 
Following the method of paper I, the expressions for B,(x,z) and B,(x,z) 
are given by 


3 
rise = (Bua) 
(4.3) Xx * x 
1/d°B\ » 
B,(x,z) roe Se 


where B = B,(x,0) as before. In equation (4.3) the coefficients of the various 
powers of z are functions of x which are determined by the requirement that 
curl B = 0, and div B = 0. 

The force in the x direction, on an element of volume is given by, keeping 
terms up to 2 only, 


a 2 eS ee 
, dxdydz 2 dx ‘2 2 dx” ' BY \ dx /\dx'/ 2B \dx/ FT 


where «x is the volume susceptibility. Let the part in the square bracket be 
denoted by f(x). Then if the center of the specimen is displaced in the z direction 
from the x-axis by an amount g cm, the force on a circular slice of thickness dx 
and diameter d is given by 


iF, . a We ee 
(4.5) (=) = gis ee (“sata ‘af (x). 
x] ' 
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It follows that the true susceptibility of a slice, provided the correction is 
small, is given by 


(4.6) Xe = Xall —G?f(x)/(dB?/dx)). 


The mean value of the second term in the square bracket over the length 
of the specimen is approximately 0.1 q? for both the low and the high fields. No 
correction was made to take account of a z positioning error. Displacements 
of the order of 0.1 cm may have occurred. 


(c) Diameter Error 

If the specimen diameter is larger by an amount Ad cm than that of the 
calibrating specimens, then it may be shown by (4.5) with g = 0, that for a 
slice of thickness dx 


(4.7) Xt = Xall —d(Ad)f(x) /(8dB?/dx)]. 


No correction to the measured results was necessary since the maximum 
value of Ad was 0.05 cm. 


(d) Force in the z Direction 

If it is desired to know the maximum susceptibility which may be measured 
before trouble through, for example, contact with the walls, is liable to arise 
from the horizontal displacement of the specimen by the field, the force in 
the z direction must be known as a function of g. It may be shown that the 
force in the z direction on a slice of thickness dx is given by 


| dF.) _= ie * Cag 
48) (4) gi *| arte a J: 


This effect is considerably more noticeable in the high field. For the present 
pole pieces, the force in the z direction on a specimen 1 cm in diameter and 
2cm long is approximately —10X10° xq mg/cm. 

The error in the measurement of the mass is negligible. Reference to the a 
curves, Fig. 3, shows that the length of the specimen need only be determined 
to within 1 mm to obtain a to better than 0.05%. 

The Aw/Amv curve shown in Fig. 6 has a scatter of about +0.05% but as 
noted in Section 3 some of this scatter was probably due to opening and 
closing the system to either add or remove weights for each point. The repro- 
ducibility of the force during an actual experiment is much better than 0.1%. 

As an indication of the long term stability and the reproducibility of the 
apparatus the results obtained on a single specimen, which was made as a 
standard, may be quoted. A silver—cadmium alloy containing approximately 
30 at.% cadmium was measured three times with an interval of about 6 
months between the measurements. The maximum deviation from the mean 
susceptibility was 0.1%. 

The limiting feature of the apparatus appears at present to be the provision 
for positioning the specimen in the sz direction. Further improvements are 
probably not worth while, however, until the purity of the specimens is 


improved. 
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6. TEMPERATURE DEPENDENCE OF THE MAGNETIC SUSCEPTIBILITIES 
(a) Copper 

The susceptibility results in this and following sections are all c.g.s. e.m.u. 
per gram X 10°. 

Seven specimens were measured. The specimens were made from four 
separate lots of copper purchased from the American Smelting and Refining 
Company (A.S. and R.) and analysis of which, according to the manufacturer 
is: Fe <0.00007, Sb <0.0001, Pb <0.0001, Sn <0.0001, Ni <0.0001, Bi 
<0.00001, Ag <0.00003, As <0.0002, Cr <0.00005, Si <0.00001, Te <0.0002, 
Se <0.0001,S <0.0001, O2 <0.0001 wt.%. The iron analysis has been checked 
in our laboratory and found to vary from 0.0001 to 0.00002 wt.%. An oxygen 
and hydrogen analysis was made by the Department of Mines and Technical 
Surveys and found to be: oxygen, 0.0001 +0.00003; hydrogen, 0.00004 wt.%. 
The specimens were cut to length on a lathe from the “as received”’ rod; they 
were then etched heavily with 10 N nitric acid to remove ferromagnetic 
contamination. The specimens were annealed at 575° C for 5 days and air- 
cooled. It was found in I that this annealing treatment did not alter the 
susceptibility; it was done in the present work to help put into solution any 
small amounts of ferromagnetic material not removed by etching. 

The a@ curves (Fig. 3) were determined from a number of specimens of 
varying length which were not necessarily used to determine the high tem- 
perature results. [t was necessary to assume that there were negligible amounts 
of ferromagnetic impurities in these specimens; some selection of specimens 
was made. 

In the high temperature experiments the a curves were used even for the 
room temperature susceptibility determinations. The average of the room 
temperature susceptibilities of the seven specimens is —0.08600, which is 
essentially the same as the calibrating value of —0.0859 . These room tem- 
perature results lend strong support to the assumption regarding negligible 
amounts of ferromagnetic impurities in the calibrating specimens. There was 
no indication of ferromagnetic impurity in the specimens used; the low and 
the high field determinations agreed within experimental error. There was 
no indication of contamination as a result of being at an elevated temperature; 
the subsequent room temperature values agreed within experimental error. 

The temperature dependence of the susceptibility of the seven copper 
specimens is shown in Fig. 8 along with the data of Ross (1934), and Bitter, 
Kaufmann, Starr, and Pan (1941). Ross used electrolytic copper of unknown 
purity; the specimens were prepared in a hydrogen atmosphere. A para- 
magnetic impurity such as iron in solution causes the slope to be smaller 
and hence this does not explain the difference in slope. Bitter et al. extrapolated 
the results of copper-iron alloys of higher concentrations to zero concentration 
to obtain their data for pure copper. Their specimens were melted in a hydrogen 


atmosphere and chill-cast in vacuum. They also checked their results against 
a specimen of A.S. and R. copper. They were not attempting to make an 
absolute measurement on copper but were interested in relative results 
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Fic. 8. The temperature dependence of the magnetic susceptibility of copper, 295°-975° K. 


(Bitter and Kaufmann 1939) and hence the difference in absolute value 
between their work and the present work is perhaps not significant. 

In order to estimate the error in an individual measurement, a smooth 
curve was drawn through the data on a large graph such that the sum of 
the deviations was about 0.0001 9. The mean error was then calculated as 
(v?/n)}, where v is the deviation and m is the number of points. The mean 
error is found to be +0.00014. Some of the error is due to the inclusion of 
some early results in which the temperature was not regulated as well as 
later on. The values in Table | were taken from the smooth curve. 


TABLE I 


Magnetic susceptibilities per gram of copper, silver, 
and gold X 108 


Gold 





Temp., °K Copper Silver 
295 0.0859, 0.18072 0.1430; 

400 0.08463 0.1792, 
500 0.08333 0.1780, 0.1428, 

600 0.08183 0.1769. 
700 0. 0802, 0.1763, 0. 1424, 

800 0.07853 0.1760; 
900 0.07675 0.1776, 0.14195 
0.0753; 0.17865 0.1417; 


975 


In Fig. 9 the low temperature data of Ross (1934), Bitter, Kaufmann, 
Starr, and Pan (1941), Bowers (1956), and Pugh, Coles, Arrott, and Goldman 
(1957) are shown with the present high temperature data. The last three 
investigators have used the same copper (A.S. and R.). The data is presented 
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Fic. 9. The temperature dependence of the magnetic susceptibility of copper above and 
below room temperature. 


as a check on the room temperature value of the slope of the susceptibility 
curve for copper and to show the degree of agreement between the investi- 
gators. The present work agrees moderately well with that of Ross, Bitter et 
al., Pugh et al., but not with Bowers. 


(6) Silver 
Three specimens were measured. The specimens were made from twice- | 
electrolyzed silver (Harrow and Henry), the analysis of which is: Fe 0.0002, i 
| 
: 


: 
i 





Si 0.0002, Pb <9.0005, Cu 0.0001, Mg <0.0001, Cd <0.0002, Mn 
<0.0001 wt.%. The specimens were melted under vacuum in quartz crucibles 


measurements agreed within experimental error. 

The temperature dependence of the susceptibility of the three specimens 
is shown in Fig. 10. The same procedure was followed as for copper to deter- 
mine the error in an individual measurement. The mean error is +0.0001:. 
The values in Table I were taken from the smooth curve. 

The minimum in the absolute value of the susceptibility at about 750° K 
was unexpected and further work was done to rule out the possibility of oxygen 
contamination being the cause. Steacie and Johnson (1926) found that the 
solubility of oxygen in silver is a maximum at about this temperature. The 
vacuum-cast specimens were found earlier (I) to have the same susceptibility i 
as degassed specimens. It was possible, however, that some oxygen con- 


and chill-cast (Craw and Henry 1956). After casting the Si content was ; 
found to be 0.0004 wt.%. The specimens were cut to length and then given 
a 2-minute treatment with concentrated nitric acid. They were then sealed : 
in Pyrex, in vacuum, and annealed for 7 days at 575° C, and then air-cooled. : 
It was found in I that the room temperature susceptibility was not influenced 
by this procedure. 

There was no evidence of ferromagnetic impurity; the low and high field 
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Fic. 10. The temperature dependence of the magnetic susceptibility of silver, 295°-975° K. 
© normal specimens; ™ 12 hours in air at this temperature; 1) degassed specimens. 


tamination of the specimen occurred during the measurement. The specimens 
marked in Fig. 10 were held overnight im situ, at the temperatures indicated 
in air and then measured. The specimen was then degassed overnight, in a 
vacuum of approximately 5X10-> mm, at 875° K and measured. It was held 
a further 24 hours at 875° K in vacuum and measured again. The two squares 
which almost superimpose represent these latter measurements. The lower 
square at about 750° K represents the result of the measurement at that 
temperature after the degassing at 875° K. The magnitude of the change as 
a result of the contamination indicates that the minimum could not be the 
result of oxygen contamination. The additional fact that indicates that there 
was little or no contamination during the experiments on the other specimens 
is that the room temperature values before and after being at temperature 
agreed within experimental error. 

Pugh and Ryan (1958) found the susceptibility of silver to vary from 
—0.183 to —0.181 from 20.4 to 295° K. The slope has the same sign as in 
the present work. 


(c) Gold 

Three specimens were measured. The specimens were made from twice- 
electrolyzed gold (Harrow and Henry), the analysis of which is: Mg <0.0001, 
Fe <0.0001, Si <0.0001, Cu <0.0001, Pb (estimate) 0.0001, Ca 0.0001, 
Ni <0.0001, Ag 0.0001 wt.%. The gold was melted under argon in quartz 
crucibles and chill-cast. The specimens were cut to length and then given a 
heavy etch with hot (60° C), concentrated aqua regia. They were then sealed 
in Pyrex, in vacuum, and annealed for 7 days at 575° C. 

There was no evidence of ferromagnetic impurities; the low and high field 
determinations agreed within experimental error. 

The results of the measurements on the three specimens are shown in 
Fig. 11. The error in an individual measurement was found to be +0.0000;. 
The values in Table I were taken from the smooth curve. 
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Fic. 11. The temperature dependence of the magnetic susceptibility of gold, 295°-975° K. 


7. DISCUSSION OF RESULTS 

In the present work the mean susceptibilities of the silver and gold specimens 
at 295° K are found to be —0.18072+0.0000, and —0.1430;+0.0000; re- 
spectively where the errors are the root mean square errors relative to copper 
determined from the measured values at 295° K. The fact that the silver 
value has decreased and the gold value has increased compared to I suggests 
that the differences are not due to some systematic error in making the measure- 
ment. In I, one specimen of twice-electrolyzed silver was included among 
the pure specimens. This specimen had the lowest absolute value of the 
group (—0.18109) but was not considered a ‘‘wild’”’ one. A measurement on 
another specimen of twice-electrolyzed silver made after the publication of I 
on the apparatus used in I yielded a value of —0.1807;. These latter measure- 
ments taken in conjunction with the results of the present work suggest that 
the increased purity of the specimens is somehow responsible for the difference 
between the value for silver in I and the present value. The difference in the 
values for gold is also probably due again to increased purity. The specimen 
of highest purity in I, made from proof plate, gave a value of —0.1429p in 
good agreement with the present value. Although the susceptibilities of the 
other three specimens were corrected in I for an iron content of 1.2 p.p.m. 
by weight, the presence of additional transition metal impurities at this level 


could account for the difference. 
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THE IDENTIFICATION AND HALF LIVES OF 
FISSION-PRODUCT Rb® AND Rb*! 
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K. FRITZE AND T. J. KENNETT 


; 
ABSTRACT : 

The existence of two new rubidium isotopes, Rb and Rb®, has been estab- t 

lished and their half lives measured. The half lives of these short-lived fission t 
products were determined using a technique of timed precipitations. The values t 
obtained for Rb and Rb* were 5.3+0.5 sec and 5.6+0.5 sec respectively. The t 

half lives of the strontium and yttrium daughters were also determined. The f 
strontium isotopes were studied by observing the decay rate of a characteristic i 

y-ray peak. For Sr®, the decay of the 1.37-Mev line gave a value of 2.71+0.02 ‘ 

hr. A y-ray peak at 590 kev, which was found to be associated with Sr, decayed i 

with a half life of 7.54+0.06 min. The half lives of the yttrium daughters were : 
determined by 8 counting. The values found for Y® and Y*% were 3.53+0.02 hr i 

and 10.1+0.1 hr respectively. ' 

i 

(A) Rb® AND Rb*® 

The existence of Rb® and Rb® was implied by the discovery of their ' 
respective parents Kr® and Kr. The latter were established to be fission ; 
: ° ° ° ° . ° t 

products through gas sweeping experiments carried out with irradiated uranium i 
solutions by Dillard et a/. (1951). Hahn and Strassmann (1940) reported an i 
80-sec activity associated with rubidium fission products which has been 
assigned to Rb® in recent nuclear data compilations. | 
To investigate the half lives of Rb® and Rb”, a method similar to that : 
described by Sugarman (1949) was used. This technique consists of dividing 
i 


a sample of fission products into portions from which certain elements are 
precipitated. The precipitations are carried out serially with a time separation 
between each. At a later time, a daughter product is extracted from the 
precipitate, purified, and counted. Since the quantity of the daughter found in 
each sample is proportional to the amount of the parent present at the time 
of precipitation, one has a method of determining the half life of the parent. 
This method has the distinct advantage of unique mass assignment, a virtue 
not shared by direct measurements of the half life. 

The usefulness of the timed-precipitation technique depends strongly upon 
the nature of the precipitation reaction. This reaction must be fast in com- 
parison to the intervals between successive precipitations in order to minimize 
timing uncertainties. Also, and of greater importance, the precipitate must 
not carry any daughter activities. For this work it was found that the pre- 
cipitation of rubidium perchlorate fulfills both requirements, provided that a 
sufficient quantity of rubidium carrier is used. The co-precipitation of the 
carrier-free daughter products which was determined by using Sr®° and Y®, 
was found to be of the order of 10-*%. 

To determine the half lives of Rb® and Rb, the rubidium perchlorate 
precipitate can be analyzed for the strontium (2.7-hr Sr®, 7.5-min Sr) or 
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yttrium (3.5-hr Y®, 10.1-hr Y%) daughters. The yttrium daughters are the 
obvious choice because of their convenient half lives and the uncomplicated 
decay to stable or long-lived zirconium isotopes. To obtain the maximum 
activity of Y®, the first strontium-yttrium separation should be performed 
4.5 hours after the irradiation. The decay curve of the purified yttrium samples 
will then be composed of the 3.5-hr Y® and 10.1-hr Y® plus a contribution 
from the 91 mass chain (50-min Y"”" and 58-day Y*). By unravelling the 
Y*” and Y* components in the decay curves and correcting for chemical 
vields, the half lives of Rb® and Rb® can be easily calculated from the time 
schedule of the successive precipitations. 

Preliminary measurements conducted for Rb® and Rb® indicated that both 
half lives were less than 10 seconds. Since the krypton precursors (3-sec Kr®, 
2-sec Kr*) should be allowed to dee:y almost completely, the time available 
for the rubidium perchlorate precipitations ranges from 15 to 30 seconds after 
the irradiation. Because of this rather short time interval only two precipi- 
tations were performed for each run. 


Procedure 

Samples containing 2 ml of uranyl nitrate solution (93% enriched in U?*) 
were irradiated in sealed polyethylene bottles using the rabbit system of the 
McMaster reactor. Immediately after delivery the irradiated solution was 
mixed with 400 mg of rubidium carrier contained in 18 ml of 2 M HCl. The 
resulting solution was divided into two parts and the first one added to approxi- 
mately 10 ml of 20% HCl1O, contained in a large millipore-filter apparatus 
of the type normally used in radiochemical work. Immediately after precipi- 
tating the rubidium perchlorate, the liquid was drawn through the filter. 
Exactly the same procedure was then used approximately 10 seconds later 
for the second part of the sample. Both precipitates were washed with 10% 
HCIO,. The two precipitates were then rinsed into 100-ml volumetric flasks, 
each of which contained 20 mg of strontium and 10 mg of yttrium carrier in 
9 ml of 9 M HCI. After the rubidium perchlorate had dissolved, one 10-ml 
aliquot was taken from each sample for chemical-yield determinations. The 
remaining solutions were then quantitatively transferred into beakers and 
after evaporation to near dryness, analyzed for yttrium in the following steps: 
Y(OH); precipitation 4.5 hours after the end of the irradiation, tributyl- 
phosphate-solvent extraction of the dissolved Y(OH); (Stanley 1954), Y(OH); 
precipitation from the back-extracted yttrium solution, passage of the re- 
dissolved yttrium in 10 M@ HCl over Dowex 1, Y(OH),: precipitation of the 
eluant, and finally an yttrium oxalate precipitation. The yttrium oxalate was 
filtered over Whatman No. 42 paper, washed with 50% C2H,OH, dried, and 
mounted on a turntable for automatic counting of two samples. Here the 
samples were alternately counted for 10-minute periods each followed by a 
0.3-minute changing interval. The counting was done using either a §-scintilla- 
tion detector or an end-window proportional counter. Figure 1 shows a typical 
decay curve as obtained in this way. After the counting was finished, the 
yttrium oxalate samples were ashed and weighed as Y:O;. The chemical 
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Fic. 1. A typical decay curve for the yttrium daughters of a fast precipitation of fission- 
product rubidium. The solid lines indicate the components obtained from a least-squares 
analysis of the data. 





yields of the initial rubidium precipitations were determined as rubidium 
perchlorate. Each initial rubidium perchlorate precipitate usually contained 
slightly more than 40% of the total rubidium (= 80-90% yield for the 
individual rubidium perchlorate precipitations) whereas the yttrium yields 


were always around 60%. 


Results 

The decay curves were analyzed with a least-squares program using a 
Bendix G-15 digital computer. The data, which were accumulated for at least 
70 hours for each run, were analyzed into four exponential components. The 
components sought for from the decay curve were 50-min Y®”, 3.5-hr Y®, 
10.1-hr Y“, and a combination of 58-day Y* and a constant background. The 
least-squares program was designed to yield the counting rate and associated 


error for each component at the time when counting began. A simple calcula- 
tion was then used to find the counting rate at the time of the first strontium 
yttrium separation. Figure 1 shows a decay curve with the components found 


using the least-squares analysis. 

Table | gives the irradiation conditions and half lives obtained for Rb” 
and Kb”. The errors quoted with each measurement are only those arising 
from the least-squares analysis of the decay curves. Errors associated with 
uncertainties in the timing of the precipitations can be inferred by the repro- 
ducibility between runs. The largest uncertainty in any half-life determination 


of this kind arises from the fact that the washing of the precipitate removes 
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TABLE I 
Flux Duration of RbCIO, Rb” Rb* 
Us (neutrons irradiation precipitation times half life half life 
Run (ug) cm? sec™!) (sec) (sec) (sec) (sec) 
Ist ppt. 2nd ppt. 
1 100 3X10" 30 23 33 5.8+.3* 5.7+.1* 
2 200 3X10" 20 20 28 5.64.2 5.74.1 
3 200 4X10" 20 19 33 4.8+.1 5.8+.1 
4 200 4X10" 20 17.5 28.5 5.14.1 5.4+.1 
Average 5.3 5.6 








*Quoted errors are only those obtained from the least-squares analysis. 


some of the daughter atoms which have grown in since the time of the initial 
precipitation. In order to estimate the magnitude of this effect, rubidium 
perchlorate was precipitated 1 minute after the end of an irradiation and 
then washed very carefully with 10% HCI1O,. It was left on the membrane 
filter in order to let the Rb® decay. The precipitate was then washed again 
with 10 ml of HCIO, for a contact time of approximately 5 seconds. Both 
the precipitate and filtrate were analyzed for Sr*. It was found that for 
even this short a contact time, 40% of the strontium was washed out of the 
precipitate.* However, the relatively small scattering of the results in Table I 
indicates that the errors due to these effects cancel out to a large extent, 
particularly in view of the fact that in runs | and 2 the washing was done 
after both precipitations were performed whereas in runs 3 and 4, the washing 
was done after each precipitation. Consideration of the small number of 
measurements and the difficulty of assigning errors to some aspects of the 
experiment lead us to suggest the following half lives: 


Rb*® = 5.3+40.5 sec 
Rb*® = 5.6+0.5 sec. 


Finally, it should be mentioned that in one of the preliminary experiments, 
four precipitations of rubidium perchlorate were carried out at intervals of 
approximately 2 minutes, starting | minute after the end of the irradiation. 
The resulting yttrium samples showed only the 50-min Y° and the 58-day 
Y*. It must therefore be concluded that the 80-sec rubidium isotope observed 
by Hahn and Strassmann (1940) is not connected with either the Rb” or the 
Rb”. A possible explanation is that this activity is the 1.67-min isomer of 
Rb* (INofoed-Hansen and Nielsen 1951). 


(B) Sr 


Two determinations of the Sr half lite have been reported in the literature. 


An early value of 7 minutes was obtained by direct measurement (Lieber 


» 


1939) while a more recent value deduced trom yttrium milking vielded 8.2 


minutes (Knight e¢ al. 1959). The convenient halt lite of Sr favers a direct 


measurement. However, since this nuclide is obtained as a fission product the 


*Since cesium is co-precipitated with the RbCIO,, the precipitate and filtrate were also 
analyzed for Ba®’, The result was the same as for the case of St 
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complexity arising from the presence of other strontium isotopes and the 
ingrowth of their respective yttrium daughters presents a disadvantage. The 
irradiation time and subsequent separation times can be optimized in order to 
reduce the yield of other isotopes. If this is done the only activities which 
contribute to any measurable degree are the ingrowing 10.1-hr Y* daughter 
and the 2.7-hr Sr® — 3.5-hr Y® system. Even then it would appear that any 
measurement involving the decay rate of gross 6 or y radiation would be 
somewhat ambiguous because of the uncertainty in the nature of the back- 
ground. In order to circumvent this background uncertainty it was decided 
to investigate the y rays associated with the decay of Sr® in the hope that more 
favorable conditions could be found. 

A series of measurements were conducted to determine the y rays emitted 
in the decay of Sr%*. This study revealed the presence of strong y rays at 
150, 250, 590, 720, and 900 kev* (see insert in Fig. 2). Since the decay of 
Sr”, Y*%, and Y%® yield no y rays at 590 kev it was decided to use this line 
for a half-life determination. 

The Sr® samples were prepared by irradiating samples containing 100 ug 
of U** for 10 seconds in a pneumatic rabbit where the flux was 5 X10” neutrons 
cm~? sec~!. The irradiated solution was then mixed with 10 mg each of stron- 
tium and barium carrier and, 30 seconds after the irradiation, strontium 
nitrate was precipitated. The strontium was then purified by several barium 
chromate and yttrium hydroxide precipitations, and finally obtained as 
strontium nitrate about 15 minutes after the irradiation. 

The system used for following the decay of the 590-kev line consisted of a 
good quality 3X3 in. Nal(TI) scintillator, a single channel analyzer, and a 
256-channel analyzer. The multichannel analyzer was modified to operate as 
a multichannel scaler. The address was advanced by an external time-scaling 
system, which in this case stepped channels at a rate of 0.4 min/channel. 
The only events accepted for storage were those whose pulse height fell 
within a 4% window centered about the 590-kev line. Typical data obtained 
in this way are shown in Fig. 2. After the analysis was completed (101.6 min) 
a strontium—yttrium separation was performed and the strontium fraction 
again counted for the same. period. This technique is quite useful since it 
gives the shape of the background caused by the Sr®-Y® system. Thus it is 
possible to reconstruct the background which is associated with the original 
measurement. These background components contributed less than 2% to the 
initial counting rate. It was also found that the contribution from ingrowing 
10.1-hr Y® was negligible. Two such experiments were conducted and the 
data were analyzed using a least-squares program. The results obtained were: 
Sample 1: 7.58 
Sample 2: 7.50 

é 


Average: 7.54+.06 min. 


(C) Sr® 
A measurement of the Sr® half life was carried out using the same technique 
as described for Sr®. The y spectrum of Sr® contains a strong line at 1.37 Mev 


*Details on the decay scheme of Sr will be published elsewhere. 
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Fic. 2. A decay curve obtained for Sr. The only events recorded were those whose pulse 
height fell within a 4% window centered about the 590-kev y-ray line. The y-ray spectrum 
associated with the decay of Sr® is shown as an insert. A least-squares analysis of the decay 
curve gave a half life of 7.58+0.06 min. 


(Heath 1957a) which is not greatly influenced by the decay of the ingrowing 
Y* and was therefore ideally suited for this measurement. 

The Sr® source was prepared by precipitating strontium nitrate from an 
irradiated uranyl nitrate solution. Following this several barium chromate 
and yttrium hydroxide precipitations were used to purify the sample. About 
75 minutes after the end of the irradiation strontium nitrate was precipitated 
and counted. The equipment employed was identical with that used for the 
Sr® half-life measurement. 

A least-squares analysis of the data gave a half life of 2.71+0.02 hr. This 
value is in good agreement with published results (Gétte 1941; Herrmann 
and Strassmann 1956) which were done using other techniques. 


(D) Y* AND ¥™ 


Although a number of well-agreeing half lives have been reported for Y® 
(Lieber 1939; Ames et a/. 1952) and Y® (Hahn and Strassmann 1943; Ballou et 
al. 1951; Heath 19576) a remeasurement seemed desirable because of the 
important role these half lives play in the Rb® and Rb® calculations. 

For the half-life determination of Y®, samples were prepared using two 
different techniques. This approach was undertaken in order to check the 
possibility of contamination. The first sample was precipitated as strontium 
nitrate 30 seconds after the end of the irradiation. The strontium nitrate was 
purified by several barium chromate and yttrium hydroxide precipitations. 
The last yttrium hydroxide precipitation was performed 2 hours after the 
irradiation to ensure the complete removal of Y*. The Y® was then allowed 
to grow in from Sr® for 4.3 hours during which time the strontium solution, 
in 9 M HCl, was passed over Dowex 1. After this, yttrium hydroxide was 
successively precipitated a total of three times, each time in the presence of 
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20 mg of strontium holdback carrier. The final precipitate was dissolved and 
reprecipitated as the oxalate. This was then counted after drying. 

The second sample was precipitated as strontium nitrate 30 seconds after 
the end of the irradiation. The strontium nitrate was then dissolved in 9 M 
HCl, passed over Dowex 1, and reprecipitated as strontium nitrate a total of 
three times, the last being 3.5 hours after the irradiation (to remove the Y®%), 
Four hours later the strontium nitrate was dissolved in water, 0.2 mg of iron 
added and ferric hydroxide precipitated. The ferric hydroxide was dissolved 
in 0.3 ml of concentrated HCl and passed over a 0.1 ml Dowex 1 column. After 
it was washed with 0.1 ml of 9 M HCI, the eluant was diluted with H:O to 
about 12 ml and the carrier-free yttrium absorbed on a few milligrams of 
Dowex 50. This resin was then transferred to the top of a Dowex 50 column 
maintained at 65°C. The yttrium was eluted with a 0.2 M solution of a- 
hydroxy-iso-butyric acid adjusted to a pH of 4.2 (Choppin and Silva 1956). 
The elution was monitored with a Geiger counter in order to determine the 
position of the yttrium peak. The carrier-free yttrium, which was eluted in 
less than 15 minutes, was mixed with 20 mg of yttrium carrier, and precipi- 
tated as the oxalate for counting. 

Both samples were counted using a {-scintillation detector. Typical data 
which were accumulated for the half-life measurement are shown in Fiz. 3. A 
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Fic. 3. The decay curve obtained for Y* by gross 8 counting. The least-squares analysis 
was conducted using the data between the two arrows. This yields a half life of 3.53+0.02 hr. 


least-squares analysis which was carried out for both samples gave the following 
results: 
Sample 1: 212.4 min 
Sample 2: 210.8 min 
Average: 211.6+0.8 min 
=3.53+0.02 hr. 


The value obtained is in good agreement with previous measurements. 
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A remeasurement of the Y® half life was also carried out. In this case, the 
objective was to obtain a sample of nearly pure Y* with a minimum amount 
of 3.53-hr Y®. To achieve this condition, a precipitation of fission product 
strontium was done 15 seconds after the irradiation. The strontium was allowed 
to decay for about one half life of the Sr® and the resulting yttrium separated. 
This yttrium is highly enriched in Y* and, at least at the beginning of the 
counting, in Y* and Y®. The yttrium was purified in a manner similar to 
the second Y® sample. 

The counting of the Y® was done in exactly the same manner as described 
for Y®. Figure 4 shows the decay curve obtained. The presence of 20-min 
Y* and 10-min Y® are indicated at the beginning of counting and a small 
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Fic. 4. The decay curve obtained for Y* by gross 8 counting. The arrows enclose the 
region used for the least-squares analysis. The half life was found to be 10.1-0.1 hr. 


contribution of 3.53-hr Y® may possibly exist over the first 100 points (16.7 
the data were analyzed from 1100 


hours). To avoid any influence of the Y® 
to 4100 minutes. A least-squares fit over this region vielded a half life of 


604+1 min = 10.08+0.02 hr. 


(E) CONCLUSION 


The work reported here can be summarized as follows: 


5.3 sec 2.71 hr 7 3.53 hr 

Rb = Se yet nae ogee 
5.6 sec 7.54 min 10.1 hr 

Rb? - pec SOR rey — Zr, 


The techniques described in this paper are, of course, useful to determine 
the half lives of other nuclides. In particular they have been applied to show 
the existence of the hitherto unknown fission products: Rb*, Cs'!, Cs, and 
Ba'. Experiments are in progress to determine their respective half lives. 
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SCATTERING BY A NARROW UNIDIRECTIONALLY 
CONDUCTING INFINITE STRIP! 


S. R. SESHADRE 


ABSTRACT 


The scattering of a plane electromagnetic wave of wave number k by a uni- 
directionally conducting infinite strip of width 2a is investigated, The problem 
is formulated in terms of an integral equation whose solution is obtained by a 
well-known procedure in the form of a series in powers of ka. Expressions for the 
far-zone ‘ields and the first two terms in the series for the total scattering cross 
section are obtained. 


INTRODUCTION 

A unidirectionally conducting screen is one which is perfectly conducting 
in a given direction and is insulating in the perpendicular direction. It is an 
idealization of a screen composed of fine parallel wires. Consider a uni- 
directionally conducting infinite strip of width 2a, illuminated by a plane 
electromagnetic wave of wave number k. It is desired to obtain a solution 
of this scattering problem for the width of the strip much smaller compared 
to wavelength, i.e. ka K< 1. The problem is formulated in terms of an integral 
equation which specifies the current induced on the strip. The integral equation 
is solved for ka K 1 and explicit expressions for the diffraction pattern in the 
far field and the first two terms in the series for the total scattering cross 
section are obtained. The method employed in the solution of the integral 
equation is similar to that used by Bouwkamp (1953) in the solution of the 
corresponding problem of scattering by a perfectly conducting strip. 


FORMULATION OF THE PROBLEM 
A unidirectionally conducting infinite strip occupies the region |x| < a, 
—-o <y< o, and z = 0, where x, y, and z form a right-hand rectangular 
co-ordinate system. It is convenient to introduce also a rotated system of 
co-ordinates £, 7, z where 


= xcosa+y sina, 


ctr 
| 


(1) n = —xsina+ycosa, 

and 0 <a < x/2. The strip is assumed to be conducting in the & direction 
only. When a field E'(r), H'(r) is incident on the strip, there arise the scattered 
fields E’(r), H8(r). The total fields E'(r)+E%(r), H'(r)+H(r) denoted by 
E(r), H(r) respectively satisfy the time harmonic Maxwell's equations 


VXE(r) = :kH(r), 


(2) ae 
—tkE(r) 


I 


Vv <XH(r) 


1Manuscript received August 5, 1960. 
Contribution from the Gordon McKay Laboratory, Harvard University, U.S.A. 
2Present address: Electronics Research and Development Establishment, Bangalore 1, India. 
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in the exterior of the strip. The harmonic time dependence e~ is implied 
for all the field components. The unidirectional conductivity of the strip is 
taken into account by imposing the following boundary conditions on the 


strip, i.e. |x| < a, 


(3) E:(x,y,0) = 0, 
(4) (Hz (x,y)] = He(x,y,0+) —H;(x,9,0-) = 0, 
(5) [E,(x,y)] = 0. 


Further, the scattered fields are required to be outgoing at infinity. In order 
to have a unique solution, it is necessary to guarantee quadratic integrability 
of E(r), H(r) at the edges of the strip and to assume the condition 


(6) [H,(x,y)] = 0, for x = a. 


On account of the unidirectional conductivity of the strip, the current 
induced on the strip and hence the vector potential A(r) has only a & com- 
ponent. It therefore follows that H;*(r) = 0. It is evident (see Karp 1957) 
that when E;'(r) = 0, the incident field is undisturbed by the strip. Hence 
it is possible to decompose the total field into two components: one component 
has the magnetic vector and the other the electric vector transverse to the 
direction in which the strip is conducting. It is obvious that for the second 
component of the incident wave, the strip has no effect and therefore it is 
enough to consider only that component of the incident wave which has its 
magnetic vector in y-z plane. The incident wave may therefore be taken as 


E'(r) = [(¢? +n2)t— bqn— pnz)e*Petantn2) | 


7 
( ) H'(r) a [nh — gle @etantne) | 
where 
pb = lcosa+m sina, 
(8) 
gq = —lsina+m cosa, 


and /, m, n are the direction cosines of the direction of propagation of the 
incident wave. 

The scattered fields are derived from A(r) with the help of (2) and the 
relation 
(9) H®(r) = VXA(r). 


Since H;8(r) = 0, condition (4) is satisfied. Further it follows from considera- 
tions of symmetry with respect to the plane z = 0 that 


sin H®)(2) = —H®,(—2), 


>(s) >(s) 
Ez ,(z) = Ex,(—2). 

As a consequence of (10), condition (5) is satisfied. The vector potential is 
obtained from the current distribution on the strip by the relation 
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eit Fel 
(11) A;(r) -f ¥ J (x’ WZ ia, ‘ace’ dy’ 


where r, @ denote the position vector of a point in space and in the x-y 
plane respectively. It is obvious from the geometry of the problem that all 
the field components depend on y only through the phase factor e*””. Hence 
(11) becomes 


2) dla) = Fe" fe HP bo V8 FET ae 


where 


= ky/(1—m*) = km. 
Together with (1), (2), and (9), (12) gives 


(14) Et(x,z) = “7 La + (cos 4+ ibm sin a) | 


"Ie HS [Row (x—x )? +2" dx? 
With the help of (3) and (7), (14) gives the following integral equation 
(15) | &4+(cos a 4 +ikm sin a) If J (x')H}” [ko|x—x’ | ]dx’ 


=4k(q’+n’)e*™ 


By solving the differential equation (15) for the integral it is easily obtained 
that 


(16) J kJ (x’)HS” (kmi\x—x"'|)dx’ = 4e“* +Ce™? + De *°* 
where 

(17) Cc, = seca—m tana, 

(18) Co = seca+m tana. 


The two arbitrary constants of integration C and D are determined in terms 
of the unknown nee by consideration of the integral and its x derivative. 
For 0 <a < 7/2, J(4a) = 0 by virtue of (6). Hence (16) can be reduced to 





tees ikcgr 
(19) J RI (x! fir iy fe Sebo emi) 
, theyz ,~ tkegz d i 
Hata xe ) pe Hs” (km,|x p | dx 


= intr (I+c2) ,ikeyz (l—c1) ~thear 
= 4¢ ates vo 4e ; 
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SOLUTION OF THE INTEGRAL EQUATION FOR ka <1 


With the substitutions x = acosu, x’ = acosv, e = ka, (19) becomes 


(20) J 2F,(0)| H8(emslcos u— COS v|) te t€c,cos 4 cie7 teen008 “\ 


H§”(em,|cos v}) 








‘ : i d 
—__—. ¢ piteicos u__ ,—tecgcos w — H§ (em,|cos »|) sin v dy 


~ de(¢1-+¢2)sin v dv 


= Ae ttt cos u cae Mites) 4, seosene © (l=c;) C1) 4 eo teeacos u 


CitCe 24a." 
where 0 < u < 7 and 
(21) eF\(v) = J(acosv). 
For small ¢, it can be shown that 

9 
(22) H$” (em,|cos u—cos v|) = = (got+gimie’+gomie'+ ...) 
where 
= p + log|2(cos u—cosv)| = p—2 > ne , 

(23) — 


(-1)" ae iE = 1 <= sone coma 
Zn = 2 nly? (cos cosv)™”"| p—1 ao [— 


1 


for nm > 1 and 


yem, 1. 
(24) p = log —<— Tt, 
log y = 0.5772157 is Euler’s constant, 
' 21 

(25) 5” (em,|cos v|) = ei [goot+gomie +goomie'+ ...] 
where 

go = p z (—1)"cos 2mv 

0 = p— ———, 
: i m 

wi a oo ) | oe ae ¥ (tentne) 

bm = 3G} “7 2g m : 
and 

dy) 2% 22 feos 

(27) yee (em,|cos v|) = = (hothimie +homie + .. -) 
where 


ho = 2 =. (—1)"sin 2mv, 
1 








i 
6. 
i 
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(—1)"*"n . 1 1 SG (—1)"cos 2mv 
(28) h, = oT aly 32 (cos v)”"" sin v| p—1 pal 2, ant ee 


(—1)"(cos v) 


—3-TGnl) > (—1)" sin 2mv. 
A solution is assumed in the form 


(29) Fy(v) = fothietfethetfett+... 


where f,, depend only on v, /, m, and log e. The above expansions are inserted 
in (20), the exponentials are expanded into their power series, and the co- 
efficients of like powers of e on both sides of (20) are equated. The result is 
the following infinite set of integral equations which enable the determination 
of f,’s iteratively. 


r 
. - . 2 
J dv sin v foKo = iBycos u, 
0 
ir 7 
: m 3 ° se 
f dv sinv fiKo = —6, cos u-f dv sin v foKy, 
0 0 
T r or ° 
. Fs ; 4 ‘ - ‘ - 
f dv sin v f2Kyo = —iB2cos u -{f dv sin ofiKi- | dv sin v foKe, 
0 0 0 
r r z= 
‘ . 5 ° - ° aca 
(30) f dv sin v fsKo = 8;cos u-f dv sin vfeKi-f dv sin v f:Ke 
0 0 0 
r 
— f dv sin v foKs, 
v7 


r r r 

« a ' 6 : a ‘ sa naa 

f dv sin v fxKo = 18, cos u-f dy sin vfsKi-f dv sin v foKe 
0 0 0 


-f dv sin bfiKs—f dv sinv foK,... 
0 0 


COs u 


Ky = ota 
° & g sin v 





Soho. 
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ll 

| 
=. 
~ 
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> 
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cos’ *u 


= Solty it dsho, 





(31) Kz = gimi—gormi+y72 cos'u goo— 


cos" 
K3 = 173 cos *u Zoo—1 si 








Bshymi-+i « Ssho, 


4 4 2 2 4 
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and 










a ct —(-c2)"*" 
, (n+1)!(ci+e2) ’ 
p) gw Saal} 





n!(Ci+¢2) ; 


8, = ~ | etc (an ea da) | , 

(n+2)! Citce Citce 
On account of (6) f,(v) = 0 for v = 0, w. Hence f, may be expressed in a 
Fourier sine series and the coefficients of the various terms in the series can 
be determined with the help of the integral equations (30). It is found that 
the solution of (30) is readily obtained in the form 








fon = Dod, sin (27+1) v 


r=) 
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TOTAL SCATTERING COEFFICIENT 
The total scattering cross section o per unit length of the strip is the total 

scattered power per unit length of the strip divided by the incident power 
flow through unit area normal to the direction of the incident wave. The 
ratio of o to the area 2a per unit length of the strip is a dimensionless quantity 
called the total scattering coefficient ¢. By integrating the real part of 
V-E* XH throughout the volume bounded by the two sides of the strip and 
a cvlindrical surface at infinity, it can be easily shown that 


, 


l ao — ik i , “— 
(35) +S Re | eo J (x ) dx’. 


=(l —@ 
With the substitution x’ acos?, € ka in (35), and » relations (21 


(29), and (33), it can be easily derived that the first few terms in the power 


series of ¢ are given as 
(36) 


where 


\ } 


Che substitution of (84) and (387) in (86) together with the 


‘ 


vives the first two terms in the total scattering coefhcient as 
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where 






d = log yem,. 







For the case when the direction of propagation of the incident wave lies in 
a plane perpendicular to the plane of the infinite strip, i.e. m = 0 and when 
a = 0, the total scattering coefficient becomes 


re | Be 8, 4 ‘| 
(39) t= 6 (1-1) 1+76 ates ‘ 


For a = 0, the strip is conducting in the direction perpendicular to its 
edges and the present problem is easily seen to be identical with that of 
scattering by a perfectly conducting infinite strip of a plane wave which has 
no component of the magnetic field in the direction perpendicular to the 
edges of the strip. This problem has been extensively treated in the literature 
(for example, Millar 1960) for the particular case m = 0, which corresponds 
to the propagation vector of the incident wave lying in a plane perpendicular 
to the plane of the strip. The first two terms in the total scattering coefficient 
for this special case are in complete agreement with (39). 



















SCATTERED FIELDS IN THE FAR ZONE 

It is desired to give the expressions for the scattered fields in the far zone. 
Setting x = psin#, z = pcos@ in (12) and taking only the leading term in 
the asymptotic expression of the Hankel function for large p, it is obtained 
that 


(40) A(r) = 12 e i. gene z T 2n+1 oo 


2kmp 0 














where 72,4; is the same as f2,;1, with / replaced by m, sin 6. With the help 
of (1), (2), (9), and (40) the far-zone scattered fields for z > 0 are obtained 


as 





H;§(r) = Q, 






H,'*(r) = —m, cos 6-4, 






H,5(r) = (mcosa—my,sin @sina)A, 


(41) a . : 
E(r) = —[1—(m sin a+my,sin 6 cos a)*|A, 









ES(r) = (m cosa—m,sin 6 sina) (m sina+m,sin 6 cosa)A, 








E,5(r) = m, cos 0(m sin a+m, sin 6 cos a) A 








where 


1 ET ils oil ~ 

‘ (km p—2 /4) 7 2n+2 
42) A=- ——— e™™ o' Te 1 € 3 
( 4 2kmyp » - 


n=0 
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SCALAR DIFFRACTION BY A PROLATE SPHEROID 
AT LOW FREQUENCIES! 
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ABSTRACT 


For the scalar problem of the diffraction of a plane wave by a prolate spheroid 
the exact solution is known, and by expanding this in ascending powers of ka, 
where k is the wave number and 2a is the interfocal distance, the Rayleigh 
series for both the “‘soft‘t and ‘“‘hard‘‘ bodies are obtained up to and including 
terms in (ka)*. The corresponding results for an oblate spheroid can be deduced 
by a trivial change of parameters. Some particular cases are examined. 










1. INTRODUCTION 

It is well known that the scalar problem of diffraction by a perfectly con- 
ducting prolate or oblate spheroid can be solved exactly. The solution is 
obtained as an infinite series of spheroidal functions with specified coefficients, 
but because of the difficulties associated with the computation of these func- 
tions very few calculations have been carried out based on this solution. This 
is particularly true at the shorter wavelengths since the slow convergence of 
the series requires a large number of terms to be considered. 

At the other end of the frequency spectrum the task is more simple, and 
it is a relatively straightforward matter either to compute the exact solution 
or to deduce a sequence of terms in the Rayleigh series. Rather surprisingly, 
the Rayleigh series for the scalar problem seems to have been neglected, and 
apart from the work of Sleator (1957) for the ‘hard’ body, a search of the 
literature has failed to produce anything beyond the first term. Although the 
Rayleigh series has a limited radius of convergence, a significant number of 
terms are required in order to obtain an insight into the way in which the 
scattered field is generated at low frequencies, and in the present paper the 
terms are calculated up to and including (ka)®, where k is the wave number 





















and 2a is the interfocal distance. 








2. THE EXACT SOLUTION 





The solutions for prolate and oblate spheroids can be deduced from one 
another by a trivial change of parameters, and for this reason the analysis 
will be given for only the first type of body. Consider, therefore, a prolate 
spheroid which is defined in terms of the prolate spheroidal co-ordinates 
(t,n,¢) by the equation & = &. Incident upon the body is a plane wave 
travelling in the (x, z) plane of a Cartesian co-ordinate system (x, y, 3) where 


x = a[(1—n?)(€2—1)]? cos ¢, 
y = a[(1—n?)(¢2—1)]? sin ¢, 










z= aré, 






‘Manuscript received August 12, 1960. 
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and if the direction of incidence makes an angle ¢ with the positive z axis, 
the field can be written as 







vi _— etk(z sin [+2 cos f) 






‘ 





A time factor e~* is employed and suppressed. 
The incident field can also be expanded in an infinite series of spheroidal 
functions in the form* 







Ms 


(1) v' = x Ymn (c, o) Von (c; nN; & ) 


0 


2 
I 





where c = ka and 









(2) Ymn(C, €) = 2(2—b0m)i"Ninn Smn(C, cos £), 
(3) VP, = Sinn(C, n)R2(c, £) cos me. 


Sinn(e,m) and RY (c, &) are the angular and radial functions respectively, 
and it will be observed that in equation (1) the radial function is of the first 
kind. For the scattered field V* a similar expansion is assumed, and this is 




















(4) Viw 2 2 Athntts6s @ 


m=() n=m 


where the amplitudes 4,,,, are to be determined from the boundary condition 






at E — Eo. 
If the body is ‘soft’ so that the boundary condition at & = & is 






V'+V* =0 










(Dirichlet condition), the orthogonality of the cos m@ and the angular functions 


gives immediately 










Rn (C, £0) 
5 Aye = —Vmn\C, )- ea 
(5) Vault, RC ) 


3 $0 





In the ‘hard’ case for which the Neumann boundary condition applies (vanish- 
ing normal derivative of V'+V®* at & = & ), the corresponding amplitudes 
differ from equation (5) in that each radial function is replaced by its first 
derivative with respect to £. It follows that by differentiation it is possible 
to deduce the results for the Neumann problem from those for the Dirichlet, 
and we shall therefore restrict attention to the latter problem. It is noted in 
passing that this similarity enables us to deduce the Rayleigh series for the 
hard body from the Rayleigh series for the soft body by differentiating each 
radial function (and hence, each Legendre function in terms of which the 
RD) (c, &) are expanded) with respect to £, and although this might suggest 
that the two series are similar in form, we shall see later that the differentiation 












serves to remove some of the powers of c. 






*The notation is that of Flammer (1957), and the reader is referred to this book for the 
definitions of the symbols and functions here used. 
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For the sake of simplicity we shall now particularize the problem to the 
case of a field which is incident in the direction of the negative z axis, so that 


V! = eg tke 


(corresponding to ¢ = 7m). Since 


lim Sya(c, cos ¢) = 0, n ~ 0, 
for 


Ymn is zero unless m = 0, and the summation over m in equations (1) and (4) 
produces only the first term. Hence 


@~ 


(6) V° = D> AcwSon(c, 0) Row (c, €) 


n=0 
with 


(7 a: ce: aegis ) Rin (6s £0) 
() Alm = Yon\C, R® 5)" 


and both V' and V* are now independent of ¢. 
In the far field c& ~ RR, where R is the distance from the center of the 


spheroid. For large values of cé, 


(3 . 1 rere ’ 
Rv (ce, &) ~w—e ies 4 +1) I 
cé 


and if this is inserted into equation (6) the coefficient of e“*/RR is the far- 
field amplitude and is 

a Son(c, —1) Rin (c, 0) 

fn) = 2 D> = — =F go 


(8) = — Son(C, 7). 
n=O Non Ron (c, £o) on(¢ 1) 


For sufficiently small values of c, f(7) can be expanded in a convergent series 
of positive increasing powers of ¢ with coefficients which are functions of & 
and 7. This is the Rayleigh series for the body, and it is the purpose of the 


paper to calculate the first few terms. 


3. THE RAYLEIGH SERIES 
The angular spheroidal function So,(c, 7) is defined as 


2 


Sin(e,n) = 30’ dM (c)P,(n) 
r=(),1 
where /,(y) is the Legendre polynomial of order r and argument n, and the 
summation extends over even or odd values of r according as n is even or 
odd. This is denoted by a prime attached to the summation sign. The radial 
functions Ri’ (c, —) are similarly defined in terms of Legendre functions of 
the first and/or second kind of argument £, and the coefficients are either 
the same d”, or can be determined from them. The d™ themselves are given 
by a recurrence relation coupled with a normalizing equation, and since each 


coefficient can be expanded in a series of powers of ¢, it is only necessary to 















(9) 








where 






















first derivative, 


Legendre function identities in simplifving the 


J@) = -¢ x an (&) GS 





405 


t 6 
a2(ko) = ; c*+0(« 


above expressions. 
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*Unless otherwise shown, the Legendre functions have argument £ 
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Qo 







insert the appropriate expansions to obtain the required Rayleigh series for 
f(n). Moreover, a trivial calculation shows that the first three terms in equation 
(8) are sufficient to specify the terms in the Rayleigh series up to and including 
c®, and this is as far as the analysis will be carried. 
Using the expansions given in the Appendix it now follows that* 








) 














and a,(&) = O(c*), 2 > 3. The corresponding expansion for the hard body 
can be obtained by replacing each Legendre function of argument & 
and to preserve this feature no attempt was made to use 
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To complete the expansion of f(y) it is necessary to expand the angular 
spheroidal functions So,(c, 7), and to achieve the simplest possible form for 
the final results we shall now sacrifice the symmetry property referred to 
above. Taking first the solution f(7) = f;(n) for the soft body, and using the 
Legendre function identities 


ee 
P2D, ic Fas 


‘ 1 
Pe 1 = gPrtghs 


we have 


(13) faln) = 6D ta(n)(—ie)" 


where 
> ae 0 
a oe, 
= 0,7 0) 
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+ r o(n), 
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18 Qo 3Qo | 225 Oo) 9\O/ ~ 2025 
For the hard body the expansion for the far-field amplitude f(y) fn (n) 
is of somewhat simpler form. Since Po’(&o) = 0, where the prime denotes 


differentiation with respect to &, the cae for do(&) deduced from 


equation (10) contains only terms in c* and c*, and if 
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The differences between the Rayleigh series for the hard and soft bodies 
are now apparent, and whereas f,(y) contains all powers of ¢ starting with 
the first, fyx(m) begins with c* and has no term in ct. In this respect the expan- 
sion for fy(m) is similar to the Rayleigh series for the vector problem (see, 
for example, Stevenson 1953), and the agreement is further strengthened by 
the fact that the first imaginary term in fy(n) is provided by v5 (and is there- 
fore of order c*). A comparison between the actual coefficients will be given 


+ & 
. ae 


“J 
qo 


at a later date. 
For the oblate spheroid the appropriate expansions can be obtained from 


equations (13) and (14) by replacing ¢ by te and & by Z&. 


tf, PARTICULAR CASES 

We shall conclude by examining the form taken by the above expansions 
when the spheroid degenerates into either a sphere or a disk. 

If & > © and a — 0 in such a way that ago tends to a finite limit p/k, the 
spheroid reduces to a sphere of this radius. The Legendre functions in the 
expressions for the “,() and v(m) can now be replaced by the leading terms 
in their asymptotic developments for large £, and in the limit cf = p the 
Rayleigh series for the soft and hard spheres are found to be 


» » 


hives » j ‘ 4 2 l . § a i a .-§ sd 
(15) I.(n) = “4 \! to; T3 lp ri? 45 ip Po(n) 


l . | ' l 6 
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o> “ 





1638 CANADIAN JOURNAL OF PHYSICS. VOL. 38, 1960 


(16) Ko Aaetde aero 


+50" ee oe. 4; ; tp tp, i(n) — = 'Px(n) +O(6 »|, 


respectively. These are in agreement with the results obtained by expanding 
the Bessel function ratios in the known scalar solutions for a sphere. 

The second case to be considered is the disk, and for this purpose it is 
necessary to go over to oblate spheroidal co-ordinates. Replacing c by —ic 
and & by 2 and then letting &) > 0, the (oblate) spheroid reduces to a disk 
of radius a, and the Rayleigh series for this body can therefore be deduced 
from equations (13) and (14). Putting & = 0 in these equations with —ic 
instead of c, we have 


- \_ eff, _2ic (4 2) ee i) 
(17) f£(@) = —-* -5)+== (4-5 
(16_ 16 ee a 
te e. Or” TT. a \n Ort 2 2025 g Po(n) 
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(18) fx(m) = ae a ytgee +2 * tp, 1(n) toe CP s(n) +06 5]. 


In equation (17) the terms up to and including c® can be compared with 
those obtained by Bouwkamp (1954) for a circular aperture in a hard screen 
and if account is taken of the difference in sign implied by Babinet’s principle, 
the expansions are seen to be in agreement. 
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The radial function of the first kind, R{”(c, &) is defined as 


RY a 
On (C, f) = RiP (c) Son(c, £), 


len TE a (n even) 
(n odd) 


Hence 


] 
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The radial function of the third kind, R{* (c, &), is 
Rin (¢, £) = Rin (c, £) +iRon (c, ) 
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These expansions are sufficient to give the far-field amplitude correct to c°. 





THE HYDROGEN ION EFFECT IN WHISTLER DISPERSION! 


R. E. BARRINGTON AND T. NISHIZAKI 


ABSTRACT 


Four low-altitude whistlers have been carefully analyzed by a specially 
developed filtering technique. In each case the low-frequency portions of the 
resulting dispersion curves show similar departures from the simple tft = const. 
law developed by Eckersley. The departures are found to be of the form and 
magnitude predicted by Storey on the assumption that the propagation is partly 
supported by the light hydrogen ions of the exosphere. If it is assumed that the 
observed deviations are due only to hydrogen ions, and that above some transi- 
tional level the positive ions are mostly hydrogen, an estimate of 1000 km for the 
height of the transitional level is obtained. 


1. INTRODUCTION 


Whistlers are electromagnetic radiations of natural origin found in the 
audio-frequency band of the radio spectrum where they normally appear as 
quasi-musical tones of descending pitch. Their interpretation as provided by 
Storey (1953) implies a much greater extent to the earth’s outer ionosphere, 
its exosphere, than was believed a decade ago, and presents the problem of 
the origin, composition, and abundance of the ionization thus revealed. 
Whistlers have provided considerable information on the exospheric electron 
density, but an understanding of its origin requires a knowledge of the positive 
ions and their distribution. Some authors (for example, Chapman 1957) have 
suggested that solar protons, injected into the earth’s magnetic field, provide 
the bulk of the positive ions. Others (for example, Johnson 1960) consider the 
exosphere to be populated by hydrogen which has diffused from the lower 


ionosphere, has been ionized by charge exchange or some other process, and 


has been trapped in the earth’s magnetic field. These conflicting ideas only 
emphasize the need of measurements of the positive-ion constituents of the 
exosphere. 

Storey (1956) has shown that the relationship between the instantaneous 
frequency f of a whistler, and the time ¢ taken by this frequency to travel from 
source to observer, should provide information about the presence of hydrogen 
ions along the whistler path. This relationship may be conveniently studied 
by determining the variation of the whistler ‘dispersion’, which will be 
defined here as ¢f?, with frequency. To a first approximation the dispersion is 
constant although at frequencies approaching the electron gyro-frequency it 
increases and gives rise to the well-known ‘nose’ effect in whistlers. There is 
a similar but smaller increase in dispersion for frequencies near the gyro- 
frequency of any abundant positive ion of the medium, and it is this variation 
which is of particular interest in the present paper. Of all the positive ions 
likely to be present in appreciable quantities in the exosphere, only hydrogen 

‘Manuscript received July 15, 1960. 
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has a gyro-frequency of sufficient magnitude to produce a measurable effect 
on whistler dispersion, in the frequency range of observation. Thus whistler 
dispersion is a potential source of information about the protons of the high 
ionosphere. 

In the work reported here, the dispersion of four selected whistlers observed 
in Japan was measured. The analysis was carried out with high accuracy 
since the maximum effect due to the hydrogen ions was 2 or 3% of the total 
dispersion at most. The method by which this was achieved is outlined in 
Section 2 and the results are presented in Section 3. A model exosphere is 
employed in Section 4 to derive dispersion curves, which are then compared 
with the observed results for the two best whistlers. In the final section the 
implications of this analysis in terms of the hydrogen ion component of the 
exosphere are discussed. 


2. METHOD OF DISPERSION MEASUREMENT 

Whistlers from Wakkanai, geomagnetic latitude 35.3° N., and Toyokawa, 
geomagnetic latitude 24.5° N., were chosen for analysis primarily because the 
protons of the exosphere are expected to account for a greater percentage of 
the dispersion at these stations than at higher latitudes. It has also been found 
that the whistlers recorded at these stations have a greater purity of tone 
than those observed at higher latitudes. Since the analysis is slow and tedious, 
only a few of the best samples were used. The examples chosen were well 
defined over a frequency range greater than 5 kc/s and extended down to at 
least 1 ke/s. 

A preliminary study of the whistler recordings was made with a sonograph. 
This provided the approximate rate of change of frequency in various parts 
of the whistler spectrum. It also gave an indication of those individual whistlers 
which offered some possibility for the identification of their initiating spheric. 
In a similar investigation Iwai and Outsu (1959) used a sonograph to measure 
the relative arrival times of the frequency components of the whistler. This 
complicates the analysis. The filtering technique used in the sonograph is such 
that its response to a gliding tone is asymmetric in time, reaching its maximum 
somewhat after resonance. This time delay is a function of filter bandwidth 
and rate of change of frequency in the gliding tone. The net result is to introduce 
an error into the measurement of the whistler dispersion which has the same 
variation with frequency as the hydrogen ion effect. 

Storey and Grierson (1958) analyzed the response of various filters to a 
gliding tone, and developed a filtering technique which gives a response 
symmetrical in time and centered on the instant of resonance. This obviates 
the need for any correction to the time of maximum filter response, and was 
used here in all cases to measure the relative time delays of the frequency 
components of the whistler. The method utilizes a parallel tuned narrow-band 
filter and time-reversal tape recording. The filter output is rectified and time 
marker pulses are added to it at intervals of 10 msec by a combined adding 
and switching circuit (Grierson 1958) which suppresses the filter output during 


the marker pulses. When displayed on an oscilloscope, displacement of the 
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trace in one direction from a fixed base line shows the envelope of the filter 
response, and displacement in the opposite direction, the timing marks. The 
waveform appearing on the oscilloscope was recorded on moving 35-mm film. 
Figure 1 shows a sample of the results obtained in this manner. 

The maximum of the filter response can usually be determined with an 
accuracy of about +1 millisecond. The resonant frequency of the filter was 
measured with a signal generator and counter, and is accurate to +1 cycle/sec. 
The bandwidth of the filter was set for optimum signal-to-noise by the criterion 


4eW? = |df/dt| 


where W = bandwidth of filter between 3-db points, 
and df/dt = rate of change of frequency in the gliding tone. 


To convert the relative time delays of the whistler frequency components 
to dispersion, the time of the initiating spheric must be known. To determine 
this, the whistler and a few seconds of the record preceding it was played 
through a 10 kc/s bandwidth filter, and displayed as before. This procedure 
makes the spherics stand out above the noise as sharp spikes. The approximate 
time of the initiating spheric was determined by extrapolating linearly from 
the mid-frequencies of a plot of relative time delays against f-*. Allowance 
was made for the time delay introduced by the propagation of the spheric to 
the observer on the assumption that the spheric occurred at the conjugate 
point. In three of the four cases studied here there was only one spheric within 
10 msec of the approximate time origin determined in this manner. This 
spheric was taken to be the initiating spheric and all times were measured from 
it in calculating the dispersion. The error involved in measuring the time of 
this spheric is less than 1 msec. In the case of the fourth whistler there were 
three spherics fairly close to the approximate time origin. Hence in computing 
the dispersion the approximate time was used, and the resulting dispersion 


probably has a time error. 


3. EXPERIMENTAL RESULTS 

Figures 2-5 show the dispersion of the four whistlers (denoted by numerals 
1-4 respectively) which were analyzed, as a function of frequency. Comparison 
of the event observed at Toyokawa, for which the time origin is not accurately 
known, with the other three examples seems to indicate that the ion effect 
is about the same size at both stations. This may be due to a small time error 
associated with the Toyokawa results. The important feature of these curves 
is their similarity for whistlers observed in different seasons, and at different 
times of day. This indicates that the observed variation of dispersion with 
frequency is due to some constant properties of the path. It has been established 
that the electronic constituent of the exosphere can only cause the dispersion 
to increase with increasing frequency. Thus it seems likely that the increase of 
dispersion at low frequencies which is observed in all cases studied, and is of 
the form predicted by Storey from considerations of the effects of protons on 


the propagation of whistlers, is indeed a hydrogen ion effect. 





I 


f=3100 


Fons 


ae: <E 
| ~ f= 2700 : 


Fic. 1. Filter response to a whistler at various frequencies as a function of time. The large 
time markers are at 50-msec intervals. 
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Fic. 2. Dispersion of whistler No. 1 observed at Wakkanai, geomagnetic latitude 35.3, 
March 9, 1958, 2-5 hr JST. 

Fic. 3. Dispersion of whistler No. 2 observed at Wakkanai, geomagnetic latitude 35.3, 
Dec. 24, 1957, 1600 JST. 
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Fic. 4. Dispersion of whistler No. 3 observed at Wakkanai, geomagnetic latitude 35.3, 
March 9, 1958, 2-5 hr JST. 

Fic. 5. Dispersion of whistler No. 4 observed at Toyokawa, geomagnetic latitude 24.5, 
Feb. 23, 1956, 1700 JST. 
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The scatter of the experimental points appears to be due largely to a residual 
multipath effect, rather than to errors introduced by the method of analysis. 
In some cases the scatter represents errors in time measurements of 5 msec. 
This is much greater than the error involved in estimating the time of maximum 
of the filter response. It was also noticed that occasionally the filter response 
was asymmetric or had a two-humped shape. This is suggestive of interference 
effects due to signals which have propagated by slightly different paths. 
Further support for such a conclusion comes from the fact that sometimes 
several of the experimental points show a systematic trend of a finer scale 
than the over-all frequency dependence of the dispersion. Such a multipath 
effect is not unexpected, since for whistlers observed at higher geomagnetic 
latitudes the errors introduced by this effect are so large that they mask the 
ion effect entirely. 


4. COMPARISON OF EXPERIMENTAL AND CALCULATED DISPERSIONS 


It is possible that ducts or inhomogeneities of the medium might have effects 
on whistler dispersion similar to those predicted for protons. In order to 
provide additional support for the view that the observed effects are due to 
protons, and also to provide some indication of the proton distribution and 
abundance, a plausible model of the exosphere was used to calculate the 
variation of disperson with frequency. Such curves are not found to be critically 
dependent on the electron-density model chosen, and do fit the observed 
dispersion curves surprisingly well. 

In making such a test, several plausible assumptions which cannot be 
entirely justified are necessary. It is assumed that the whistler propagated 
along the line of force on which it was observed. The gyro-frequencies along 
this path are considered to be adequately described by the centered-dipole 
approximation. The longitudinal refractive index is used to calculate the 
dispersion. This is a good approximation if propagation is by a magneto-ionic 
duct as suggested by Helliwell (1960). The purity of tone of the whistlers 
studied here supports the idea of propagation by such a duct. 

The dispersion of a whistler when the motion of the positive ions is taken 
into account can be approximated by combining the asymptotic expressions 
given by Storey (1956) to form the equation 


oe ee 14+2%fa/f_ V4, 
OC path 2NG-f/ ny [1+ a /f a, ~~ 


where fo = plasma goa 


(1) 


Su electron gyro-frequency, 
fm = proton gyro-frequency, 
f = wave frequency. 


This applies only for a medium in which the numbers of protons and electrons 
are equal. If the medium contains more than one type of positive ion, of which 
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only the protons have a gyro-frequency of sufficient magnitude to affect the 
observed dispersion, equation (1) must be replaced by, 


_ ad al 


72 


af ds 
1 


2c path Su 


(2) ii? ~ 


C—fifa 


i 
he 
where € = ratio of number of protons to number of electrons. 

The first term within the brackets of equation (2) is due to the electrons 
along the whistler path, and may be expanded as a binomial series in f/fy. 
Since the whistlers examined in the present study originate in low geomagnetic 
latitudes, the ratio f/f” was always much less than unity along the whole 
length of the path even for the highest frequencies observed. The terms in the 
power series above that in f/f are negligible and the dispersion due to the 
electrons may be written 


ae 

(3) Da = es fe, as+|3 J Jo, as|p = Dot Dif 
<C J path tu 4c Ju 

say. The first term here is independent of frequency and approximately equal 

to the mean whistler dispersion. The second term is small and shows that in 

the frequency interval of interest, the electrons of the medium cause the 

dispersion to increase linearly with frequency. 

The frequency dependence of the other terms in equation (2), due to the 
protons on the path of propagation, is not as simply expressed as that due to 
the electrons. Over most of the whistler path, €fi/(fia + f) is small compared 
with unity. Under these circumstances the contribution of the protons to the 
observed dispersion is given approximately by 


1 . fo j faltmfiey® 
36 Jia? Vfl So 


The bracketed part of the integrand of this expression is plotted as a function 


(4) 


of fia/f in Fig. 6. Since this curve passes through the origin and is a single- 
valued function of fi/f with a single broad maximum in the interval 
0<fim/f<.6 it can be approximated in this range by a function of the form 


(5) Df? + Duif”. 


Thus the frequency dependence of the dispersion of a whistler may then be 


approximated by 
(6) if’? ~ Df? + Dif + Dy + Dif. 
In this expression Dy and D,; depend on the electron distribution along the 


path while D 
values of these constants for whistlers | and 2 were determined by a least- 


» and D-; depend on the proton and electron distribution. The 
squares fit of the experimental points to a function of the form (6). In such a 
fit Do is determined with the greatest accuracy, while D_s and D_, have the 


great est errors. 
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Fic. 6.. Plot of 3 fin(f—fin)/(fia +f)? as a function of fin/f. 


* The ratio of Do/D, has been shown by Storey (1957) to provide a limit to 
the range of latitudes over which the whistler might have entered the iono- 
sphere, and also information about the electron distribution along the path. 
It was found to be much more sensitive to the path of propagation than to the 
electron distribution so it was not employed as a check on the latter. Instead, 
it was used to confirm the assumption that the whistlers studied here actually 


propagated along the line of force on which the observing station is situated. 
In the case of whistlers 1 and 2 the ratio Do/D,; was determined fairly 
accurately, and in both cases showed that the whistler emerged from the 
ionosphere within 2 er 3 degrees of the observing station. 

The electron distribution was determined by using the value of Do to 
normalize models of the exosphere suggested by other workers. Johnson 
(1960), from theoretical considerations, has arrived at an electron density 


distribution given by 
(7) N(h) = No exp [goRo?’m/(Ro + 4) KT] exp [m 0°(R + h)? cos’0/2KT} 
where N(h) = number density at altitude /, 

No normalization constant, 

£0 acceleration of gravity at the earth’s surface, 

Ro = radius of the earth, 

m half the ion mass, 

Q the earth’s rotational velocity, 


6 the latitude of the point under consideration. 
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At the Naval Research Laboratory a model ionosphere has been derived by a 
compilation of sounding, rocket, satellite, and whistler data. These two models 
were adjusted by the choice of an arbitrary factor to give the observed values 
of Dy and are compared in Fig. 7. While the derived distributions differ 
somewhat, the discrepancies are not large compared with the errors of measure- 


—— JOHNSONS DISTRIBUTION 
--- NRL'S DISTRIBUTION 


nN 
°o 
o 
o 


€ 
x 
- 
x 
S 
uw 
= 


103 10* 
ELECTRONS/cm-3 


Fic. 7. Comparison of electron densities derived for Johnson’s model and the Naval 
Research Laboratory model. 


ment encountered elsewhere in this work. (The differences are. somewhat 
larger at greater heights.) In subsequent calculations the Naval Research 
Laboratory's model is used, although Johnson’s model of electron densities 
leads to essentially the same results. 

D_» and D_, contain information about the number and distribution of 
protons along the whistler path, but they are not sufficient to determine the 
distribution. They also have errors of measurement associated with them, as 
well as inaccuracies due to the assumption that equation (2) can be approxi- 
mated by equation (6). Rather than employ these quantities to determine 
something about the proton distribution, an idealized model of the exosphere 
was used. In the absence of any measurements of the proton distribution it 
was assumed that, from the top of the whistler path down to some transition 
level hy, the number of protons and electrons is equal. Below this level the 
number of protons is assumed to be zero. On this basis the approximations of 
equations (2) and (3) can be avoided, and (1) applied with the aid of the 
electron-density distribution described previously to calculate dispersion 
curves for various heights of the transition level. The relationship 


h h, 
1/2 _ tf ae tf ” fo { 1 ie 1+2f m/f_ \ 
(8) f C v0 fa” saad” or fa” (—f/fay* +f /fP? " 
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_Fic. 8. Comparison of experimental points with calculated dispersion curves for whistler 
No. 1. Dispersion curves are shown for transition heights of 750, 1000, and 1250 km, and for 


Johnson’s model of the exosphere. ; . , , 
Fic. 9. Comparison of experimental points with calculated dispersion curves for whistler 


No. 2. Dispersion curves are shown for transition heights of 750, 1000, and 1250 km, and for 
Johnson’s model of the exosphere. 
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with h,, = greatest height of whistler path, was calculated for various values 
of hy the assumed transition height and the experimental points were then 
fitted to the calculated curves as shown in Figs. 8 and 9 for whistlers 1 and 2. 
Also shown on these figures are dispersion curves calculated from Johnson’s 
predicted proton and electron distributions. 


5. DISCUSSION AND CONCLUSIONS 


On comparing the experimental points with the curves calculated on the 
basis of the preceding assumptions, the agreement is found to be quite good, 
especially in the case of the curves calculated for a transition level of 1000 km. 
The fact that the analysis of two whistlers observed at different times and 
seasons leads to the same transition level for the protons, even though the 
values of Dy and hence the eiectroii densities of the exosphere are quite different 
at the times of observation, lends further support to the conclusion that the 
observed behavior of the dispersion at low frequencies is due to protons. 

The proton distribution used in calculating the dispersion is at best a gross 
approximation to the actual distribution. It is useful in demonstrating that the 
ion effects observed under different conditions fit a single model for the proton 
distribution. Johnson’s model of this distribution, which is derived from 
considerations of the escape and ionization of hydrogen diffusing from the 
lower atmosphere, does not fit the observations as well. While observations of 
the ion effect in whistler dispersion do not make it possible at present to derive 
a unique model of the exospheric ionization, they do provide an observational 
test for any model which is suggested by other considerations. Whistler dis- 
persion does show tnat protons form a major constituent of the exosphere, 
down to heights of 1000 km or below. 

More accurate information about the path of propagation and the electron 
densities of the ionosphere could be obtained if the dispersion were measured 
at higher frequencies (Storey 1957). This would increase the accuracy with 
which the dispersion due to the electrons of the medium could be separated 
from that due to the ions. Any further increase in accuracy of measurements 
of the ion dispersion requires a reduction in the multipath effects. Some pro- 
gress in this direction might be made by more careful selection of the whistlers 
used in the analysis. Satellite observations of whistlers offer the possibility of 
a much greater reduction of this effect as well as better information on the 
path of propagation. Probably the largest uncertainty of this work is in the 
location of the initiating spheric. This could be overcome by timing and direc- 
tion-finding equipment near the point conjugate to that at which a satellite 
observed a whistler. 
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ELECTRON —- HYDROGEN ATOM SCATTERING 
IN BORN APPROXIMATION! 


Ta-You Wu 


ABSTRACT 


The elastic (1s—1s) and the inelastic (1s—2s, ls-2p) scattering cross sections in 
the Born approximation at energies of 1, 4, 9, 16 rydbergs have been calculated 
exactly from the closed formulas of the matrix elements for these transitions. Both 
the differential and the total cross sections are given here. 


The theory of the scattering of an electron by a hydrogen atom in the 
Born approximation has been known since 1926, and with exchange, since 
1928. The fact remains, however, that no exact calculations of this approxima- 
tion for even such simple transitions as 1s—2s, 1s-2p have been carried out. 
The only calculation in the literature seems to be the early work of Massey 
and Mohr (1931), in which the “exchange’”’ scattering amplitude is evaluated 
by means of some approximations (namely, the neglect of higher order terms 
in an infinite series and the replacement of Bessel functions in the integrals 
by its asymptotic forms). Some years ago, such matrix elements as those 
occurring in the Born approximation have been evaluated in this laboratory 
in closed analytic forms (Corinaldesi and Trainor 1952). The actual calcula- 
tions have now been carried out and the results are presented below. 

In atomic units (length in Bohr radius ag, and energy in rydbergs), the 
scattering amplitudes f(@), g(@) for “‘direct’’ and ‘“exchange’’ scattering are, 
for the transition k, do — Rn’, dn; 

) = 2 f fein 2 | sei 
(1) f(@) = ins J ; é b4(r2)| - i hal e go(r2) dridro, 

e iJ , 
(2) g(6) = a j fe* oy(r)| -24+ | e" bo(re) dridrs. 
4nJ. re 


[ri—ro| 


The total cross section is given by 


(3) ° an ae 4 JO) +e00)°-+2/0)—eoI* sin 6 dé 


=2r ky f\r@+e@-eoye} sin 0 dé. 


ex 
(A). 1s —1s (ELASTIC SCATTERING) 
1 ] 
f(0) = 2P?~ OP AFP)’ 
OT cite a nll lade bi Vacs ght 
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V(1+P") 
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where ? is the energy of the incident electron in units of rydbergs, and 
P = ksin (6/2). 


From these expressions, one obtains, in the forward direction, 





: _ i a ‘|. 
(B). Ils ~2s (INELASTIC SCATTERING) 
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(C). 1s +2p, Am = 0 (INELASTIC SCATTERING) 
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536>/2 ; 
f(6) = 15360/ [2k—x~/4k?— 1), 


“FEE 


1024 \/2, 8 /2k 
g(0) = a [K Co+32xC3+96C4] 


2k—x/4k?—3 
— [15x?D1+-48«D2+128D3]. 
K 


(D). 1s > 2p (Am = +1) (INELASTIC SCATTERING) 


(0) = 15364/(1—x?) ao V/ (4k? —3) e™, 
(8) 


g(0) = 0. 

The differential cross sections f?+ g? —fg for the 1s—1s, 1s—2s, 1s-2p (Am = 0), 
ls-2p (Am = +1) transitions at energies 1, 4, 9, 16 rydbergs (13.6, 54.4, 
122.4, 217.6 ev respectively) are shown in logarithmic scale in Figs. 1-4. (The 
ordinates in Figs. 1-4 give the differential cross section f?+g?—fg as defined 
in eqs. (1)-(3) in the text, in ag”. k = 1, 2, 3, 4 correspond to incident electron 
energy of 1, 4, 9, 16 rydbergs respectively.) The total cross sections are given 
below. 


Energy 1s—2p 1s-2p 
in rydbergs  1s-Is 1s—2s Am=0 Am= +1 








1 2.104 1.486 1.351 0.0748 
4 0.574 0.0862 0.3981 0.2549 
9 0.256 0.0406 0.2039 0.1945 
16 0.144 0.0236 0.1197 0.1437 


Total cross sections in units of rap?. 


There are very meager data with which these exact Born approximation 
results can be compared. The total cross sections of elastic scattering at 
k? = 1,4 can be compared with the variational calculations of Massey and 
Moiseiwitsch (1951). Those for the inelastic scatterings 1s—2s, 1s-2p can be 
compared with the approximate values obtained by Massey and Mohr 
(1931). The comparison is shown below. 








Born approx’n 
k? energy Born Massey—Mohr Variational 
in approx’n (from their Massey— 
rydbergs here calc. Fig. 4) Moiseiwitsch 
1s—1s .60 10.5 
.80 1.55 
1s—2s 65 ~5.3 
‘ ~2.1 
27 
1s-2p 71* ~6.2 
3.5 ~2.8 
2.83* 





Total cross sections in ap?. The values * are the sums of the cross sections for 
Am = 0, +1. 
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Fic. 1, Elastic scattering 1s — 1s. 
Fic. 2, Inelastic scattering 1s — 2s. 
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Inelastic scattering 1s + 2p, m 
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Fic. 4. Inelastic scattering ls + 2p, m 
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It is seen that at 13 ev, the Born approximation for 1s—ls is not satis- 
factory, as expected, but at 54 ev, it is in quite satisfactory agreement with 
the variational calculations of Massey and Moiseiwitsch, which is the ‘“‘best”’ 
theoretical result available at present. 


The writer appreciates the help of Mrs. M. Helfenstein in the computations 
involved. 
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A SOLUTION OF THE EINSTEIN FIELD EQUATIONS! 


PETER RASTALL 


ABSTRACT 


An exact, cylindrically symmetric, time-dependent solution of the Einstein 
gravitational field equations for empty space is derived. A particular case of 
the solution has singularities only on the axis of symmetry and may represent a 
number of particles in an otherwise empty universe. 


INTRODUCTION 


Static, cylindrically symmetric solutions of the Einstein gravitational field 
equations for empty space were considered by Weyl] (1917) and Levi-Civita 
(1919). They showed that no such solutions exist which are pseudo-Euclidean 
at infinity and which have more than one singularity on the axis of symmetry. 
This must be so if the singularities are interpreted as the positions of particles 
and the field equations are to imply the proper equations of motion. Later 
work on cylindrically symmetric solutions is summarized by Marder (1958, 
1959). 

We shall look for a class of solutions which are cylindrically symmetric but 
not static. We assume a line element of the form 


(1) ds? = h(x®) {u(x?, x2) (dx! +dx2") + v(x!, x2) dx3?} + A(x!, x2) dx”, 
where hk: is not a constant. 
SOLUTION OF THE FIELD EQUATIONS 


When the metric tensor is diagonal, g,, = 0 for  ¥ p, the most convenient 
form for the contracted curvature tensor R,, is 


(2) Re m (1/4) = 12(€rrw/ Bre) a t2Gee0/ Br). 


THE 


- (Zax ul Gee) a_i + (ee, shee) 


— (gix,r/Err8er) + > tel teGieltala ‘ 


on 


) J. ) 
(3) Res = (1/4) bs 2 Err /B er) 0 Gee.n/ Bex) (Ber x/ Ber) 


THT ,p 


— (Bo0,.7/Bo0) (Srr.0/Brr) + (See (Ber) (Brr.e/Brr) § , 


In these two equations the summation convention has been suspended and 
w ¥ p. 


‘Manuscript received June 10, 1960. 
Contribution from the Department of Physics, University of British Columbia, Vancouver 8, 
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Substitute the metric tensor of (1) into (3). The field equations Rip = 0 

and Rey = 0 imply 
Ayho/Ah = 0, Aoho/AA = 0, 

where to save writing commas we put fy = dh/dx®, \y = 0\/dx', etc. Since 
hy ¥ 0, except perhaps at certain instanis, it follows that \ is a constant, and 
we may take A = —1. 

One can separate variables in the remaining equations R,, = 0 and find 
for h the two equations 


Zhhoo — ho? = 0, 2hhoo os ho? = 2 Ah. 


A is a separation constant. The solution, provided that the origin of time be 
suitably chosen, is 


(4) h = (A/4)x”. 


The assumption that / is not constant implies A # 0. 
uw and » must satisfy the four equations 


(5) 2(ui1 + wee)/m — 2 (ur? + wo?) /w? = WA, 

(6) (viy — vo2)/v — (vy? — vo?) /2v? — (uni — weve)/uv = O, 
(7) 2y12/v — (uve + wovr1)/mv — vive/v? = 0, 

(8) (vin + v22)/v — (ve + 2?) /2v? = pA. 


Equation (5) involves only u and is of a type which has been discussed by the 
mathematicians. It is best to introduce new dependent variables u and v by 


io 
2 /ay2 
Us 


u 

and new independent variables ~ and n by 

& = (x! — 1x7)/2, 9 = (x! + t£*)/2. 
The set of equations (5), (6), (7), (8), is then equivalent to 
(9) —VU,_ + Uv; = A/4, 
(10) t/t = 
(11) Ure/U 
(12) Ulyt + UVgt — 

The solution of (9) is 

(13) v = Ky*(m)Le4(é) {RK(n)L(E) + 1 + mK(m) + nL(é)}, 


where K and L are arbitrary analytic functions, and k, /, m, n are constants 


satisfying 


(14) mn — kl = A/A4. 
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The constant 4/4 which appears in ds? through the h of equation (4) may be 
absorbed into the arbitrary function v. One then finds, without any loss of 






generality, 
(4)’ | x0” 











(14)’ mn — kl = 1. 






Equations (10) and (11) are satisfied by a u of the same form as 2, 
(15) u = Ky *(n)Le4(8) {pK ()L(8) + q + rK(n) + LQ}. 


The K and L are the same as in (13). The constants , g, 7, s are not restricted 
by a relation like (14), but (12) is satisfied if and only if one has 


(16) ms + nr = kq + lp. 


If the line element is to be real and to have the correct signature, u? and v? 
must both be real and positive. Thus “ and v must be real, and we take 





















K(n) = L*(é), m = n*, r= s*, 






and assume R, 1, p, g to be real. 
Provided that k ¥ 0, one can write uw and vy in simpler forms. Define a new 

function H by 

(17) H = k(K + n/k). 







Then substituting in (13) and (15), 
(18) uw = 1/0 = H,A,*/(|HH* — 1), 
(19) vy = uw2/y = (P(HH* + 1) + QH + QO*H*|?/[HH* — 1). 







P is an arbitrary real and Q an arbitrary complex constant. Restrictive con- 
ditions like (14) and (16) are no longer necessary. The solution (18) is given 
by Forsyth (1906, Art. 225). 







A PARTICULAR SOLUTION 





Now interpret x! as the radial co-ordinate, x? as the z co-ordinate, and x* 
as the angular co-ordinate of a system of cylindrical polar co-ordinates. Choose 
the H in (18) and (19) to be 






N 
(20) H = > a,/(2n—ib,), 
j=l 


where a, and 6; are real constants, and take P = 0, Q = 1. The resulting u 
and v are singular only at the points x! = 0, x? = b,, and on the surfaces where 
HH* = 1 which surround each of these points. One should compare this with 
the spherically symmetric Schwarzschild solution, where, too, a point singu- 
larity is surrounded by a singular surface. 







CONCLUSION 
It is possible that the solution given in the last section describes a set of 
N particles alone in an otherwise empty universe. It satisfies the correct 
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boundary condition: at great distances from the particles both uw and » tend 
to zero, which is what one expects if Mach’s principle is valid. However, it is 
arguable that Mach’s principle would imply that displacements in the x* 
direction should make no sense (v = 0). This is a difficult point which will 
perhaps be discussed in a later paper. 

The distances between the singularities of the solution increase or decrease 
linearly with time. If we do interpret the singularities as representing particles, 
then in Newtonian language we must say that they exert no forces on one 
another. The usual gravitational attraction depends on the existence of 
matter off the line joining the attracting particles. 
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BACKSCATTERING FROM A CONDUCTING CYLINDER 
WITH A SURROUNDING SHELL! 


M. A. PLoNuS 


ABSTRACT 


Far-field backscattering from a perfectly conducting cylinder with a surround- 
ing shell has been investigated. The spacing of the shell from the cylinder and 
thickness of the shell are arbitrary. The material in the shell is also arbitrary and 
is characterized by the propagation constant h. The incident plane wave is at 
right angles to the cylinder, and is either horizontally or perpendicularly polar- 
ized. When the shell is thin in units of wavelength a much simpler expression for 
the backscattered field coefficient is obtained. It was possible to express this 
coefficient in a form which resembles the coefficient from the conducting cylinder 
alone plus a perturbation term due to the shell. Another simplification resulted 
when the propagation constant h of the shell is much larger than the free-space 
propagation constant k. 

It was desirable to see what scattering properties a cylinder with a surrounding 
shell exhibits. The cylinder was chosen to be large with respect to wavelength 
and the shell spaced a resonant distance from the cylinder. The scattering cross 
section, for this particular combination of parameters was then given by a 
slowly converging series which proved too lengthy for hand-computation, and 
was then programmed for and computed by the IBM 704. The scattering cross 
section versus shell spacing is shown in graphical form. 


FORMULATION OF THE PROBLEM 


The problem is a two-dimensional one, and can best be carried out in 
cylindrical co-ordinates. Figure 1 shows a cross section of the cylindrical 


structure. 


oP (r, 8) 


Incident 
Wave 


Fic. 1. Cross section of the cylindrical structure. 


A cylindrical shell surrounds a perfectly conducting cylinder of radius a. 
The lossy shell which is assumed to be thin is spaced a distance 6—a@ from 
‘Manuscript received August 26, 1960. h : 
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the conducting cylinder. Region I is free space with a propagation constant 
k, region II which is the shell has a propagation constant hf, and region III 
has the same propagation constant as free space. The propagation constant h 
can have a loss component and arbitrary permittivity «, however, u in region 


I and region II are identical. 


INCIDENT E VECTOR PARALLEL TO AXIS OF CYLINDER 
The solution of this problem is carried out by matching the solutions of 
the various regions at the boundaries. When the time dependence is e**, 
the direction of incidence is normal to the axis, and the £& vector of the incident 
field is parallel to the axis of the cylinder the solutions in the three regions 


are: 


ik be 2 ind 
E} es et? cos (8 +> A,H; (kr)e™ : 


n=—c 


P= > (BiH (hr) +CaJn (hr) ye 


oo 


St = SS (DH (kr) +EnJIn(kr) oe. 


“42 


n=—oD 


The incident field can now be expanded in terms of Bessel functions, thus 


D 


kr cos (8-6) . in(B—@) 
e! r cos - z i" J, (kr)e™ 


Bs 
The angle of incidence 6 can be chosen as zero without loss of generality. At 
the boundaries of the various regions the electric and magnetic fields must 
be continuous and the electric field must vanish on the surface of the per- 
fectly conducting cylinder. These boundary conditions will provide enough 
equations to solve for the unknown coefficients in the various regions. The 
two-dimensional scattering cross section is usually stated as 
| rs | 2 
o = lim 2ar| a 


4 


T3900 


Thus the scattering cross section for this problem can be written as 
Rl ee 2, | 
c= — "A, 
Tw ln=-o 
The scattered far field will be known if the coefficient .1, is determined. The 
boundary condition at r = a is 
E;"*(a,¢) = 0. 
Continuity of fields at r = 0 yields 
E:""(b, 0) = E:"(b, 6) 
0F:'"(b,0) dE," 


Hi" (b, 0) = HE"(b, 6) = ee 
or or 
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Continuity of fields at r = ¢ yields 






Ei(c, 0) = Et*(e, 0) 
a1/ 
Ho(c, 0) = He"(c, @) or OEM) a Es! — (¢, 0). 
or or 
This gives five equations in five unknowns. After considerable manipulations 
these can be solved for the coefficient A,. The above shorthand notations 
should be interpreted as 








GE(c, 0) _ dE(r, 6) 
_— = or 
Further abbreviations that will be used are 


I'(kr) = OT (kr) 



















H'(hc) = 










Ay 
i.e. all subscripts m will be left off. The resulting five simultaneous equations 


are thus: 






DH(ka)+EJ(ka) = 0 
DH(kb)+EJ(kb) = BH(hb)+CJ(hbd), 
UJ (kce) +AH(ke) = BH(hc)+Cl(he), 
DH' (kb) +EJ'(kb) = BH'(hb)+CH'(hb), 
i"J'(ke) + AH" (ke) = BH'(he)+CJ'(he). 




















When these are solved for .1 we get the following result 


(he) +& t|H (he LE: (ke) — H"( (hc) J(ke)) 














n J( (ke) J” (hc) — J' (ke). 


Awa H(kc) J’ (hc) —H' (ke) J (he \+¢|H (he)H (ke) —-H’ (hc)H (ke) ] 







where & is 






(2) §= 
J (ka)|J"(hb)H (kb) — J (hb) H"(kb) | +H (ka) | J (hb) J" (kb) — J"(ib) J (Rb) 
Tika) |H(kb)H (hb) ~H(hb)H (kb) | -+H (ka) |H (ib) J’ (kb) wer T(kb)] 







Since the lossy shell was assumed to be thin, i.e. 






hé = h(c—b) < 1 





we can expand all functions with argument ¢ about 6. Thus to a first approxi- 


mation we can write 







9] (hb) 
T(he) = I(hb)—h(b—c) oe) = (hb) +87"(hd) 


(3) y —_ Ohb 
H'(kc) = H'(kb)+65H’' (kb). 
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These are then substituted in the expressions for A and £. The complicated 
expressions which result can be simplified when the various terms are‘identified 
with the Wronskian and the derivative of the Wronskian, i.e. 


J (kx)H' (kx) — J'(kx)H (kx) = at 


as 
J (kx)H" (kx) — J" (kx)H (kx) = = 


J'(kx)H" (kx) — J" (kx) H’ (kx) = te) Qi 


x TX 


The last expression can be derived from Bessel’s equation 
x?J"" (kx) +xJ' (kx) +[(kx)?—n?|J (kx) = 0. 
The final result for A is then 
Ro me 
J (ka) +6 $= ((hj) J (ka) — (jj)H (ka)] 


(4) —¢% =b 
H(ka)+6 oF [ (hh) J (ka) — (jh)H(ka)]} 


where 


(hj) = H(kb)J''(kb)+J(kb)H' (kb) it J (kb)H (kb) (w—% ), 


(ji) = J (eb) J” (kb) + J (kb) J’ (kb) ++ F°(Rb) =) 


(hh) = H(kb)H''(kb) +H (kb)H’ (kb) *+H*(Rb) (x-#%) ; 


(jh) = J(kb)H"' (kb) +H (kb) J'(kb) 5+ J (kb) H (Rb) (v4) ‘ 


The answer which has been obtained for A can be easily checked. When the 
thickness of the lossy layer goes to zero, i.e. 6 goes to zero, or when the propa- 
gation constant h of the lossy layer approaches the free-space propagation 
constant k, A should approach the scattered field coefficient for a perfectly 
conducting cylinder of radius a in free space. Indeed when letting 6 = 0, or 
h = k the above A simplifies to 


_» J(ka) 


4a =<, 
H (ka) 
which is the scattering coefficient for a perfectly conducting cylinder above. 
A different kind of check can be performed by letting the propagation con- 
stant h of the shell approach infinity (short wavelength limit), thus making 
the shell a dielectric with an infinite dielectric constant which is a perfect 
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reflector and should behave like a perfectly conducting cylinder of radius 6. 
Again, when the limit operation is performed A simplifies to 


_ x Jeb) 
H (kb) 

If the free-space scattering coefficients A are to be computed it would be 
desirable to avoid Bessel and Hankel function terms which have to be differ- 
entiated. When the coefficients in A are re-expressed using Bessel’s equation 
and the Wronskian of argument kd a further simplification results 


A= 


(Aj) — 5+ J (kb)H (kb) (i —#), 


(jj) = J°(kb) (h°—k’), 
(hh) = H*(kb)(h?—k’), 


bi) t+ J (Rb)H (Rb) (i? 8’). 


Dividing the numerator of A by its denominator and neglecting terms of 
order higher than 6 one obtains 


i ites ek Pe a ‘) 
A= (ka) 1 O35; T(ka)H (kay (7 (#0)# (ka) J (ka)H (kb) | 


This is a simple enough result which would lend itself easily for computation. 
The material of the thin shell was assumed to be non-magnetic, of arbitrary 
dielectric constant e, and conductivity a, i.e. 


Me = Kh 
k wv (uo€o) 
h = wv (ue*) 


where 


e* = @e,—(to/w). 


An interesting observation can be made now. In the above expression for 
A,, when 6 = a, the thin shell is in contact with the surface of the perfectly 
conducting cylinder. The terms that multiply 6 vanish now. Thus, when the 
E vector of the incident field is parallel to the axis of the cylinder one can 
say that to a first-order correction the scattered field of a perfectly conducting 
cylinder coated with a thin dielectric or lossy shell is the same as the scattered 
field of a perfectly conducting cylinder alone. 

INCIDENT E VECTOR PERPENDICULAR TO AXIS OF CYLINDER 
The mathematical manipulations for this polarization are even more 
involved than.in the previous case. The initial equations in the various 


regions are now 
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Hi= > i), (krle™+ D>. AH (kre, 


n=—c n=—co 


co 


Hy" = Q0 (Buln (hr) +C,Ja(hr) )e™, 


r=—a 


oc 


HP? = ¥ (DHS (kr) +E,Ja(kr))e™. 


n=—co 


The boundary condition on the surface of the perfectly conducting cylinder is 
El(ka, 0) = 0. 
Matching the fields at r = 0 
Hi""(kb, 0) = H}" (hb, @) 


III II 
E}'" (kb, @) = Ej'(hb, 6) or 2 oH, (kb, 0) = 2 oH, (hb, 0) 
4 or h or 
Matching at r =c 
H}(ke, 0) = H}" (he, @) 


Whe 9) = EE 1 aH: (ke, 0) _ 1 aH" (he, 4) 
Ea(kc, 0) = Eo (he, @) or B rae a. ay 


This again gives five equations in five unknowns, which are 
DH'(ka)+EJ'(ka) = 0, 
DH(kb)+EJ(kb) = BH(hb)+CJ/(hb), 
"J(ke) +A (ke) = BH(hc)+CJ(he), 


7: (DH'(kb) +EJ'(kb)) = * (BH'(hb)+CJ'(hb)), 


re (i"J'(kc) +AH'(kc)) = ‘ (BH' (he) +CJ'(he)). 


Solving these for 41 we get 


ia a ia RS (ke) J" (he )— h*J'(ke) I (he) +&(h"H (he) J'(ke) — kJ (ke)H'(he)} 


k°H (kc) J’ (he) —h°H (ke) J (he) +&[h7H (he )H' (ke) —k°H (ke) H (he) 
where & is 
(7) ¢-= 


J'(ka) {kJ (hb) H (kb) —h* J (hb) H (kb) |4+-H (ka) (hh J (hb) J'(kb) 
—k°J' (hb) J(kb)| 


J'(ka)(k°H (kb)H (hb) —h"H (hb) H'(kb)|+H' (ka) [AH (hb) J'(kb) 
k H' (hb) J (kb) | 
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All functions of argument c will again be expanded in a Taylor series about 
b. Since 






= k(c—b) <1 






only the first two terms in the expansion are kept and substituted in the 
above expression for A. Simplifying the resulting expression, we get the 
following result for the scattered magnetic field coefficient 


n J” (ka) +6[ (hj) J’ (ea) — (99) 1" (ka)] 




















(8) - OT (ka) +6 (hh) J'(ka) — (jh)H (ka)} 

where 
(hj) = +n (k°—h? 1G J'(kb)H' (kb) +e S | (kb) 1(88)) é' 
(77) = aE — (k°— wy (W (Rb) + J° *as)) 
(hh) = ay 2 (Reh? (42 "HT! (kb)? ()) 






(jh) = stop sb (eH) (way Ca) +8 IBY) 





Dividing the numerator of A by its denominator and neglecting terms of 


order higher than 6 one obtains 


Et .~2___ 2 -¥ 
’H (ka) 






( = 
(9) A By (kh) ‘te (ka)H’ (ka) 





[ie (H' (ka) J'(kb) — J’(ka)H'(kb))? “tS, (H'(ka) J (kb) ~7'(kayH8))*}). 





This expression is more complicated than the corresponding one for the other 
polarization, but still adaptable for hand computation, at least for properly 






chosen ka and kd. 

When the corresponding expressions for the two different polarizations are 
compared, one finds that for the perpendicular polarization, 4 has an additional 
term which does not vanish when the shell is in contact with the surface of 
the cylinder. Letting b = a, we get the scattered field coefficient for a coated 






cylinder 


ee F(a) ( Zin" ( h “—k " as 
(10) A= i H (ka) I+6— nah? J’ (ka)H (ka) 


Thus for perpendicular polarization one should expect different scattered fields 
for a coated cylinder and a bare cylinder, whereas for the case of parallel 
polarization the scattered fields are the same to a first-order approximation. 
However, if second-order approximations would have been considered even 
the parallel polarization case would have yielded a second-order correction 
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term. This should be expected since a coated cylinder and a bare cylinder are 
different physically and this difference must show up in the mathematics 
eventually, even if only as high-order correction terms. 

Another expression for A which might be useful is obtained when eggs. (6) 
and (7) are so rearranged that terms which contain the propagation constant 
h and terms which contain the propagation constant k are grouped together. 
Such an arrangement might prove to be useful when one of the propagation 
constants is much larger than the other one and asymptotic expressions are 
to be used. Thus, eqs. (6) and (7) for perpendicular polarization can be written 

in J (ke)a+ J" (ke) B+p(I(ke)y+J"(ke)r 


(11) A= —" Fi (ke)a-+H (ke)B-+p(H (ke)y-+H" (en 


where 


_ b'(J"(ka)H" (kb) — J (kb)H" (ka)) 
~ (J (kb)H (ka) —H(kb)H (ka)) ’ 


k°(J'(hb)H'(he) — J'(hc)H'(hb)), 
h*(J(hc)H' (hb) —H (he) J'(hb)), 
k* (J (hb)H' (he) — J'(he)H (hb)), 
h*(J(he)H (hb) — J(hb)H (he)). 





When 
h>k 


ka > 1 


we can use asymptotic expressions for the Bessel and Hankel functions of 
argument hb and he. If we substitute 


J,(hb) = i 2 = Z coal ib—n =" =) 


(2) as 2 
Fe the): ace 


t(—he+n5+5) 


in the above expression we obtain 


k’h® J (ke) sin hi+ h° J'(ke) cos hé + p[k* hI (ke ) cos hi—h'J' (ke) ) sin hd) 


k’h’H (ke) sin hb-+h°H (ke) cos hi + pik’ “hH (ke) cos hs —h°H’ (kc) sin hé]} 


where 


hé = he—hb. 
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AN APPLICATION 
It was desirable to see what scattering properties a cylinder with a surround- 
ing shell exhibits. To make a sample calculation for a particular cylinder the 
defining parameters are chosen as 


reflection coefficient |R|? = .8, 
dielectric constant for shell (€/¢o9) = 100, 
ka = 12. 


The shell was spaced at distances 4/2, A, 3/2, but to catch any resonance 
effects a few [increments] were taken on each side of the spacing distances. 
Stratton gives formulas which relate the reflection coefficients R to the 
thickness 6 -i the shell. He considers a thin plane dielectric sheet separating 
two semi-infinite dielectric media. 
The reflection coefficient R for the shell is then given as 
(nietres)*— 4rjore3 sin "Rod 


R= eo 
(1+ riers)” —4riere3 sin”k2d 


where the small lettered r’s are the reflection coefficients at each boundary. 
When the specified values of R and (€/e) are substituted in the above ex- 
pression one can solve for the thickness 6 of the shell in units of wavelength, 
thus 


hé 416 


and since 
= 10k 


then also 
ké = .0416. 


Since for the particular value of ka chosen, hand-computation proved too 
lengthy, the problem was programmed for the IBM 704. The computer cal- 
culated eq. (5) as two sums, i.e. 


co 


« n as 10 - bel n J, (ka) 
(13) do "An = 2 | 1) H, (ka) TG - I. 


nes 


The advantage of this calculation is that the thickness and material of the 
shell could be specified after the calculation by the constant 6;, where 4; is 


given by eq. (16). 
For perpendicular polarization the computer calculated eq. (9) as three 


sums, 1.e. 


y Ji (ka) 
"“A,=- —1)"=4 os) +be(...) 
(14) 2, tA, ant 1) i'(ka) b+ +82(...)] 


rn=—c n= 


where 52 is given by eq. (19). 
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For parallel polarization eq. (5) the result was presented in the form 
(15) a+ib+6,(c+id) 


where 












997 


21 





(16) 5, = 882 gt_ yt) = BE asus). 


;: 


We are interested in the scattering cross section, which now becomes 





2r : a 
eee la+ib+6,(c+-7d)|*. 





Carrying out this operation and dropping terms of order 6; we get 





_ = (a°-+b°+26;(ac+bd)). 





This result can then be interpreted as the cross section of a perfectly con- 
ducting cylinder alone and the contribution to the total cross section from 






the shell, thus 


: 2r 2 2 2r 
(17) o = Teyt +O shen = ‘ee ed r= 26; (ac+bd) 





where oa., is well known, i.e. 





eyl 





ore eee 
ee € aes T dX (-1) H\” (ka) | ° 


n=—co 








For perpendicular polarization, eq. (9), the results are given as 


(18) at+ib+6,(c+71d)+62(e+7/) 


where 














997 


(h°—k’) = — (kd) (kb) 


‘ae 
aig ae 


21 








and 





’ og RHR _ 99 rhb 
(19) s=6o i “ina 











The cross section is then again given as 









9 x ! 

(20) ea = (a?-+-b?+25, (ad—bc) +252 (af —be)) 
here 

BN aegis: nema 

Toy = - (a +6 ) =< r z ( 1) H, (ka) . 






No significance could be attached to the cross section obtained by this 






method when 





hé = .416. 
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Equations (5) and (9) were obtained by expanding the exact expressions 
in powers of 4é and retaining terms of order hé only. For hé = .416 this 
would clearly not be valid. However, after the computer calculated eqs. (13) 
and (14) and presented the results in the form of eqs. (15) and (18) a signi- 
ficant conclusion could be drawn: the constants c and d in eq. (15) and the 
constants c, d, e, and f in eq. (18) were comparable in magnitude to the con- 
stants a and 3b, in most cases smaller by three or four magnitudes. Thus, when 
the shell is thin (i.e. h6 <1) the backscattering cross section for either 
polarization is approximately equal to that from the perfectly conducting 
cylinder alone, the shell having only a small perturbation effect. 

To calculate the cross section when 


hé = .416 
and ((€/€9) = 100) the IBM 704 was programmed to calculate eq. (12). Equa- 
tion (12) demands only that 

i> Rs 
In this calculation 2 = 10k. The computer calculated a set of values for 


shells spaced from .46\ to .64\ from the perfectly conducting cylinder. The 
results are shown in Fig. 2 in graphical form. This can be compared to a 


ieee 


For Solid Cylinder 
of ka = 15.14 


For Solid Cylinder 
ka = 12 


(b—a)/X 
Fic. 2. Variation of the radar cross section as a function of the shell spacing from the 
cylinder. 
perfectly conducting cylinder of ka = 12 with a perpendicularly polarized 
incident wave which is 
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and to a perfectly conducting cylinder which is as large as the shell, i.e. 
kb = 15.2 which gives 


Ji (15.2) |? 
in 1)" G52) 2) = 11.9. 


n=—oco 


This particular set of parameters shows some resonance phenomena which 
Fig. 2 displays. The average gain of the shell plus cylinder combination over 
a cylinder of the size of the shell alone is approximately 


10 log i ll 


whereas the peak gain is 


or 
10 log aa = 2.02 db. 


The values computed by the IBM 704 are tabulated in Table I. 


TABLE | 





Shell spacing Shell spacing 
(b—a)/d (x/2d)o | (6— a)/X 








.46 3.923 08 
47 31 . 584 
.48 21 .588 
49 .98 .59 
.50 .O1 .592 
51 9.65 595 
.52 .02 .60 
.53 .057 .61 
54 .405 615 
55 11.52 .618 
.56 12.885 .62 
.564 14.72 . 622 
.568 16.73 .625 
.O7 16.88 63 
.572 16.25 . 64 
575 14.385 
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THE CANADIAN STANDARD FREE-AIR CHAMBER FOR 
MEDIUM QUALITY X-RAYS! 


W. H. HENRY AND C. GARRETT 


ABSTRACT 


A description is given of the Canadian standard free-air chamber for measure- 
ment of medium quality X-rays, including measurements of the contribution 
from radiation scattered from the air and from the diaphragm, and a method for 
the accurate alignment of the chamber and X-ray source. 


INTRODUCTION 

The free-air ionization chamber to be described is the Canadian standard 
instrument for measuring X-radiation of medium quality, and forms part of 
the X-ray standardization equipment at the National Research Council 
Laboratories. 

The use of X-rays in such fields as medicine and biology requires measure- 
ments of exposure dose. Exposure dose is essentially a measure of the ionizing 
ability of X-radiation in air. The practical measurement of exposure dose is 
complicated by the fact that one is required to collect all the ion pairs pro- 
duced by the electrons which are themselves released in a known mass of air 
through the interaction of X-radiation in that mass. In the range of X-ray 
qualities below 500 kv the best instrument for the absolute measurement of 
exposure dose is the free-air chamber. This is essentially an ionization chamber 
with air walls of a thickness roughly equal to the range of the electrons released 
by the X-radiation. Owing to its large size it is usually used only in the labora- 
tory, and as the standard against which more portable instruments can be 
calibrated. 

The accuracy of a free-air chamber depends upon the fulfillment of a number 
of requirements. Data relevant to the design of a chamber have been published 
by Attix and DeLaVergne (1954a, 6) and Wyckoff and Attix (1957). Briefly, 
the chamber must permit an accurate knowledge of the mass of air from which 
the electrons are released by the X-radiation, ensure the collection of ions 
equal in number to the total number produced by these electrons, and have 
a means of accurately measuring the resulting ionization current. The volume 
of the mass of air is defined by the effective area of the entrant aperture and 
by the limiting lines of the electric field ending on the edges of the collecting 
plate. The electric field in the collecting region must therefore be uniform. Its 
uniformity depends upon the smoothness and alignment of the collector, guard 
and high potential plates, and in a parallel plate chamber it is increased by the 
use of some system of guard leops surrounding the space between the plates. 
The electric field must also be of sufficient strength to reduce ion recombination 
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to a low value. The dimensions of the chamber at right angles to the beam 
must be sufficiently large for most of the energy of the released electrons to 
be expended within the air before striking the plates or the guard loops. The 
ionization produced by radiation scattered out of the main beam should be 
low, and radiation other than that passing through the entrant aperture must 
be reduced to a negligible value by adequate shielding around the chamber. 


GENERAL ARRANGEMENT 


The X-ray standardization equipment is housed in an air-conditioned room 
where the temperature variations have been reduced to about +0.5° C. Figure 
1 is a general view of the setup. One side wall of the standard air chamber 


Fic. 1. General view of X-ray standardization setup, showing the standard chamber near 
one end of the lathe bed and the two X-ray tubeheads mounted on rails perpendicular to 
the lathe bed. 


has been removed, exposing the rear surface of the guard plate. The mounting 
of the chamber and its supporting carriage on a lathe bed parallel to the X-ray 
beam permits accurately reproducible positioning of the chamber at distances 
of 50 to 300 cm from the focal spot. Adjustments are available for varying 
the vertical position and angular alignment of the chamber. The standard 
chamber can be moved off the carriage, out of the beam, and onto a table at 
the end of the lathe bed. Secondary chambers to be calibrated can then be 
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mounted on the same carriage and placed at the position occupied previously 
by the standard chamber diaphragm. A removable micrometer pointer assists 
in accurate positioning. 

The box enclosing the plate system is lined with ,$-in. lead sheet to 
attenuate stray radiation. On the front wall a thickness of 4 in. of lead has 
proved adequate under all normal conditions to reduce ionization due to 
primary beam penetration to a few hundredths of one per cent of the total 
ion current. The inside dimensions of the box give a clearance of at least 
10 cm between the open sides of the plate system and the box. 

Two interchangeable sources of X-rays are available: a 250-kv constant 
potential generator, and a 100-kv unit with a beryllium-windowed tube for 
which half-wave or full-wave rectification or constant potential operation can 
be chosen. The nominal focal spot sizes of the tubes are respectively 3.5 mm 
square and 3.5mm by 1.5mm. The two tubeheads are mounted on rails at 
right angles to the lathe bed so that either source can be placed in position. In 
order to facilitate alignment, the tubeheads are adjustable within the support- 
ing framework. 

In an X-ray standardization setup, either the exposure dose rate of the 
X-ray beam must be very well stabilized, or the fluctuations in it must be 
measured and suitable corrections made. In the present equipment each X-ray 
generator incorporates a tube current stabilizer and a line voltage stabilizer 
which reduce fluctuations and drift in the X-ray exposure dose rate to few 
tenths of one per cent. Greater stabilization than this would be costly and is 
made unnecessary by the use of a monitoring ionization chamber, permanently 
mounted in the path of the X-ray beam at the end of the lathe bed near the 
tubehead. The chamber consists of three parallel plates perpendicular to the 
beam: two outer, 1 mm thick, beryllium plates at a potential of 300 v, separated 
about 1cm from the central collector of aluminized mylar film. This adds 
very little extra filtration to the X-ray beam. The fluctuations in the ratio 
of the standard chamber ionization current to monitoring chamber ionization 
current are of the order of 0.1%. 


PLATE SYSTEM 

Figure 2 is a horizontal sectional schematic view of the chamber, showing 
the main dimensions of the parallel plate system. The 26-cm separation of the 
plates, which are 35cm in height, maintains the loss of ionization due to 
inadequate size of the chamber at generally less than 0.25%, depending upon 
the quality of the X-radiation. Each plate is made from a ribbed aluminum 
casting which has been aged by heat treating and machined smooth and 
plane to within 0.001 in. The central collector plate is supported in the well 
of the guard plate and co-planar with it by polystyrene insulators. The gaps 
between the edges of the collector and guard plates are approximately 0.25 mm 
wide and were made as narrow as possible to reduce uncertainties in the 
collector length, consistent with the retention of high insulation. Since the 
plates have a thickness through the strengthening ribs at the edges of about 
4cm, the gaps have been tapered a few degrees toward the outer surface of 
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Fic. 2. Schematic horizontal sectional view of standard chamber. 


the plates. This enabled conducting dust particles present after assembly to 
be removed easily. 

Since the effective collector length, being determined by the limiting lines 
of force ending on the collector, will be altered by any distortion of the electric 
field, the collector and guard plates have been made co-planar by aligning 
the adjacent edges to within 0.0005 in. along their entire lengths. In the 


central region where most of the ions are collected the alignment along the 
edges is within 0.0001 in. The high potential plate is parallel to the guard 
and collector plates to within 0.001 in. The high potential and guard plates 
are held in this fixed relative position by four micalex insulators located at 
the corners. In this way the three plates form an independent assembly, shown 
in Fig. 3, whose mechanical stability has proved to be excellent. 

The plate system as first constructed had two separate guard plates. This 
required four additional micalex supporting insulators, between the high 
potential plate and the corners of the guard plates adjacent to the collector. 
The movement of charges within these insulators induced currents in the 
collector which were not greatly reduced by guard wires, but which were very 
nearly eliminated by the guard bar system described below. In order to 
remove not only the remainder of this effect, but also the possibility of field 
distortion in the collecting region due to charge distribution on the insulators, 
the present plate system was constructed with the single guard plate casting, 
thus eliminating the four insulators nearest the collecting region. A comparison 
of the ionization currents measured with the two plate systems showed no 
actual difference greater than 0.1%, the limit of accuracy of the measurements. 

The chamber has normally been operated with an electric field of 192 v 
cm~!, supplied by a highly stable, 5000-v source. An estimate of the loss of 
ions due to lack of voltage saturation has been made by plotting the reciprocal 
of the measured ionization current versus the reciprocal of the electric field 
strength, for fields up to 346 vcm~. The curve becomes nearly linear for 
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Fic. 3. Photograph of the parallel plate system. 


fields above 100 v cm~', and an extrapolation to infinite field strength indicates 
a loss in collected ions at 192 vcm™ of 0.2% for an exposure dose rate of 


20 r min—. 


GUARDING 

In order to make the electric field uniform and perpendicular to the collector, 
a system of guard bars surrounds the volume within the plate system. A 
photograph of the guard bars in place within the plate system is shown in 
Fig. 4. Wires (Taylor and Singer 1930) and strips (Wyckoff and Kirn 1957) 
are commonly employed, and Failla and Marinelli (1937) used thin rods. Bars 
offer the twofold advantage of high electrostatic shielding from the effects 
of external disturbances, while requiring a minimum use of insulating spacers. 
Twenty-four parallel rectangular loops, made from 3-in. square aluminum 
stock, are equally spaced at approximately 1.1 cm, center line to center 
line. This leaves a small gap of approximately 0.13cm between adjacent 
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Fic. 4. Guard bar system mounted within the plate system. 


bars. The inherent stiffness of the bars allows the use of supports and insulating 
spacers only at the four corners of the system. This forms an independent 
assembly, shown in Fig. 5, which is easily removed as a unit from the plate 
system, making possible the insertion of other trial systems of guarding. 
Holes have been made in the guard bars at the front and rear to allow 
passage of the X-ray beam, while electrical continuity is maintained by offset 
wires shunting the holes. The possibility of the holes distorting the electric 
field in the collecting region was investigated by making comparison measure- 
ments of charge collected with fine aluminum wires stretched across the larger 
exit gap. No effect could be detected. A potential divider sets the potentials 
of the guard bar loops to give a linear gradient between the plates. Since the 
spacing of the center line of each outside loop from the surface of the plate 
is one-half the spacing between center lines of adjacent loops, the potential 
difference between each outside loop and plate is made one-half the potential 
difference between loops. The wire-wound resistors of the potential divider 
chain are mounted outside the chamber to avoid heating of the chamber air. 
The electric field close to the guard bars is, of course, distorted by the 
presence of the bars. If the distortion extends into the collecting volume, then 
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Fic. 5. Guard bar system. 


the length from which ions are collected will be unknown. Measurements of 
ion current were made with the guard bars connected in a way that would 
considerably increase this distortion. In the first case, loops were shorted 
together in pairs, thus effectively forming half as many loops of approximately 
twice the width. In the second case, groups of four were shorted together, 
effectively forming only six guard loops, each of about four times the width. 
Table I shows the values of the measured ion current, normalized to a fixed 


TABLE I 


Jon current with different guard loop arrangements 


Measured ion current, 
Guard loop arrangement arbitrary units 





Normal, connected singly .851+0.002 
Shorted in pairs .852+0.002 
Shorted in groups of four .855+0.002 


X-ray exposure dose rate by means of a monitoring chamber. It is evident 
that the error due to field distortion from the bars themselves is probably 


considerably less than 0.1%. 
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The proximity of the grounded outer box also can cause field distortion. 
The resulting error in collecting volume is usually measured by comparing 
the ion currents collected with the box at ground potential and at the potential 
of the high voltage supply. One-half the difference is taken as the error. Since 
it was not convenient in this case to operate the box at high potential, the 
same purpose was accomplished by temporarily mounting aluminum sheets 
on insulators inside the box about 1 in. from the wall. The sheets could be 
set at either ground or high potential, and the indicated error from this type 
of field distortion was less than 0.05%. 


DIAPHRAGMS 


The volume of the mass of air in which X-ray interactions take place is 
determined not only by the length of the collector but also by the area of 
the entrant diaphragm. Cylindrical apertures of 7 mm and 10 mm diameter 
have been machined in an alloy of tungsten (General Electric Hevimet), that 
are round and smooth within 0.5 uw throughout their lengths. The diameters 
have been measured within 1.5 4 by the Electricity and Mechanics Section 
of the National Research Council. 

The cylindrical bore of each diaphragm was initially 10 mm long, followed 
by a 3-mm section at the rear (i.e. away from the X-ray source) tapered by 
45°. In order to remove slight unevenness in the bore of the 7-mm diaphragm, 
No. 7A, shown in Fig. 6(a), the front section was tapered by about 3.5°, 
leaving only a short cylindrical bore of satisfactory roundness and smoothness 
about 1 mm in length. An added advantage of a short cylindrical bore is the 
greatly eased problem of alignment (Aitken 1958), since the error introduced 
by a misalignment is approximately proportional to the length of bore. 

Scattered and fluorescent photons from this diaphragm contributed signi- 
ficantly to the total ion current. Such a possibility has been pointed out by 


(a) Shape of diaphragm No, 7A. (b) Shape of a low scattering diaphragm. 
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Greening (1960), and the measurement of this contribution is described in the 
next section. To reduce this effect to a negligibly low value, the diaphragm 
design has since been improved following Kénig (1939), as shown in Fig. 6(8). 
The cylindrical bore of only several millimeters length is preceded on the side 
toward the X-ray source by a short section tapered at 45°, and followed by a 
section of about 8-mm length tapered at 5°. Five degrees is sufficient to avoid 
incidence of the X-ray beam on the rear conical surface at the smallest used 
source to chamber distance. The surface of the 45° tapered section can not 
‘‘see’’ into the chamber and so secondary photons from this surface are 
attenuated by the walls of the diaphragm. This leaves only the short cylindrical 
section as a possible source of scattered radiation. The rear conical section and 
the cylindrical bore are concentric within 3 » and blend with a radius no greater 
than 0.7 mm, thus giving a well-defined plane of definition. 

Although diaphragm 7A itself, when reversed, is similar to the low scattering 
diaphragm, it is not used in that position owing to the intersection of the 
cylindrical bore and the added 3.5° taper not giving a well-known defining 
plane. 

A diaphragm on the tube housing limits the emerging X-ray beam to a 
convenient size. The larger the beam cross section, the greater is the contri- 
bution from penetration of the front wall. On the other hand, an aperture 
smaller than the focal spot would shade part, or all, of the entrant diaphragm 
from some points on the focal spot. The resulting change in effective area of 
the entrant diaphragm would not be in proportion to the change in exposure 
dose rate at its center. 

























MEASUREMENT OF SCATTERED AND FLUORESCENT PHOTONS 










The contribution to the total ion current from scattered and fluorescent 
photons was measured by a method described by Attix and DeLaVergne 
(1954a, 6). A plastic tube was mounted inside the chamber, between the 
entrant and exit diaphragms, and of such diameter to just clear the X-ray 
beam. The tube was coated with graphite and operated at one-half the potential 
of the high tension plate in order to minimize distortion of the electric field. 
The wall was sufficiently thick to stop all secondary electrons ejected from 
the air within the tube. Using entrant diaphragm 7A the observed ionization, 
due only to scattered and fluorescent photons and corrected for the shadowing 
effect of the tube and for the estimated attenuation of the photons in the 
walls of the tube, is shown in the third column of Table II. 














TABLE II 


Scattered and fluorescent radiation 












Contribution from scattered and fluorescent photons, 
% of total ion current 






Measured 







Constant potential BV .ba contribution from 
(kv) (mm) Measured Calculated diaphragm 7A 
80 1.6 Al 0.72 0.59 0.03 
140 0.7 Cu 0.83 0.46 0.18 
0.75 0.40 0.19 





250 2.0 Cu 
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The calculated values in the fourth column of Table II are based on the data 
of Wyckoff and Attix (1957) for the contribution from secondary photons 
arising from the interaction of the X-ray beam and the air within the chamber. 
Their X-ray qualities differed slightly from those in Table II, but since the 
contribution varies slowly with quality, the necessary interpolation introduces 
only a small error. Part of the difference between the measured and calculated 
values is due to scattered and fluorescent photons from the entrant diaphragm, 
7A, and also to the length of the chamber differing from that used by Wyckoff 
and Attix. 

The scattered radiation was further investigated by studying its variation 
along the length of the chamber. Lead sleeves, of 10- to 12-cm length and 
1.6-mm wall thickness, were fitted over the entire length of the plastic tube 
except for a 2-cm gap at the front guard bars where there was not sufficient 
available space. The removal of a sleeve allowed unattenuated passage of the 
scattered photons through that particular section. Figure 7 shows the results 
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Fic. 7. Contribution to total ion current from scattered and fluorescent radiation passing 
through sections of the plastic tube along the chamber length; at 140 kv, 0.7 mm Cu H.V.L., 
with entrant diaphragm 7A. The gap between 12 and 14cm corresponds to the position of 
the:front guard bars, and the position of the collector is from 34 to 44 cm. 


obtained at 140 kv. As expected, the contribution of scattered photons is 
largest from the central region, opposite the collector, and falls off toward 
each end of the chamber. The large increase at the front of the chamber is 
due to the scattered and fluorescent radiation from entrant diaphragm 7A. 
An estimate of the scattered radiation (including the fluorescent radiation) 
from the diaphragm can be made from Fig. 7 by subtracting the extrapolated 
contribution from the air. However, more accuracy is obtained by comparing 
the total scattered radiation, from the air and diaphragm, measured with 
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diaphragm 7A, with the total when measured using a low scattering diaphragm 
of the type described above. The difference is then equal to the difference in 
scattered radiation from the diaphragms, and so approximately equal to that 
from 7A. Column 5 of Table II shows the values so obtained, representing the 
contribution from scattered radiation from diaphragm 7A transmitted through 
the front of the guard bar system. This transmission has been estimated to be 
a few per cent at 80 kv and 60 to 70% at the higher kilovoltages. 

The scattered radiation from a cylindrical diaphragm would be approxi- 
mately proportional to the length of the bore. A low scattering diaphragm 
was compared with a diaphragm having a cylindrical bore of the same diameter 
but of about six times the length, loaned by the National Bureau of Standards. 
The difference was below the limit of measurement at all X-ray qualities, 
indicating that the contribution from scattered and fluorescent photons from 
either diaphragm is no greater than a few hundredths of one per cent of the 
total ion current. 
CURRENT MEASUREMENT 
The ionization current is measured by a standard Townsend balance type 
of circuit in which the charge is collected on a capacitor, the resulting voltage 
balanced by a compensating potential, and the null indicated by a vibrating 
reed electrometer. The capacitor is a specially designed, air dielectric, quartz- 
insulated, three terminal type of about 2000-yuf capacitance. Precise measure- 
ments by the Electricity and Mechanics Section ef the National Research 
Council have shown that the capacitance has not changed by more than 
0.01% over many months. The compensating potential is supplied by dry 
cells and measured with a differential voltmeter calibrated against a built-in 
standard cell. The feedback of the electrometer maintains the input at very 
nearly ground potential, thus reducing to a negligible value loss of charge by 
leakage from the input circuit to ground through the capacitor insulators and 
the supporting insulators of the collector. The input being held near ground 
potential by the feedback also reduces distortion in the electric field due to 
the collector and guard plates not being at the same potential, and eliminates 
the necessity of continuously maintaining a null. The capacitor and electro- 
meter preamplifier are located on top of the standard chamber inside a }-in. 
lead shield which was found necessary to reduce ionization from stray radiation. 

The ionization current from the monitor chamber is measured in the same 
way. The much larger current, however, necessitates a much larger capaci- 
tance. A commercial solid dielectric capacitor of about 1-uf capacitance is 
used and is entirely satisfactory since the resistance is adequately high, 
dielectric absorption effects are not noticeable, and long term stability is not 
a consideration. 
ALIGNMENT 


The alignment of the standard air chamber requires the placement of the 
center of the focal spot of the X-ray tube on the axis of the entrant diaphragm, 
and the maintenance of this condition for all positions of the chamber along 


the lathe bed. 
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The latter requirement is conveniently satisfied first. As the chamber moves 
along the lathe bed, the locus of the center of the defining plane of the entrant 
diaphragm is, of course, a straight line, and the focal spot must be positioned 
on this locus. Using a short entrant diaphragm of 1-mm diameter, a film is 
placed over the exit diaphragm and a double exposure made, one with the 
chamber at each limit of travel on the lathe bed. The exposure times are 
adjusted to give a small dark image within a larger, less dense image. From 
the relative positions of the images and the known geometry, the displacement 
can easily be calculated by which the focal spot is off the locus. When correctly 
positioned, the double exposure gives two concentric images. 

The axis of the entrant diaphragm must then be made to pass as nearly as 
possible through the center of the focal spot. A misalignment of a few minutes 
of arc can produce an error of 0.1% (Wyckoff and Attix 1957). The angular 
alignment is most easily accomplished by installing a long, narrow entrant 
diaphragm in order to exaggerate the error for any particular angular mis- 
alignment. The diaphragm used has a length of 17cm and a diameter of 
1 mm, made by accurately aligning small holes at the ends of a tube of larger 
diameter. The radiation transmitted by this diaphragm is conveniently 
measured with a survey meter mounted at the exit diaphragm. From a plot 
of the meter reading against angular rotation of the entrant diaphragm, the 
correct alignment,is readily determined within 1 to 2 minutes of arc. Since it 
is not possible in this chamber to rotate the entrant diaphragm independently, 
the same purpose is accomplished by moving the focal spot. This gives an 
effective .ngular rotation of the diaphragm, and in Fig. 8 is shown the measured 
variation in transmitted radiation. Then, having determined the angle by 
which the diaphragm is out of alignment, the entire chamber can be rotated 
by the indicated amount. 

The X-ray beam within the chamber is parallel to the plane of the collector 
and perpendicular to the plane of the electric lines of force ending on the 
collector edge, within 10-20 minutes of arc. Hence, cosine errors from this 


misalignment are less than 0.01%. 
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Fic. 8. Transmitted radiation versus effective angular rotation about the aligned position 
of the 1-mm diameter aperture, of 17-cm length. 
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ACCURACY 


The maximum error in the determination of exposure dose by means of a 
standard free-air chamber has been estimated in the Report of the Inter- 
national Commission on Radiological Units and Measurements (1956) to be 
likely no greater than 1.1%. This is the sum of the estimated maximum errors 
in the individual quantities involved. 

In an intercomparison between two standard chambers, since some of the 
quantities are common to both chambers, it is estimated in the same report 
that the agreement should be within 0.9%. When the Canadian standard 
chamber and the standard chamber of the National Bureau of Standards 
were intercompared, the agreement was well within the estimated 0.9%. A 
forthcoming paper will describe this intercomparison. 
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THE ‘‘REACTOR”’”’ NEUTRON CAPTURE CROSS SECTION OF Y"! 


GWEN M. MILTON AND W. E. GRUMMITT 


ABSTRACT 


A high specific activity source of Y, freed from natural uranium and rare 
earths, was irradiated in a high flux position in the NRX reactor. After further 
chemical purification, 3.6-hour Y® was identified in the sample by measurement 
of the 3.6-Mev £-particle associated with its decay, using a scintillation spectro- 
meter. Two such measurements have established the reactor neutron capture 
cross section of Y® to be 1.4+0.3 barns. 


INTRODUCTION 


As a result of its relatively high fission yield (5.8%) 58-day Y* is a potential 
reactor poison of some importance, and hence the measurement of its thermal 
neutron cross section is of more than routine interest. A study of the nuclear 
reactions involved suggests a possible method for this measurement: 


B 1.54 Mev 
y"*——--—_-—--——+ Zr" (stable) 


58 d 


.6 Mev 


8 { 2.7 Mev 
.3 Mev 
y2 ———— 7r® (stable) 
3.6 hr 


\ 


Assuming that a reasonable fraction of the disintegrations of Y” proceeds by 
emission of the 3.6-Mev 8, a small quantity of this nuclide should be easily 
identifiable in the presence of much larger amounts of Y* by using a scin- 
tillation spectrometer to select those 6-energies above 3.1 Mev (thus eliminating 
chance coincidences of two Y* 8-particles being counted as a single 3.1-Mev 
particle). 

Preliminary experiments indicated, however, that this cross section is 
<10 barns (Eastwood et al. 1958) and that with the available counting 
equipment it could not be measured in the presence of appreciable amounts 
of uranium and natura! rare earth contamination. Thermal neutron fission in 
uranium produces three yttrium isotopes which could interfere with this 
experiment, Y%, Y%, and the 3.6-hr Y® under study. (See top of p. 1691.) 
Also, Y% and Dy'® are formed by neutron capture in stable rare earths 
present in the irradiated sample of Y*. These nuclides can interfere with 
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the determination if the energy discrimination is poor, or if the coincidence 
rate is sufficiently high to produce “‘pile-up’’. This latter effect is caused by 
the numerical addition of the energy of two 8-particles within the amplifier 
and is a function of the disintegration rate of Y®, of the contaminants, and 
the resolving time of the instrument (1 usec). 

In the initial experiments, using carrier-free Y" (4X10-" g = 10° d.p.m.) 
prepared from fission product Sr®, a limit of <7 barns was set for the Y* 
cross section. As Y® and Y® from fission, and Dy'® from capture were limiting 
the measurement, a considerable effort was required to reduce the quantity 
of uranium and rare earths present in the sample prior to the irradiation. 












EXPERIMENTAL 









Chemical 

(a) Preparation of a High Specific Activity Source of Y™ 

Carrier-free Y*! can be prepared in a radiochemically pure state by milking 
it from Sr“. Unfortunately, Sr®! obtained by separation from irradiated 
uranium necessarily contains small amounts of other strontium activities and 
hence two other yttrium isotopes are collected in this milking process, Y** 
(stable) and Y*®°. However, the amounts are very small if the irradiation time 
is short (<1 day), and the time allowed for growth of Y* from the separated 
strontium is a minimum. Y°*? contamination can be reduced further by decay, 











as its half-life is short relative to that of Y™. 

Approximately 25 millicuries (5.6X10!d.p.m.) of carrier-free Sr® were 
separated from 50 mg of natural uranium and mixed fission products imme- 
diately after irradiation in the NRX reactor. The first step in the separation 
was the removal of anion complexes by their adsorption on Dowex 1 from 
10 M HCI, or on a second column, from 0.5 17 HCl. The cations which had 
passed freely through both these columns were adsorbed on a Dowex 50 
column in the acid form and subjected to chromatographic elution with 
1.5 M ammonium lactate pH 7. The strontium fraction was identified by its 
radiation characteristics and by its position in the elution sequence. It was 
found necessary first to purify all reagents by redistilling, or by pouring them 
through a resin column at a pH suitable for the removal of rare earths. Ammo- 
nium lactate was chosen as eluant because of the relative ease of purification. 

The separated strontium was set aside overnight to allow the 9.7-hr Sr® 
to decay to Y", which was subsequently recovered free from strontium and 
rare earths by elution with 0.4 1, pH 4 ammonium isobutyrate from Dowex 
59. From this step onwards all procedures were carried out in polythene, 
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avoiding glass as a possible source of natural rare earths. Yttrium was recovered 
from the ammonium isobutyrate by metathesis to the chloride on a very 
small Dowex 50 column in the acid form. By washing this column carefully 
with water all traces of the isobutyrate were removed before the Y was 
eluted with 6 J HCl. This Y eluate was evaporated to dryness in a platinum 
boat and ignited; the residue was then taken up in 3 WM HCl. The uranium 
content was further reduced by selective adsorption on a 1/10th ml Dowex 1 
column. This was repeated in a glove box to minimize contamination, and 
the final eluant was stored in polythene for several weeks to allow Y® to 
decay. 

(b) Preparation of the Sample for Irradiation 

A preliminary study was made of the fission product contamination resulting 
from the irradiation of empty nylon, quartz, and polythene containers of 
roughly equal size. This study indicated that the activity to be expected as 
contamination in future irradiations would be least if polythene were used, 
the activities produced in the quartz and nylon containers being 3.5 and 6.5 
times as high respectively. From the absolute magnitude of the contamination, 
it was estimated that a cross section of 1 barn could be measured by the 
irradiation of 3.5X10°d.p.m. Y* in polythene. However, the polythene did 
not withstand irradiation satisfactorily and it was decided to use very fine 
quartz tubing, reducing the surface area 10-fold and the contamination 
contribution below that of the polythene. 

Samples were prepared for irradiation as follows: A 10. sample of HCl 
solution containing ~3X10°d.p.m. Y*! was drawn into a warmed capillary 
vial and centrifuged to transfer all the solution to the sealed end. The open 
end was sealed in a flame and the capillary heated for several hours at 110° C 
to check that it was leak-free. The tube, wrapped in superpure aluminum, was 
then put in an iron irradiation capsule along with a cobalt wire flux monitor 


carrier to a high flux position in the NRX reactor for 3 hours. 

(c) Procedures Following Irradiation 

Immediately on delivery from the reactor the quartz sample tube was 
washed in boiling aqua regia, broken into a beaker containing 5 mg of yttrium 
carrier in concentrated HCl, and boiled to facilitate solution and exchange. 
Radioactive contamination was removed from the Y by the following treat- 
ment: passage through a column of Dowex 1 with 10 4 HCl as complexing 
agent, three hydroxide precipitations alternated with three fluoride precipi- 
tations, and elution from Dowex 50 with 0.3 M isobutyrate, pH 4.0. The 


separated yttrium was recovered as oxalate, washed with water, alcohol and 
ether, mounted on an aluminum counting tray, and covered with a mylar 


film. 

In several preliminary runs a search was made for such nuclides as Np’, 
Ba, and Sr, in an effort to set an upper limit on fission contamination 
present (i.e. a measure of the total uranium in the sample plus capsule). 

(d) Preparation of a Source of Yttrium-92 

In order to standardize the counting equipment, and also to estimate the 
relative abundances of the three 6-particles associated with Y® decay, it was 
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necessary to prepare a relatively pure source of this isotope. This was obtained 
from the decay of 2.7-hr Sr® during the preparation of Sr® (see Chemical 
(a)), by milking yttrium from the cation column following a 1-hour growth 
period, after the removal of gross rare earth activities, and before the strontium 
elution. 


Counting 

An anthracene crystal 5cm in diameter and 1cm thick, attached to a 
6342 RCA photomultiplier tube, was used for this measurement. The crystal 
and tube, wrapped with aluminized mylar film, rested in a lucite ring to the 
underside of which a slide holder was attached. This was connected through 
a transistor preamplifier (AEP 1403) to a fast high gain amplifier (AEP 1440, 
the Moody amplifier). An anticoincidence shield was provided by a layer of 
nine 12-in. long, 1-in. diameter Geiger counters directly under the sample 
holder and along two sides of the crystal. Outer shielding consisted of 8 inches 
of iron plate on the sides and underneath, with 2 inches of lead surrounding 
the photomultiplier. 

Two methods of data collection were utilized during the experiment: 

(1) Pulses from the amplifier were fed through a discriminator before 
entering the anticoincidence circuit and scaler. By suitable adjustment of 
the discriminator, counting efficiencies of 7% for Ru'*-Rh!% 8-spectrum 
(end-point energy 3.5 Mev) were obtained at bias settings of about 3 Mev, 
for which the efficiency of counting Y* was 10-*°%. Background readings at 
these settings averaged 0.2 counts per minute. 

(2) Pulses from the Moody amplifier were fed through a 1-ysec delay line 
to a 100-channel pulse height analyzer (AEP 2230) gated by the anticoincident 
circuit. By preparation of a Fermi-Kurie plot to establish the cutoff channels 
for the various 8-energies present, and by comparison with similar data 
collected for a relatively pure source of Y®, it was possible to determine the 
contribution of the 3.6-Mev component and plot its decay. 

A standard Ru!*-Rh' source was used to estimate efficiencies for both 
methods. 


RESULTS 
A. Calibration of the Counter 


The spectrometer was calibrated for these experiments by the preparation 
of a standard Fermi-Kurie plot of the spectrum of a Ru!%—Rh! source. Sub- 


stituting the approximation V[|N(¢)]/for the quantity [N(e)}[f(Z,e)e(@—1)4}-? 
(where N(e) is the total counts per unit time in an energy interval e¢, € is the 
total energy in units of moc?, and f(Z,e) is the fermi function) and using the 
known end-point energy (the 3.53-Mev 8 of Rh!) plus an estimated cutoff, 
it was possible by iteration of this plot to determine the cutoff channel for 
this energy. The error introduced by equating the electron momentum to the 
electron energy and neglecting the charge on the Rh nucleus (i.e. putting 
f(Z,e) = 1) and by using a thin crystal probably totalled <3% in deviations 
from the straight line plot and is negligible in determining the end point. 
The energy/channel value obtained from the Rh! cutoff was then used 
in the preparation of similar plots for Y", and in an estimation of the cutoff 
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for Y%. In the case of the Y* the source used was the sample itself, imme- 
diately prior to irradiation. Figures 1(a@) and 4(@) show the Rh! cutoff at 
the start of the two experiments outlined in Sections C and D; Figs. 1(0) 
and 4(b) show the Y* cutoff on the same dates. Using the allowed shape for 
the Y first-forbidden transition also introduced a negligible error into that 
calibration. Using the end-point energy of 1.54 Mev quoted in the most recent 





ae 
€ 


ARBITRARY 
UNITS 











10 20 30 
CHANNEL NUMBER 


Fic. 1. (a) @ Rh’ spectrum at time of first expt.; cutoff channel 70.8; energy/channel 
= 50 kev. (b) @ Y® spectrum prior to first expt.; cutoff channel 31; energy/channel = 49.7 
kev (estimated cutoff for Y® = 45.7). (c) @ Scan of Y® prior to first expt. with instrumental 
background subtracted; Y* cutoff channel 73; energy/channel = 49.4 kev. (d) First expt.: 
X original scan, Y*" cutoff channel 31, energy/channel = 50 kev; @ original scan X10, 
minus Y* contribution, Y® cutoff channel 46; J original scanX100, minus Y® and Y® 
contributions, Y® cutoff channel 72. 
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table of isotopes (Strominger ef al. 1958) the energy/channel values obtained 
for the Y and Rh plots are seen to be in good agreement. 

The preparation of a sample of Y® for calibration purposes has already 
been described, and Fig. 2 substantiates the purification procedures. The half- 
life of the decay observed in channels 46-60 over a 3-day period was 3.5 
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counts /HR 





Ol 


° 20 40 60 80 100 120 
TIME — WR 


Fic. 2. Y® calibration: X total counts/hour in channels 46-60; @ longer-lived com- 
ponent subtracted. 
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hours, with less than 0.13% longer-lived contamination seen. Three 8-particles 
have been identified in the decay of Y® (Ames, Bunker, and Starner 1952). 
Little information is available on the relative abundances of the 6-rays and 
as seen in Fig. l(c), the Fermi-Kurie plot obtained from the scan of. the 
spectrum showed only one end point above 1.3 Mev. However, Dzelepov and 
Peker (1958) list the 3.6-Mev 8-branch from Y® as 77% per disintegration 
and in the absence of further information it was decided to use this value 
for the purpose of calculation of the cross section. 


B. Determination of the Uranium Content of the Irradiated Material 

As pointed out earlier, any uranium present in the Y® at the time of 
irradiation will undergo fission and contribute to the Y® observed in this 
experiment. Attempts were made to determine this contribution by separating 
and counting Sr*!, Ba, and Np”? formed during the irradiation. Np** gave 
the most unambiguous results, permitting the estimation of as little as 10-" g 
of uranium in the sample plus capsule. The quantity of uranium which would 
produce the observed activity of neptunium (42 dis./hr) would also give rise 
to 90 dis./hr Y® at the time of its separation from the parent Sr%. This 
corresponds to ~2% of the total Y® observed. 


C. First Y* Cross-section Measurement 

A Fermi-Kurie plot prepared from the first 100-channel scan of the source 
(approximately 4 hours from the end of the 2-hour irradiation), using e« values 
determined from the Y* calibration, Fig. 1(0), is shown in Fig. 1(d). Once 
the Y* contribution is subtracted the remainder of the activity is easily 
separable into two components with cutoffs at channels 46 and 72, identifiable 
with reasonable certainty as Y°° and Y® respectively. 

In order to determine the decay of the component with end point in channel 
72, that portion of the total counting rate appearing in channels 50-60 was 
summed for successive scans and plotted as a function of time, as shown in 
Fig. 3. (Channels 60-70 were not included in the totals because counting 
rates in that region were so low as to make unusual background fluctuations 
assume too great an importance.) As can be seen from the graph there was 
an appreciable long-lived component contributing to the counting rate in 
these channels, probably as a result of pile-up of Y* to give an apparent 
3.1-Mev particle, but when this component was subtracted a 3.7-hour activity 
remained. Extrapolating this curve to fo (210 counts/hr) and assuming a 
counting efficiency equal to that for the standard rhodium source (68% of 
the rhodium disintegrations by emission of a 3.53-Mev particle and 77% 
of the yttrium by a 3.6-Mev particle) in these same channels, a value of 
6800 dis./hr Y*% at to) was obtained. 

Substituting in the formula: 


(J ) 091 = (NX) 92/No1 -F 

where F = average flux as determined by production of Co during the same 
period in an accurately weighed Co wire (Jervis 1957), relative to o2200 = 37.0 
barns, 

(NA)oo = dis./hr of Y® at fo, 
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Ny, = total number of atoms of Y* determined by comparison with a 
smaller source of the same sample whose activity was measured in a calibrated 
Geiger counter. 





Then 
ey et 6800 /0.327 ee eda 4b 
791 = 703310 %<6.45xX10"xX3600 ~ 127X10% cm’, 
where F = 6.45 K 10 neutrons/cm? sec, 


Noi = 7.03 X10'* atoms, 


(NX) 92 = 6800/0.327 dis. per hr where the factor 0.327 represents the 
fraction of saturation reached during irradiation. 


ARBITRARY 
UNITS 
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Fic. 4. (a) @ Rh" at time of second expt.; cutoff channel 66.6; energy/channeil = 53 kev. 
(b) @ Y* spectrum prior to second expt.; cutoff channel 27.8; energy/channel = 55.4 kev 
(estimated cutof! for Y® = 41). (c) Second expt.: X original scan, Y® cutoff channel 28, 
energy/channel =: 55 kev; @ original scanX10, minus Y®" contribution, Y® cutoff channel 
42; If original scanX100, minus Y" and Y® contributions, Y® cutoff channel 68. 
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Unfortunately, as a result of the choice of too high a discriminator setting 
the counting rates obtained by the method 1 outlined earlier were so low 
that analysis of the data was impossible. A lower setting, yielding higher 
counting rates, but possibly including some longer-lived low energy com- 
ponents, was planned for the second experiment. 


D. Second Y* Cross-section Measurement 

Figure 4(c) is a Fermi-Kurie plot of the first scan of the source (approxi- 
mately 4.5 hours from the end of the 7-hour irradiation): As shown in Figs. 
4(a) and (0b) a shift in gain of the amplifier between experiments had been 
sufficient to make an appreciable change in the cutoff channels for the energies 
checked. Since the Y* cutoff had moved down to channel 42 it was decided 
that channels 46-60 could safely be summed for the hard-component decay 
curve, thus partially compensating for the decreased counting efficiency 
measured in channels 50-60. Figure 5 shows the decay curve so obtained, and 
the 3.7-hour activity remaining after the subtraction of a 40-hour activity. 
This longer-lived component probably results from Y°* pile-up, plus Y* 
(1.54 Mev), and Y°* (2.27 Mev) adding to give an apparent 3.8-Mev @-particle. 
(Since the sample of Y* used was considerably smaller than that of the 
previous experiment the irradiation time had been increased by a factor of 
3.5; hence the Y° contamination was much larger in this instance.) By extra- 
polating the 3.7-hour component of the activity to fo (270 counts/hr) as in 
the previous case, and using the Rh efficiency obtained from this experiment, 
a value of 9200 dis./hr Y® at to was obtained. 

Figure 6 shows the decay observed using a discriminator as outlined pre- 











Fic. 6. Second expt.: X decay at discriminator setting yielding 2.35% efficiency for 3.6 
Mev; @ 64-hour component subtracted. 
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viously. Subtraction of a 64-hour background (Y%) left a component of 
~4-hour half-life, which may be inferred from the data to be Y%. Extra- 
polation of this component activity to tf) (198 counts/hr), and use of the Rh 
efficiency measured for this discriminator setting has led to a value of 11000 
dis./hr Y* at to. This value is in reasonable agreement with that obtained 
by method 2, their average being 10,100 dis./hr. 

Substituting this in equation (1), yields 


fe _10,100/0.74 ns al a ar 
72 = 47x10" x6.12x10"x3600 ~ 11X10 em 
where F = 6.1210" neutrons/cm? sec, 


Noi = 4.1X10!° atoms, 
(NA)o2 = 10,100/0.74. 


Hence, for the two experiments the average value of the cross section of 
Y® = 1.4+0.3 barns. 
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NOTES 
A NOTE ON HANSEN’S VECTOR WAVE FUNCTIONS 


T. B. A. SENIOR 


Almost a quarter of a century ago, Hansen (1935) proposed a method for 
generating solutions of the vector wave equation in any co-ordinate system 
in which the scalar wave equation is separable. The solutions are, essentially, 
a development of those employed by Mie (1998) in a study of scattering by a 
sphere, and were used by Hansen (1935, 1936, 1937) to analyze various prob- 
lems dealing with the radiation from antennas. Because of these applications, 
Hansen was primarily concerned with spherical polar co-ordinates, and for 
the extension of the method to more general co-ordinate systems the credit is, 
perhaps, due to Stratton (1941). 

The three vectors which Hansen proposed are 


(1) L =Vy 
(2) M =V A (oy) 


where y is a soliition of the scalar wave equation 
(4) Vy+ ky =0 


in the appropi ate co-ordinate system, and ¢ is either a radial vector or, more 
generally, is such that M and N satisfy the vector wave equation 


(5) VAVAA = RPA. 


It will be observed that L has zero curl, and accordingly L is suitable for the 
representation of an electrostatic or magnetostatic field. On the other hand, 
V.L = —ky, and therefore L cannot appear in the representation of a sole- 
noidal function except in the limit k = 0, whereas M and N can. Indeed, from 
equations (3) and (5), it follows that 


(6) M=<=VAN 


and since M and N are now related by curl operations, they are ideally suited 
to the representation of an electromagnetic field (E, H). 

The above vectors have formed the basis for most attempts at the solution 
of vector scattering problems in electromagnetic theory. The ease with which 
a solution can be obtained in any given case is vitally effected by the choice 
of the vector 6 in terms of which M and N are defined, and it is conceivable 
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that in some instances an unwise choice of 6 could prevent the determination 
of a solution. When the co-ordinate system involved is more complicated than 
spherical, it therefore pays to devote some attention to the selection of 6 
before attempting the solution, and it is now of interest to know the possibilities 
which are available. Fortunately it is not difficult to find the most general 


form for 6. 
For this purpose it is only necessary to consider the vector M. From equation 


(2), 


(7) M = LAé+y¥VAz4, 

giving 

(8) VAM = Fey + P, 

where 

(9) P =V(6.L) — 2.L.V)é+LV.6 +¥VAVAG, 


and hence 


VAVAM = PM +V AP. 


It now follows that if M is to satisfy the vector wave equation,V AP must be 
zero, and since it is not identically zero for all 6, this condition, coupled with 
the fact that é6 must be independent* of y, serves to determine the most 
general 6. 

Since the vector P is defined in terms of vector operators, it is sufficient to 
carry out the analysis using a particular set of co-ordinates and then use the 
invariance of these operators with respect to transformations of co-ordinates 
to give the specification of ¢ in general. For simplicity, Cartesian co-ordinates 
(x1, X2, X3) are therefore chosen, and the independence of é and y is employed 
by choosing in turn three simple forms for y. 

If ~ = (x), the condition V A P = 0 leads to six differential equations 
involving the rectangular components a; of 6. Similarly if y = ¥(x2) or y = 
¥(x3), and the resulting set of 18 differential equations is compatible if and 


only if 
01 = o;(%}), etc. 
and 
V?6 = 0. 
Hence 


a1 = aX; + Bi, 
2 = aXe + Ba, 

o3 = a3xX3 + B3, 
where the a; and 6; «re constants, and each component of ¢ is now a first 
order linear function of the corresponding variable. Substituting into equation 
(9), we have 


*Were 6 not independent of y, the vectors M and N would have completely different forms 
for incident and scattered fields, and this would destroy much of their usefulness. 
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(10) P=V(6.L)+)) a(L—2L%,) 


and since this is only irrotational if a1 = a2, = a3(= a, say), 
(11) é = (1X1 + XoXo +- 3X3) + BiX1 + BoXe + B3X3. 


This represents the necessary and sufficient condition that M (and therefore 
N) satisfy the vector wave equation in Cartesian co-ordinates, and by using 
the invariance under transformation of co-ordinates it now follows that the 
most general form for 6 in any co-ordinate system consists of a radial vector R, 
together with a set of constant (unit) vectors. It has been known for a long 
time that both radial and constant vectors are possible, but the above proof 
is valuable in showing that no other forms for ¢6 are permissible. It will be 
observed that only with Cartesian co-ordinates are the unit vectors in the 
direction of the co-ordinate axes. 

If equation (11) is interpreted in terms of a general co-ordinate system, 
é is such that 


VAé=0, 

V.¢é = 3a 
and 

(LV )é = aL, 
so that 

P =V(ay + 6.L). 

Hence 
(12) M =LA<ge, 
(13) kN = Rey +V (ay + 6.L), 


and the fact that the second group of terms on the right-hand side of (13) can 
be written as a single gradient is necessary if M is to satisfy equation (5). 

In any practical application of Hansen’s vectors it is unlikely that ¢ will be 
required in its most general form (11), and it is usual to take one of a, Bi, 
62, and 83 unity with all the others zero. One of the most important cases is 
obtained by taking a = 1 so that 6 = R. This is the standard choice for the 
scattering of an electromagnetic wave by a sphere (see, for example, Stratton 
1941), and also proves convenient in the corresponding problem for an ellipsoid 
(Senior 1961). 

A final point concerns the behavior of the vectors in the static limit. If 
vy — gas k— 0, where V*¢ = 0, then 


L-S 
where S is the electrostatic (or magnetostatic) field defined as S =V¢. Also 


M—SAé 


and 


kN — o&S +V (6.8), 





1 an aan nals Ae ani eae hatte aN Lhe sls eee amare iae 


AeA a aS NR cin NS eat NGL SHY 8. 








3 
4 
4 
3 
7 
: 
3 
\ 
¥ 
4 
4 
q 
5 
q 
{ 
4 


i ran i 


a ANN 8 


NOTES 1705 


so that in the limit of zero frequency the field cannot be represented in terms 
of the vectors N alone. This has important consequences in any attempt to 
determine the low frequency (or Rayleigh) expansion for the field, and implies 
that in a general co-ordinate system both M and N vectors are necessary to 
match the static field. On the other hand, if spherical polar co-ordinates are 
employed with 6 = R, 


0 
y+é.L= ar (RY) 
and since y is a function of R only in the combination RR, 


= 
OR 
as k > 0. In this co-ordinate system therefore, 


kN — const. XS 


(Ry) — const. X @ 


and the static field can be represented entirely by the vectors N. 
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THE FRANCK-CONDON FACTOR (q,.) ARRAY TO HIGH VIBRATIONAL 
QUANTUM NUMBERS FOR THE 0.(B*=z — X*z}) SCHUMANN-RUNGE 
BAND SYSTEM 


R. W. NIcHOLLS* 


The Schumann—Runge band system of Oz is unusual in view of the relatively 
large (~0.4 A) difference between the equilibrium internuclear separations (7,) 
of the two electronic states concerned. The Franck—Condon principle thus 
implies that for this system the strongest bands lie on a very wide quasi 
parabola in the v’—v’’ plane. The limbs of the curve lie mainly at large values 
of v’ and v” along the v’’ = 0 and v’ = 0 progressions respectively, and the 

*Temporarily on leave of absence from the University of Western Ontario. 
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curve avoids the (0,0) region. This primary ‘‘Condon parabola” and the 
subsidiary parabolae are indicated in Table II. 

The accurate computation of arrays of Franck—Condon factors (vibrational 
overlap integral squares (|fy»y, dr|?)) for such transitions which involve 
relatively large Ar, and thus high values of v’ and v” presents two serious 
problems. The first, which is very common in quantum mechanics, is the 
severe cancellation between the often almost equal positive and negative 
parts of the integrand. The second is the inability of the empirical analytic 
models of molecular potentials, such as the Morse potential, to represent true 
molecular potentials particularly at the regions of very high vibrational quan- 
tum number. The form of the wave functions involved in the overlap integral 
is, of course, very much determined by the potential used. 

The first of these problems has been overcome in some cases by direct 
computation (using desk calculators) carrying as many significant figures as is 
possible (Jarmain and Nicholls 1954) and in others vy the development of 
approximate methods (which are most accurate at low values of v’ + v’’) 
in which the remainder left after cancellation is evaluated (Fraser and Jarmain 
1953; Jarmain and Fraser 1953; Fraser 1954). A large number of important 
band systems have been treated by such methods (Jarmain, Fraser, and 
Nicholls 1953, 1955; Fraser, Jarmain, and Nicholls 1954; Nicholls, Fraser, 
and Jarmain 1959; Nicholls, Fraser, Jarmain, and McEachran 1960). The 
second problem is harder to overcome as all of the current methods for evaluat- 
ing Franck-Condon factors require analytic wave functions specific to some 
form of potential. Recent numerical studies of the form of molecular potentials 
using the Klein-Dunham method (Jarmain 1959, 1960) or the very similar 
Rydberg—Klein—Rees method (Vanderslice, Mason, Maisch, and Lippincott 
1959; Vanderslice, Mason, and Lippincott 1959; Vanderslice, Mason, and 
Maisch 1959, 1960) to delineate the profile of vibrational oscillator turning 
points show that the Morse potential is a very fair representation of molecular 
potentials at low quantum numbers in many cases, but that each case has to 
be considered on its own merits. 

Ultimately, of course, it will be necessary to obtain numerically the wave 
functions appropriate to the numerical potentials for computation of exact 
Franck—Condon factors, and work on this is in progress. 

Arrays of Franck—Condon factors for the O: Schumann—Runge system have 
been computed with various degrees of accuracy and completeness by a 
number of workers over the past few years using a variety of approximate 
methods. This work is summarized in Table I. Each of these calculations is 
subject to some inaccuracy arising from limitations of the method used, from 
inappropriateness of the potential employed, and incompleteness in the 
number of bands studied. 

The facilities of the N.B.S. Computation Laboratory including an IBM-704 
computer have therefore been used to compute a complete array (0< v’ < 21; 
O0<v’’ < 21) of Franck—Condon factors for all the possible bands using a 
Morse model for potentials in each case. The straightforward method used 
has previously been used by Jarmain and Nicholls (1954) for the N>» first and 
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second positive systems. The zeroth vibrational wave function ¥> was computed 
at 0.01 A intervals for each of the two potentials using the analytic form of 
this function given by Morse (1929). A recursion relation (Nicholls, Jarmain, 
and Fraser 1953; Fraser, Jarmain, and Henderson 1959) was then used to 
build up, step by step, numerical wave functions at 0.01 A intervals for all 
known higher vibrational levels. Overlap integrals of all pairs of wave func- 
tions between the two states were then evaluated numerically using the 
trapezoidal and Simpson rules. Agreement between the results of these two 
methods of integration was excellent, which showed that the interval of 
integration was sufficiently small. The results of this computation are shown 
in Table II. 

It is of course realized that the Morse potential is not a particularly realistic 
representation for the B*Z_z state of O2 specifically at high vibrational quantum 
numbers, in view of the relatively large value of w,y, (Brix and Herzberg 
1952), some evidence of predissociation in it (Carroll 1959), and a comparison 
between the Morse and “‘true’”’ numerical potential (Vanderslice, Mason, and 
Maisch 1960; Jarmain 1960). Nevertheless when used in conjunction with 
intensity measurements in emission (for which the v’ = 0, 1, 2, 3 progressions 
are the most important) the data of Table II are reliable and realistic. They 
have been recently used by Hébert (1960) for the interpretation of intensity 
measurements in emission for this system. Full details of this work will be 
published shortly. 

These calculations represent the first of a series for which the IBM-704 
program is being used to evaluate arrays of Franck-Condon factors and 
r-centroids for important molecular transitions. Vibrational wave functions 
calculated during this work are on file on magnetic tape. 


The author’s sincere thanks are due to Miss I. Stegun and Miss R. Zucker 
of the N.B.S. Computation Laboratory for writing and testing the program 
and to Mr. W. R. Jarmain for providing details of the method of computation 
in a form suitable for programming. 


Brix, P. and HERZBERG, G. 1954. Can.-J. Phys. 32, 110. 

CaRROLL, P. K. 1959. Astrophys. J. 129, 794. 

Fraser, P. A. 1984. Proc. Phys. Soc. A, 67, 939. 

Fraser, P. A. and JARMAIN, W. R. 1953. Proc. Phys. Soc. A, 66, 1145. 

Fraser, P. A., JARMAIN, W. R., and NICHOLLS, R. W. 1954. Astrophys. j. 119, 286. 

FRASER, P. A., JARMAIN, W. R., and HENDERSON, G. 1959. Scientific Report No. 2, Contract 
AF 19(604)-4560, University of Western Ontario. 

HésBert, G. R. 1960. Ph.D. Thesis, University of Western Ontario. 

JarMAIN, W. R. 1959. J. Chem. Phys. 31, 1137. 

—— 1960. Can. J. Phys. 38, 217. 

Jarmatn, W. R. and Fraser, P. A. 1953. Proc. Phys. Soc. A, 66, 1153. 

JARMAIN, W. R. and NicHo.tis, R. W. 1954. Can. J. Phys. 32, 201. 

JARMAIN, W. R., FRASER, P. A., and NICHOLLS, R. W. 1953. Astrophys. J. 118, 228. 

1955. Astrophys. J. 122, 55. 





KEck, J. C., Camm, J. C., KIveL, B., and WENTINCK, T., Jk. 1959. Ann. Phys. 7, 1. 

Morse, P. M. 1929. Phys. Rev. 34, 57. 

NICHOLLS, R. W., JARMAIN, W. R., and Fraser, P. A. 1953. Can. J. Phys. 31, 1019. 

NICHOLLS, R. W., FRASER, P. A., and JARMAIN, W. R. 1959. Combustion and Flame, 3, 13. 

NICHOLLS, R. W., FRASER, P. A., JARMAIN, W. R., and McEacuran, R. P. 1960. Astrophys. 
J. 131, 399. 








snc, 


NOTES 1711 


Pittow, M. E. 1950. Proc. Phys. Soc. A, 63, 940. 

1952. Proc. Phys. Soc. A, 65, 858. 

VANDERSLICE, J. T., MASON, E. A., Matscu, W. G., and Lippincott, E. R. 1959. J. Mol. 
Spectroscopy, 3, 17. 

VANDERSLICE, J. T., Mason, E. A., and Lippincott, E. R. 1959. J. Chem. Phys. 30, 129. 

VANDERSLICE, J. T., MASON, E. A., and Maiscn, W.G. 1959. J. Chem. Phys. 31, 738. 

1960. J. Chem. Phys. 32, 515. 








RECEIVED AuGusT 12, 1960. 
NATIONAL BUREAU OF STANDARDS, 
WASHINGTON 25, D.C. 








ii Aorta 





LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not laier than six weeks | 
previous fo the first day of the month of issue. No proof will be sent to the authors. \ 


Frequency Measurement of Standard Frequency Transmissions" ? 


Measurements are made at Ottawa, Canada, using N.R.C. caesium-beam frequency 
resonator as reference standard (with an assumed frequency of 9 192 631 770 c.p.s.). Frequency 
deviations from nominal are quoted in parts per 10". A negative sign indicates that the 
frequency is below nominal. 


GBR, 16 ke/s 











Date, —_—_—_—OOOoO a 
September MSF, 5-hour 24-hour WWVB, 
1960 60 kc/s average* average 60 kc/s 

1 —159 —158 N.M. —150 

2 —161 N.M. N.M. —149 

3 — 159 — 158 N.M. N.M. 

4 N.M. N.M. N.M. N.M. 

5 N.M. N.M. N.M. N.M. 

6 — 160 N.M. N.M. —148 

7 — 150 —154 —155 — 148 

8 — 157 —154 — 156 — 146 

y —154 —151 —155 N.M. 

10 — 156 — 155 — 155 N.M. 

11 — 159 — 157 — 157 N.M. 

12 —154 —158 —151 —144 

13 — 153 N.M. N.M. —144 

14 — 162 — 146 — 144 — 146 

15 — 147 —148 —146 — 143 

16 — 148 — 153 —149 —144 

17 —150 — 148 —151 N.M. 

18 —161 —i54 — 154 N.M. 

19 —147 N.M. N.M. — 140 

20 — 146 —151 —150 —141 

21 —148 —151 —151 —140 

22 —145 — 147 — 147 —139 

23 — 148 —154 —150 —140 

24 N.M. —152 N.M. N.M. 

25 N.M. N.M. N.M. N.M. 

26 —149 N.M. N.M. — 138 

27 —145 —149 —149 —137 

28 — 148 —149 N.M. — 137 

29 — 146 — 142 N.M. — 136 

30 —145 — 146 — 148 — 137 

Midmonthly 
mean — 153 — 152 — 152 —142 
Midmonthly 


mean of WWV — 142 


Note: N.M. no measurement. 
*Time of observations: 0.100 to 04.30 and 11.30 to 13.00 U.T. 


RECEIVED OCTOBER 17, 1960. S. N. KALRA 
Division OF APPLIED Puysics, 

NATIONAL RESEARCH COUNCIL, 

OTTAWA, CANADA. 


1Issued as N.R.C. No. 6021. 
2Cf. Kalra, S. N. 1959. Can. J. Phys. 37, 1328. 
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Spin Rate of the Satellite Echo I as Determined by a Tracking Radar 


Echo I has been observed on nine occasions at the Prince Albert Radar Laboratory,} 
(PARL), since its launch on August 12, 1960. The radar signal reflected from the satellite is 
detected in two orthogonal channels, right and left circular polarization. Transmitted polari- 
zation is left circular. Since reflection of a circularly polarized radio signal changes the direction 
of rotation of the polarization, the main component of the received signal has right circular 
polarization. However, because of the finite isolation between the two channels, and some 
depolarization of the signal upon reflection, power is detected in the left circular channel 
as well. 

On all passes tracked at PARL (Table I) a regular fast scintillation has been observed in 
the received signal reflected from the satellite. The scintillation rate remained const:t at 
134 fades per minute to pass number 131. On the next pass observed, number 217, the scin- 
tillation rate had changed to 154 fades per minute. During the early revolutions of Echo I 
the fades were small in amplitude (0-3 db). On pass 131, fades as great as 10 db were observed 
in both channels. However, on pass 217 the fading in the right circular channel had decreased 
to about 5 or 6 db. No fast fading was observed in the left circular channel. Instead slow 
deep fading with a period of about 30 seconds was observed. This period of fading is thought 
to correspond to the period of rotation of the satellite. 








TABLE I 
Scintillation rate Rotation rate 
Date Revolution No. (per minute) (r.p.m.) 
August 12 1 134 1.63 

bs 12 2 = a 

13 v ay 

13 10 

13 11 

13 12 
Ps 16 46 
- 23 131 , 
P= 30 217 154 1.88 


The rotational period, or spin rate, of the satellite may be determined by studying the 
scintillations of the received signal if the method of fabricating the balloon is taken into 
account. The balloon is made up of 82 gores. If it is assumed that 

(a) the satellite is rotating about an axis which permits the sequential observation of all 
82 gores each rotation, 

(6) the scintillations in the received radar signal are caused by the motion of the gores, 

(c) there is a one-to-one correspondence between each observed fade and the opposition 
of a gore, 
then the spin rate of the satellite can be calculated. Scintillation rates of 134 fades per minute 
(pass 131) and 154 fades per minute (pass 217) correspond to periods of rotation of 36.8 
seconds and 31.9 seconds respectively. The agreement between the calculated period of rotation 
for pass 217 and the observed period of deep fading in the right circular channel indicates 
that the 2..umptions made are valid. 

Radar «b:ervations of Echo | indicate that between passes 131 and 217 the spin rate 
increased from 1.68 r.p.m. to 1.88 r.p.m. It is interesting to note that between these two 
observations the satellite first started passing into the earth’s shadow. This may have caused 
the internal pressure of the gases inside the balloon to decrease, resulting in some shrinkage 
and(or) deformation of the balloon. If it is assumed that ar~ lar momentum was conserved 
between the two observations, a decrease in radius of 7% required to account for the 


increased spin. 
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Simultaneous Observations of Pulsations in the Geomagnetic Field 
and in Ionospheric Absorption 


During the latter part of September 1959, a cosmic noise absorption monitor (sometimes 
called a riometer) and a magnetometer were operated at Saskatoon with the recording meters 

running at a sufficiently high chart speed (12 inches/hour) to permit fluctuations of a few 
seconds’ duration to be resolved. The riometer, operating on a frequency of 34 Mc/s, is broadly 
similar to that described by Reid and Collins (1959) and the magnetometer is of the standard 
three-component fluxgate type described by Serson (1957) but only one (the H) component 
was recorded at high speed. 

On the 20th and 21st of September, shortly after 7 a.m. local time, regular pulsations were 
recorded simultaneously on both of the records. The magnetic pulsations of amplitude 30-50 
gamma, together with the simultaneously recorded pulsations in cosmic noise absorption, 
of amplitude 2.5-3 db, are shown in Figs. 1 and 2. The riometer record indicates the amount 
of absorption suffered by extraterrestrial radio waves as they traverse the ionosphere while 
the magnetometer record indicates the variations of the earth’s magnetic field due, presumably, 
to currents located either in the lower ionosphere or at great distances from the earth. Cosmic 
noise absorption is usually attributed to abnormal increases in the electron density, collision 
frequency, or both at heights below about 95 km in the atmosphere. 

It can be seen from the figures that for many cycles the geomagnetic pulsations remain in 
phase with the cosmic noise absorption fluctuations but there are occasions when the two 
phenomena appear to be exactly 180° out of phase. This occasional sudden change in phase 
would appear to rule out the possibility that the geomagnetic pulsations are a purely local 
phenomenon brought about by a local periodic increase of the conductivity of the lower 
ionosphere. 

Jacobs and Obayashi (1958) have found experimental evidence in support of the suggestion 
that some magnetic pulsations are due to a toroidal hydromagnetic oscillation of the earth’s 
outer atmosphere of the kind proposed by Dungey (1954). According to this picture the 
period of oscillation should change with geomagnetic latitude and for the geomagnetic latitude 
of Saskatoon (60.5° N.) should be about 2.5 minutes. The observed periods, ranging from 
2.5 to 3.0 minutes, are in good agreement with the predicted value and suggest that the 
pulsations are, in fact, due to toroidal hydromagnetic oscillations of the outer atmosphere. 

Akasofu (1959) has concluded on theoretical grounds that the energy dissipated in the 
ionosphere by hydromagnetic waves would cause only a negligible rise in temperature, so 
that it is necessary to look for some indirect connection between these two simultaneously 
occurring phenomena. One possibility that comes to mind is that extra ionization in the lower 
ionosphere (observed by the riometer) is being caused by the periodic influx of particles. These 
particles could, ir turn, be accelerated during some part of each cycle of the hydromagnetic 
oscillation of the outer atmosphere. 





The author is grateful to Dr. P. A. Forsyth for helpful suggestions and guidance in the 
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